exit(0);

}

int DK0;

double S0,arfa,R,lh,dl;

fscanf(f1,"%d%lf%lf%lf%lf%lf",&DK0,&S0,&arfa,&R,&lh,

&dl);

arfa=dms_arc(arfa);

fclose(f1);

double m,p,c,Bh; //计算m,p,c,Bh

c=R*lh;

m=lh*0.5-lh*lh*lh/R/R/240.0;

p=lh*lh/R/24.0;

Bh=0.5*lh/R;

// printf("%lf\n%lf\n%lf\n%lf\n",c,m,p,Bh);

double T,Ly,L,E,q; //计算T,Ly,L,E,q

T=m+(R+p)*tan(arfa/2.0);

S0=S0-T;

Ly=(arfa-2*Bh)*R;

L=(arfa-2*Bh)*R+2*lh;

E=(p+R)/cos(arfa/2.0)-R;