Right hand method***********
>       r * (o * v - p * e * M),
>       -r * (e * sin(i)),
>       r * (p * v + o * e * M)



Left (POV) method:

M = planets[j].DailyMotion * days;
M += planets[j].MeanLongitude - planets[j].Perihelion;
e = planets[j].Eccentricity;
E = SolveKepler(M, e);
r = planets[j].MeanDistance * (1.0 - e * cos(E));
v = 2.0L * atan(sqrt((1.0L + e) / (1.0L - e)) * tan(E / 2.0L));
o = planets[j].AscendingNode;
p = planets[j].Perihelion;
i = planets[j].Inclination;

r * (cos(o) * cos(v + p - o) - sin(o) * sin(v + p - o) * cos(i)),
r * (sin(v + p - o) * sin(i)),
-r * (sin(o) * cos(v + p - o) + cos(o) * sin(v + p - o) * cos(i))


computers and me dont get along today

Post a reply to this message