On 9/16/2018 3:00 PM, Mike Horvath wrote:
> https://ssd.jpl.nasa.gov/txt/aprx_pos_planets.pdf
> 
> I am a little confused about the "Solution of Kepler's Equation" on Page 
> 2 of this document. Have I implemented it correctly do you think?
> 
>      #local Orrery_Temp_EStar                        =
180/pi * 
> Orrery_Temp_Eccentricity;
>      #local Orrery_Temp_EccentricAnomaly                = 
> Orrery_Temp_MeanAnomaly + Orrery_Temp_EStar * 
> sind(Orrery_Temp_MeanAnomaly);
>      #local Orrery_Temp_EccentricAnomalyTolerance    = 10e-6;
>      #local Orrery_Temp_EccentricAnomalyDelta        = 0;
>      #while (abs(Orrery_Temp_EccentricAnomalyDelta) > 
> Orrery_Temp_EccentricAnomalyTolerance)
>          #local Orrery_Temp_MeanAnomalyDelta                = 
> Orrery_Temp_MeanAnomaly - (Orrery_Temp_EccentricAnomaly - 
> Orrery_Temp_EStar * sind(Orrery_Temp_EccentricAnomaly));
>          #local Orrery_Temp_EccentricAnomalyDelta        = 
> Orrery_Temp_MeanAnomalyDelta/(1 - Orrery_Temp_Eccentricity * 
> cosd(Orrery_Temp_EccentricAnomaly));
>          #local Orrery_Temp_EccentricAnomaly                = 
> Orrery_Temp_EccentricAnomaly + Orrery_Temp_EccentricAnomalyDelta;
>      #end
> 
> 
> 
> Mike


I did end up finding a bug in that code, but I also discovered that 
there's not a visible difference with the code working or not working. 
So the major problems I am having lie elsewhere.

:(

Mike

Post a reply to this message