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

Post a reply to this message