q=2*T-L;

// printf("\n%lf\n%lf\n%lf\n%lf\n%lf\n",T,L,E,q,Ly);

S0=402.377;

T=46.763;

Ly=50.356;

c=2000;

m=9.997;

p=0.167;

double

dk[13]={410,420,422.377,430,440,447.555,450,460,470,472.733,480,490,492.733};

int n;

n=13;

double *X=new double[n];

double *Y=new double[n];

double *B=new double[n];

double *s=new double[n];

double cita=dms_arc(1800000)-arfa; //坐标变换时