#include
#include
#define a 6378245
#define b 6356863
void main()
{
int M,N,r,x;
double y,z,v,w;
double f;
printf(" 子午圈曲率半径 卯酉圈曲率半径 纬线圈半径\n");
f=sqrt((a*a-b*b)/a*a);
y=f*f;
for(x=0;x<=60;x++)
{
z=sin(x)*sin(x);
v=sqrt(1-y*z);
w=(1-y*z);
M=a*(1-y)/(w*v);
N=a/v;
r=a*cos(x)/v;
printf("纬度=%d, %d, %d, %d\n",x,M,N,r);
}
getchar();
}