Commit a124f8db authored by Alexander Lapshin's avatar Alexander Lapshin

.

parent 50524635
...@@ -285,18 +285,20 @@ void ODEPhaseSystem::EarthNonCentralAcc(Vect3 &acc, const Vect3 &pos) ...@@ -285,18 +285,20 @@ void ODEPhaseSystem::EarthNonCentralAcc(Vect3 &acc, const Vect3 &pos)
x12 += t1; x12 += t1;
x01 += t2 * (n + 1); x01 += t2 * (n + 1);
x02 += t1 * (n + 1); x02 += t1 * (n + 1);
if (n > m + 2)
if (n > m + 2)
{ {
Pm1[n] = ((2 * n - 1) * zr * Pm1[n - 1] - (n + m) * Pm1[n - 2]) * (*pn); Pm1[n] = ((2 * n - 1) * zr * Pm1[n - 1] - (n + m) * Pm1[n - 2]) * (*pn);
pn++; pn++;
}; }
if (n != m) if (n != m)
{ {
t1 = Pm1[n]; t1 = Pm1[n];
x21 += cnm * t1; x21 += cnm * t1;
x22 += dnm * t1; x22 += dnm * t1;
}; }
}; }
if (m == 1) if (m == 1)
{ {
...@@ -323,7 +325,7 @@ void ODEPhaseSystem::EarthNonCentralAcc(Vect3 &acc, const Vect3 &pos) ...@@ -323,7 +325,7 @@ void ODEPhaseSystem::EarthNonCentralAcc(Vect3 &acc, const Vect3 &pos)
double *Pmt = Pm; double *Pmt = Pm;
Pm = Pm1; Pm = Pm1;
Pm1 = Pmt; Pm1 = Pmt;
}; }
double x01 = z0 + z1 + zr * z2; double x01 = z0 + z1 + zr * z2;
double x11 = EarthGravity->mu / (r * r); double x11 = EarthGravity->mu / (r * r);
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment