ceres-sees-all/Uni-Computing
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Add Tex and update scripts

Browse source
This commit is contained in:
Ceres 2025-10-21 08:53:34 +01:00
parent 791d876347
commit b5b7b31a87
Signed by: ceres-sees-all
GPG key ID: 9814758436430045
7 changed files with 147 additions and 9 deletions
Download patch file Download diff file Expand all files Collapse all files

18
Exercise 1/exercise1.py
View file

@ -61,19 +61,19 @@ def part1(velocityFactor, orbits):
plt.show()
def part2(velocityFactor, orbits):
#Part one Diffeq
def f_part1(t, state, Me, Mm, G):
xm, ym, vx, vy, xp, yp, vpx, vpy = state
#Part two Diffeq
def f_part2(t, state, Me, Mm, G):
xm, ym, vx, vy, xp, yp, vpx, vpy = state #all input values
xpm = xp-xm
ypm = yp-ym
dxmdt = vx
dxmdt = vx #moon diffeqs
dymdt = vy
dvxdt = -(Me*G*xm)/((xm**2+ym**2)**(3/2))
dvxdt = -(Me*G*xm)/((xm**2+ym**2)**(3/2))
dvydt = -(Me*G*ym)/((xm**2+ym**2)**(3/2))
dxpdt = vpx
dxpdt = vpx #probe diffeqs
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)))
@ -102,8 +102,8 @@ def part2(velocityFactor, orbits):
vx0 = 0
vy0 = vm
xp0 = rpm+rm
yp0 = 0
xp0 = rpm+rm #as position is relative to the earth
yp0 = 0 #not to the moon
vpx0 = 0
vpy0 = vm+vpm
@ -111,7 +111,7 @@ def part2(velocityFactor, orbits):
atol = 1e-9
#Solver
results = solve_ivp(f_part1, (t_min,t_max), (xm0, ym0, vx0, vy0, xp0, yp0, vpx0, vpy0), args=(Me, Mm, G), t_eval=t, atol=atol, rtol=rtol)
results = solve_ivp(f_part2, (t_min,t_max), (xm0, ym0, vx0, vy0, xp0, yp0, vpx0, vpy0), args=(Me, Mm, G), t_eval=t, atol=atol, rtol=rtol)
#Graph plotting