#include
#include
double?fun1（double?adouble?bdouble?c）
{
	double?Ugh;
	h=a/2;g=（1-a*sin（b/180*c））/（1+a*sin（b/180*c））;
????U=tan（c/4+（b/180*c）/2）*pow（gh）;
????return?U;
}
double?fun2（double?adouble?bdouble?cdouble?e）
{
	double?r;
	r=（a*cos（b/180*c））/（sqrt（1-e*e*sin（b/180*c）*sin（b/180*c）））;
	return?r;
}
int?main（）
{
????double?Ur0e=0.081819191a=6378137xyrupt=3.1415926glx1y1;/*t=πu=μ*/
????double?BB0=30LL0=120;
????for