Commit 50524635 authored by Alexander Lapshin's avatar Alexander Lapshin

.

parent 9c684b5b
......@@ -226,8 +226,8 @@ void ODEPhaseSystem::EarthNonCentralAcc(Vect3 &acc, const Vect3 &pos)
{
ro1 *= ro;
Ron[n] = ro1;
//cn0 = EarthGravity->c(n, 0) * ro1;
cn0 = CC[0][n] * ro1;
cn0 = EarthGravity->c(n, 0) * ro1;
//cn0 = CC[0][n] * ro1;
x0 += cn0 * (n + 1) * p0;
x2 += cn0 * Pm[n];
if (n == Ng)
......@@ -271,11 +271,11 @@ void ODEPhaseSystem::EarthNonCentralAcc(Vect3 &acc, const Vect3 &pos)
{
t1 = Ron[n];
//double cnm = EarthGravity->c(n, m) * t1;
double cnm = CC[m][n] * t1;
double cnm = EarthGravity->c(n, m) * t1;
//double cnm = CC[m][n] * t1;
//double dnm = EarthGravity->d(n, m) * t1;
double dnm = DD[m][n] * t1;
double dnm = EarthGravity->d(n, m) * t1;
//double dnm = DD[m][n] * t1;
t1 = Pm[n];
......
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