gldhdryden
(Gordon Dryden)
February 7, 2017, 9:58pm
1
I get this kind of output:
In[1]:= ae = 6378137.0
eps2 = 0.00669437999014
eps = Sqrt[eps2]
phivec = 45
lam = 0
h = 1000.
phi = phivec*pi/180 // N
Q = ae/Sqrt[1 - eps2* Sin[phi]^2] // N
Out[1]= 6.37814*10^6
Out[2]= 0.00669438
Out[3]= 0.0818192
Out[4]= 45
Out[5]= 0
Out[6]= 1000.
Out[7]= 0.25 pi
Out[8]= 6.37814*10^6/Sqrt[1. - 0.00669438 Sin[0.25 pi]^2]
mthiel
(Marco Thiel)
February 7, 2017, 10:11pm
2
Dear Gordon,
I am not really sure what your question is and your code is kind of difficult to decipher. Also I need to make assumption such as that your pi is actually Pi. But if I do that something like this comes out:
ae = 6378137.0 ;
eps2 = 0.00669437999014;
eps = Sqrt[eps2];
phivec = 45;
lam = 0;
h = 1000.;
phi = phivec*Pi/180 // N;
Q = ae/Sqrt[1 - eps2*Sin[phi]^2] // N
That gives the “single number”: 6.38884*10^6.
Cheers,
Marco