l0=S0+L-l0;

double

x1=l0-pow(l0,5)/pow(c,2)/40.0+pow(l0,9)/pow(c,4)/3450.0;

double

y1=pow(l0,3)/c/6.0-pow(l0,7)/pow(c,3)/336.0+pow(l0,11)/pow(c,5)/42240.0;

X[i]=(x1-T)*cos(cita)-y1*sin(cita)+T;

Y[i]=(T-x1)*sin(cita)-y1*cos(cita);

}

}

f2=fopen("Result.txt","w"); //输出放样点

的坐标值

if (f2==NULL)

{ printf("cannot open this file\n");

exit(0);

}

for (i=0;i<n;i++)

{

fprintf(f2,"%.3lf\t%.3lf\n",X[i],Y[i]);

}