一定经纬范围的墨卡托投影.cpp

#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