matlab有关地球投影问题clearall;closeall;a=6378245;b=6356863;e=sqrt((a^2-b^2)/a^2);%y=(sin(x*pi/180)).^2;%p=(1-e^2*y).^1.5;%M=(a*(1-e^2))./p;%N=a./(1-e^2*y).^0.5;F=@(x,y)a.^2.*(1-e.^2).*(cos(x)/(1-e.^2.*(sin(x).^2)).^2);%T=[];%forc=1:5;%T(c)=dblquad(F,0,c,0,1);%endT=dblquad(F,0,2,0,1);