diff --git a/Exercise 1/exercise1.py b/Exercise 1/exercise1.py
index 66228e2..ad23992 100644
--- a/Exercise 1/exercise1.py	
+++ b/Exercise 1/exercise1.py	
@@ -9,10 +9,10 @@ def part1(radFactor, orbits):
     def f_part1(t, state, Me, Mm, G):
         xm, ym, vx, vy = state
 
-        dvxdt = -(Me*G*xm)/((xm**2+ym**2)**(3/2))
         dxmdt = vx
-        dvydt = -(Me*G*ym)/((xm**2+ym**2)**(3/2))
         dymdt = vy
+        dvxdt = -(Me*G*xm)/((xm**2+ym**2)**(3/2))
+        dvydt = -(Me*G*ym)/((xm**2+ym**2)**(3/2))
 
         return (dxmdt, dymdt, dvxdt, dvydt)
 
@@ -70,12 +70,12 @@ def part2(radFactor, 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)
 
 
@@ -94,7 +94,7 @@ def part2(radFactor, orbits):
 
     rpm = 10000000
     vpm = sqrt((G*Mm)/rpm)
-    
+
     xm0 = rm
     ym0 = 0
     vx0 = 0