Click to focus on the app, then click again to add a particle, and click again to set its velocity.

p=M2k e=1+2EM2k2 r(φ)=p1+ecos(φφ0) φ˙=Mr2 φ0=φArcTanFull(r˙pMe,(pr1)1e)

for a unit mass particle at angle ϕ with respect to some reference axis, angular momentum M, distance r from the center, orbital energy E, where p and e are the parameter and eccentricity of the ellipse, and the particle is in a potential k/r