This commit is contained in:
Ceres 2025-10-20 15:58:03 +01:00
commit 791d876347
Signed by: ceres-sees-all
GPG key ID: 9814758436430045

View file

@ -72,12 +72,12 @@ def part2(velocityFactor, orbits):
dymdt = vy
dvxdt = -(Me*G*xm)/((xm**2+ym**2)**(3/2))
dvydt = -(Me*G*ym)/((xm**2+ym**2)**(3/2))
dxpdt = vpx
dypdt = vpy
dvpxdt = -((Me*G*xp)/((xp**2+yp**2)**(3/2)))-((Mm*G*xpm)/((xpm**2+ypm**2)**(3/2)))
dvpydt = -((Me*G*yp)/((xp**2+yp**2)**(3/2)))-((Mm*G*ypm)/((xpm**2+ypm**2)**(3/2)))
return (dxmdt, dymdt, dvxdt, dvydt, dxpdt, dypdt, dvpxdt, dvpydt)
@ -96,7 +96,7 @@ def part2(velocityFactor, orbits):
rpm = 10000000
vpm = sqrt((G*Mm)/rpm)
xm0 = rm
ym0 = 0
vx0 = 0