Solve some ordinary differential equations (ODE).
Winston is used for plotting, but using Winston
takes some time.
Github: https://github.com/nolta/Winston.jl
Documents: https://winston.readthedocs.io/
Comparison of some geometric integrators, Euler, symplectic-Euler, velocity-Verlet, position-Verlet and leapfrog methods.
Package: GeometricIntegrators.jl
examples of usage:
Symplectic-Euler, velocity-Verlet, position-Verlet and leapfrog are so-called symplectic integrators.
Test them with h=0.4
and n=10000
.
Runge-Kutta: rk1.jl
Runge-Kutta method is not a symplectic integrator.
Confirm that energy is not coserved with h=0.4
and n=10000
.
Runge-Kutta method is good for simulations of acceleration: acc.jl
$ LANG=C a2ps --header='Printed by t-nissie' --media=a4 --prologue=color\
--portrait --columns=1 --margin=3 --borders=off -f 11.0 --pretty-print=julia\
-o RungeKutta.ps ~/.julia/v0.4/RungeKutta/src/RungeKutta.jl
$ ~/scripts/PsDuplex RungeKutta.ps | lpr -l