The showcase

Four simulations, zero compromises.

Open to everyone, no account required. Every simulation in the full school library is held to this standard: accurate numerics, live graphs, real mathematics.

Classical Mechanics · GravitationAdvanced

CBSE 11UniversityExplore

N-Body Orbital Mechanics

A full gravitational N-body sandbox integrated with a symplectic velocity-Verlet scheme. Fly through Kepler orbits, the figure-8 three-body choreography, binary stars and gravitational slingshots while total energy and angular momentum are tracked live to numerical precision.

\ddot{\vec{r}}_i = \sum_{j \neq i} \frac{G m_j (\vec{r}_j - \vec{r}_i)}{|\vec{r}_j - \vec{r}_i|^3}

Symplectic velocity-Verlet integrator
Live energy & angular momentum conservation graphs
Per-body position, velocity and force readouts
Figure-8 three-body choreography preset

Waves & OscillationsIntermediate

CBSE 12University

Ripple Tank: Interference & Diffraction

The two-dimensional wave equation solved live on a 63,000-cell finite-difference grid — no scripted ripples. Double-slit fringes land exactly where the path-difference condition puts them, measured on an intensity screen against the Huygens–Fresnel prediction. Pour your own waves with the pointer, push a Doppler source past the wave speed, and watch a Mach cone form at arcsin(c/v).

\frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u

Leapfrog FDTD solver, dispersion verified to 0.02%
Measured fringes overlaid on the Huygens–Fresnel prediction
Double slit, single slit, two-source and Doppler presets
Click & drag to pour waves — supersonic wakes included

ElectromagnetismIntermediate

CBSE 12ICSE 12University

Electric Fields & Field Lines

Place positive and negative charges and watch the electric field draw itself — lines streaming out of the positives and into the negatives, exactly along the direction a small test charge would feel. Everything is a true superposition of Coulomb fields: the density of the lines shows the field strength, glowing markers flow along them from + to −, and the field vector and potential update live under your cursor. The classic dipole, like-charge and quadrupole patterns all emerge from the same inverse-square law.

\vec{E}(\vec{r}) = \sum_i \frac{k\,q_i\,(\vec{r}-\vec{r}_i)}{|\vec{r}-\vec{r}_i|^3}

Real superposition of Coulomb fields — E = kΣ q r̂ / r²
Field lines traced by integrating along E (line density ∝ field strength)
Drag charges live; dipole, like-charge & quadrupole presets
Field vector and potential read out live at the cursor

Classical Mechanics · ChaosAdvanced

UniversityExplore

Double Pendulum Chaos

The full nonlinear equations of motion derived from the Lagrangian, integrated with RK4. A twin system starting 10⁻⁴ rad away demonstrates sensitive dependence on initial conditions, with exponential trajectory divergence plotted on a log scale alongside the phase-space portrait.

\mathcal{L} = \tfrac{1}{2}(m_1+m_2)l_1^2\dot{\theta}_1^2 + \tfrac{1}{2}m_2 l_2^2\dot{\theta}_2^2 + m_2 l_1 l_2 \dot{\theta}_1\dot{\theta}_2\cos(\theta_1-\theta_2) + \cdots

Exact nonlinear equations from the Lagrangian
Twin-system chaos demonstration (Δθ₀ = 10⁻⁴)
Phase-space portrait and log-scale divergence plot
Energy conservation monitor

N-Body Orbital Mechanics

\ddot{\vec{r}}_i = \sum_{j \neq i} \frac{G m_j (\vec{r}_j - \vec{r}_i)}{|\vec{r}_j - \vec{r}_i|^3}

Symplectic velocity-Verlet integrator

Live energy & angular momentum conservation graphs

Per-body position, velocity and force readouts

Figure-8 three-body choreography preset

Ripple Tank: Interference & Diffraction

\frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u

Leapfrog FDTD solver, dispersion verified to 0.02%

Measured fringes overlaid on the Huygens–Fresnel prediction

Double slit, single slit, two-source and Doppler presets

Click & drag to pour waves — supersonic wakes included

Electric Fields & Field Lines

\vec{E}(\vec{r}) = \sum_i \frac{k\,q_i\,(\vec{r}-\vec{r}_i)}{|\vec{r}-\vec{r}_i|^3}

Real superposition of Coulomb fields — E = kΣ q r̂ / r²

Field lines traced by integrating along E (line density ∝ field strength)

Drag charges live; dipole, like-charge & quadrupole presets

Field vector and potential read out live at the cursor

Double Pendulum Chaos

\mathcal{L} = \tfrac{1}{2}(m_1+m_2)l_1^2\dot{\theta}_1^2 + \tfrac{1}{2}m_2 l_2^2\dot{\theta}_2^2 + m_2 l_1 l_2 \dot{\theta}_1\dot{\theta}_2\cos(\theta_1-\theta_2) + \cdots

Exact nonlinear equations from the Lagrangian

Twin-system chaos demonstration (Δθ₀ = 10⁻⁴)

Phase-space portrait and log-scale divergence plot

Energy conservation monitor