hi prof,

not sure why but total energy is not conserved in this model, can help debug?

ke = 0.5*mtest*vxtest*vxtest;
//pe = -G*mtest*m2/Math.abs(x2-xtest)-G*mtest*m1/Math.abs(x1-xtest);
pe = Vy*mtest;
te = ke+pe;

i have been looking at this on and off, but cannot figure out what is wrong. te should be a constant value

the total energy of the mtest, te decays which is not suppose to happen. weird

:)

appreciate any help
thanks!