mirror of
https://codeberg.org/Yael-II/MSc2-Project-Chaos.git
synced 2026-03-15 04:16:26 +01:00
exam version
This commit is contained in:
BIN
Bhat_Moussouni_Numerical_Simulation_Report.pdf
Normal file
BIN
Bhat_Moussouni_Numerical_Simulation_Report.pdf
Normal file
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
BIN
Figs/pot_hh.png
BIN
Figs/pot_hh.png
Binary file not shown.
|
Before Width: | Height: | Size: 116 KiB After Width: | Height: | Size: 118 KiB |
Binary file not shown.
@@ -140,7 +140,7 @@ for h in h_range:
|
|||||||
color=colors['RK4'],
|
color=colors['RK4'],
|
||||||
label="RK4")
|
label="RK4")
|
||||||
|
|
||||||
ax.set_title("$\\Var{{t}} = {:.4f}$".format(h))
|
ax.set_title("$h = {:.4f}$".format(h))
|
||||||
ax.set_xlabel("$x$")
|
ax.set_xlabel("$x$")
|
||||||
ax.set_ylabel("$y$")
|
ax.set_ylabel("$y$")
|
||||||
ax.set_aspect("equal")
|
ax.set_aspect("equal")
|
||||||
@@ -177,7 +177,7 @@ for h in h_range:
|
|||||||
|
|
||||||
axs[i].set_aspect("equal")
|
axs[i].set_aspect("equal")
|
||||||
|
|
||||||
fig.suptitle("$\\Var{{t}} = {:.4f}$".format(h))
|
#fig.suptitle("$\\Var{{t}} = {:.4f}$".format(h))
|
||||||
axs[0].set_xlabel("$x$")
|
axs[0].set_xlabel("$x$")
|
||||||
axs[0].set_ylabel("$y$")
|
axs[0].set_ylabel("$y$")
|
||||||
axs[0].legend(loc="upper left")
|
axs[0].legend(loc="upper left")
|
||||||
|
|||||||
@@ -64,7 +64,7 @@ def hh(W):
|
|||||||
Y = W[1]
|
Y = W[1]
|
||||||
POT = pot.hh_potential(W, position_only=True)
|
POT = pot.hh_potential(W, position_only=True)
|
||||||
fig, ax = plt.subplots(1)
|
fig, ax = plt.subplots(1)
|
||||||
ax.set_title("Hénon–Heils potential")
|
ax.set_title("Hénon--Heils potential")
|
||||||
pcm = ax.pcolormesh(X, Y, POT)
|
pcm = ax.pcolormesh(X, Y, POT)
|
||||||
fig.colorbar(pcm, label="potential")
|
fig.colorbar(pcm, label="potential")
|
||||||
ax.set_aspect("equal")
|
ax.set_aspect("equal")
|
||||||
@@ -73,7 +73,6 @@ def hh(W):
|
|||||||
|
|
||||||
fig.savefig("Figs/pot_hh.png")
|
fig.savefig("Figs/pot_hh.png")
|
||||||
return 0
|
return 0
|
||||||
|
|
||||||
if __name__ == "__main__":
|
if __name__ == "__main__":
|
||||||
kepler(W)
|
kepler(W)
|
||||||
hh(W)
|
hh(W)
|
||||||
|
|||||||
Reference in New Issue
Block a user