While playing around with a number of simple integrators I wrote a quick 1 dimensional spring demo. Four point spring unit masses are animated using various methods of evaluation (simple Euler, Newton-Stormer-Verlet, 4th order Runge-Kutta and a simple explicit Verlet integrator.
All four methods have the tendency to explode over time, although the Euler integrator does that the easiest. Somewhat to my surprise both Verlet integrators are more stable than the 4th order Runge Kutta integrator.
Source code can be viewed online. Compile with
ocamlopt unix.cmxa graphics.cmxa physics.ml -o physics or
ocamlc unix.cma graphics.cma physics.ml -o physics
No comments yet.
Leave a comment
Sorry, the comment form is closed at this time.