用matlab计算M，N，R，r。

在matlab的命令窗口中计算子午圈曲率半径M，卯酉圈曲率半径N，平均曲率半径R，纬圈的半径r。椭球选用的为克拉索夫斯基椭球。角度φ选择为0，π/18，π/9，π/6，2π/9，5π/18，π/3，7π/18，4π/9，π/2，将所需参数以及公式输入命令窗口。

>> a=6378245;

>> b=6356863;
>> e=((a^2-b^2)/a^2)^0.5;
>> phi=0:pi/18:pi/2;
>> N=a./power((1-power(e,2)*power(sin(phi),2)),0.5);
>> M=(a.*(1-power(e,2)))./power((1-power(e,2).*power(sin(phi),2)),1.5);
>> R=sqrt(M.*N);
>> r=N.*cos(phi);
>> log10(M);
>> log10(N);
>> log10(R);
>> log10(r);
>> M

tip：关于.*和*，./和/的问题，前者是矩阵中对应元素相乘，后者是两个矩阵相乘，要求维数一致。