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
|