У параболы для расчета истинной аномалии v и радиус-вектора r используют среднюю аномалию M и расст. перицентра q ( а не большую полуось a, которая у параболы не имеет смысла).
И второе - ср. аномалию нельзя приводить к интервалу 0 .. 360 гр ( 0 .. 2пи) . Если она получилась при расчете 400, так и использовать, а не 400-360 = 40 , как можно бы было у эллипса.
Итераций, как в формуле Кеплера не требуется, решается через корень кубического уравнения.
Вот как я делал. Правда, на практике редко пользовался, на практике эксцентриситет либо чуть больше 1, либо чуть меньше. Разве что в качестве приближения.
procedure keTrueAnomRadVectParab(q,M:Extended;var v,r:Extended); // M, v - радианы
var S,x,si,co,S2:Extended;
begin
x:=Power(12*M+4*Sqrt(9*M*M+4),1/3)/2;
S:=x-1/x;
S2:=S*S;
co:=(1-S2)/(1+S2);
si:=2*S/(1+S2);
v:=ArcTan2(si,co);
r:=2*q/(1+Cos(v));
end;
Расчет среднего движения n для параболы
function QtoMeanMotionParabol(q,k:Extended):Extended;
begin
Result:=k/(Sqrt(2*q)*q);
end;
k постоянная Гаусса. Если единицы - а.е., масса Солнца, сутки, то k = 0,017209895 рад/сутки.
Если в системе СИ, то k = корень из GM Солнца.
Средняя аномалия = среднее движение * время от перигелия.
