diff --git a/Exercise 1/exercise1.py b/Exercise 1/exercise1.py index 5c0e797..c297513 100644 --- a/Exercise 1/exercise1.py +++ b/Exercise 1/exercise1.py @@ -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