Live JS port of the Fortran superellipsoid_contact program — 4 anisotropic particles, gravity + wall/pair spring contacts, velocity-Verlet integration, forces via numerical gradient. Shapes drawn with the same radius function the physics uses.
a=0.4 · b=0.25 · ε=0.5
a=0.3 · b=0.3 · ε=1
a=0.5 · b=0.3 · ε=2
a=0.35 · b=0.2 · ε=3
Energy conservation (velocity Verlet)
The gold Total line is the energy-conservation check — a stable Total means the integrator is behaving. Drift indicates contact activity (springs storing/releasing energy) and numerical error.