代写代做R、R代写、R语言代写代做、代做R程序设计

日期：2018-06-17 03:34

a=6378245

b=6356863

e=(a^2-b^2)/a^2

N1=a/sqrt(1-e^2*(sind(25))^2)

r1=N1*cosd(25)

N2=a/sqrt(1-e^2*(sind(47))^2)

r2=N2*cosd(47)

s1=asind(e*sind(25))/2

U1=tand(45+25/2)/(tand(45+s1))^e

s2=asind(e*sind(47))/2

U2=tand(45+47/2)/(tand(45+s2))^e

c=(log10(r2)-log10(r1))/(log10(U1)-log10(U2))

K=r1*(U1)^c/c

s3=asind(e*sind(15))/2

Us=tand(45+15/2)/(tand(45+s3))^e

Ps=K/(Us)^c

Q1=[15:5:55]';

Q2=70:5:145;

V=c*Q2;

x=(Ps-K./(tand(45+Q1./2)./(tand(45+asind(e*sind(Q1))./2)).^e).^c*cosd(V))/10000000;

y=K./(tand(45+Q1./2)./(tand(45+asind(e*sind(Q1))./2)).^e).^c*sind(V)/10000000;

z=0*Q1*Q2/10000000;

mesh(x,y,z)

view([0,90])