Phase Portraits
Tool 1 of 28
Direction Field
\(\dfrac{dy}{dx} = f(x,y)\)
Plot the slope field of a first-order ODE and overlay solution curves traced from any initial point. Type \(f(x,y)\) in standard math syntax; the canvas redraws a grid of tick marks aligned with the local slope and integrates trajectories on click.
Explore →
Tool 2 of 28
Phase Line Diagram
\(\dfrac{dy}{dt} = f(y)\)
Analyze equilibria and stability for autonomous first-order ODEs. Enter \(f(y)\), read off the fixed points where \(f(y)=0\), and classify each as stable, unstable, or semi-stable from the sign of \(f\) on either side. Solution curves and arrows on the line update in real time.
Explore →
Tool 3 of 28
Linear Phase Portraits
\(\dot{\mathbf{x}} = A\mathbf{x},\quad A\in\mathbb{R}^{2\times 2}\)
Classify equilibria of 2×2 linear systems — spirals, centers, nodes, and saddle points — by specifying the matrix \(A\) and reading off eigenvalues. Click anywhere on the canvas to integrate a trajectory from that initial condition.
Explore →
Tool 4 of 28
Nonlinear Phase Portraits
\(\dot{x} = f(x,y),\quad \dot{y} = g(x,y)\)
Plot direction fields and trajectories for arbitrary planar nonlinear systems. Visualize fixed points, limit cycles, and separatrices with RK4 and SDIRK2 hybrid solvers for stiff regimes. Click the canvas to integrate from any initial condition.
Explore →
Pendulum Systems
Tool 5 of 28
Damped Pendulum Simulator
\(\ddot{\theta} + \dfrac{b}{m}\dot{\theta} + \dfrac{g}{L}\sin\theta = 0\)
Explore the motion and phase portrait of a damped nonlinear pendulum. Compare small-angle (linear) and large-angle (nonlinear) dynamics.
Explore →
Tool 6 of 28
Elastic (Spring) Pendulum
\(\mathcal{L} = \tfrac{1}{2}m(\dot{r}^2 + r^2\dot{\theta}^2) + mgr\cos\theta - \tfrac{1}{2}k(r-L_0)^2\)
Simulate a mass on a spring that can both stretch and swing. Explore energy transfer between radial and angular modes, autoparametric resonance, and quasi-periodic orbits.
Explore →
Tool 7 of 28
Coupled Pendula with Spring
\(\ddot{\theta}_1 = -\dfrac{g}{L}\sin\theta_1 - \kappa(\theta_1 - \theta_2)\)
Simulate two pendulums coupled by a spring. Explore normal modes, energy transfer, and beat phenomena.
Explore →
Tool 8 of 28
Multi-Pendulum Simulator
\(\dfrac{d}{dt}\dfrac{\partial \mathcal{L}}{\partial \dot{\theta}_i} - \dfrac{\partial \mathcal{L}}{\partial \theta_i} = 0\)
Simulate double and triple pendulums with chaotic trajectory visualization. Observe sensitivity to initial conditions in Lagrangian mechanics.
Explore →
Tool 9 of 28
Pendulum Phase Array
\(L_k = \dfrac{g\,T^2}{4\pi^2(n_0+k)^2}, \quad k = 0,\ldots,N{-}1\)
\(N\) independent pendulums tuned so pendulum \(k\) completes exactly \((n_0+k)\) oscillations in period \(T\). The pendulums are not coupled — apparent collective motion (snake, butterfly, in-phase, wave) emerges purely from differential phase. Features a 3D perspective view, horizontal displacement wave strip, real-time pattern detection, dark/light/high-contrast modes, and a preset library.
Explore →
Tool 10 of 28
Swinging Atwood Machine
\(\ddot r = \dfrac{m\,r\,\dot\theta^2 - g(M-m\cos\theta)}{M+m}\)
A heavy counterweight \(M\) shares an inextensible string with a smaller bob \(m\) over an idealized pulley. The bob swings as a pendulum of variable length \(r\). Closed cusped orbits at \(M/m=3\); chaos otherwise.
Explore →
Mass-Spring Systems
Tool 11 of 28
Mass-Spring-Damping Simulator
\(m\ddot{x} + c\dot{x} + kx = 0\)
Simulate underdamped, critically damped, and overdamped motion. Observe displacement and phase-plane trajectories in real time.
Explore →
Tool 12 of 28
Multiple Degree of Freedom Mass Spring
\(M\ddot{\mathbf{x}} + C\dot{\mathbf{x}} + K\mathbf{x} = \mathbf{0}\)
Simulate coupled mass-spring systems and explore normal modes of vibration. Visualize how energy transfers between coupled oscillators.
Explore →
Tool 13 of 28
Forced Vibrations and Resonance
\(m\ddot{x} + c\dot{x} + kx = F_0\cos(\omega t)\)
Observe transient and steady-state behavior, including resonance. Explore how driving frequency affects amplitude response.
Explore →
Fourier & Ecology
Tool 14 of 28
Fourier Series Approximation
\(f(x) = \dfrac{a_0}{2} + \displaystyle\sum_{n=1}^{N}\bigl(a_n\cos nx + b_n\sin nx\bigr)\)
See how partial sums of sine and cosine terms approximate periodic functions. Adjust the number of terms and observe convergence.
Explore →
Tool 15 of 28
Lotka-Volterra Model
\(\dfrac{dx}{dt} = \alpha x - \beta xy, \quad \dfrac{dy}{dt} = \delta xy - \gamma y\)
Simulate classic predator-prey population cycles in phase space. Observe periodic orbits and the relationship between predator and prey densities.
Explore →
Tool 16 of 28
Lotka-Volterra with Carrying Capacity
\(\dfrac{dx}{dt} = \alpha x\!\left(1 - \dfrac{x}{K}\right) - \beta xy\)
Extend the predator-prey model with logistic growth limits. Explore how carrying capacity stabilizes or changes system dynamics.
Explore →
Chaotic Systems & Nonlinear Oscillators
Tool 17 of 28
Van der Pol Oscillator
\(\dfrac{d^2x}{dt^2} + \mu(x^2 - 1)\dfrac{dx}{dt} + x = 0\)
Explore limit cycles and relaxation oscillations of the Van der Pol oscillator. Visualize phase portraits and time series with RK4/SDIRK2 hybrid solver for stiff regimes.
Explore →
Tool 18 of 28
Duffing Oscillator
\(\dfrac{d^2x}{dt^2} + \delta\dfrac{dx}{dt} + \alpha x + \beta x^3 = \gamma\cos(\omega t)\)
Explore the Duffing oscillator — a nonlinear second-order ODE with double-well potential, hardening/softening spring regimes, and chaotic forcing. Visualize phase portraits, time series, and potential energy with an adaptive RK45 solver.
Explore →
Tool 19 of 28
Lorenz System
\(\dot{x} = \sigma(y-x),\quad \dot{y} = x(\rho-z)-y,\quad \dot{z} = xy - \beta z\)
Visualize the famous chaotic attractor and sensitive dependence on initial conditions. Explore how parameters \(\sigma\), \(\rho\), \(\beta\) shape the dynamics.
Explore →
Tool 20 of 28
Rössler System
\(\dot{x} = -y - z,\quad \dot{y} = x + ay,\quad \dot{z} = b + z(x - c)\)
Explore the Rössler attractor — a 3D chaotic system with a single-scroll topology and a clean period-doubling route to chaos. Sweep the parameter \(c\) to watch periodic orbits split into period-2, period-4 cycles and finally the strange attractor.
Explore →
Tool 21 of 28
Malkus Waterwheel
\(\dot{a}_1 = \omega b_1 - K a_1 + q_1,\ \ \dot{b}_1 = -\omega a_1 - K b_1,\ \ I\dot{\omega} = \pi g R\, b_1 - \nu\omega\)
Mechanical analogue of the Lorenz system: a tilted wheel of leaky buckets whose angular velocity reverses chaotically. Animated rotating wheel with leaky buckets, time-series, and a 3D phase-space view of the butterfly attractor.
Explore →
Tool 22 of 28
Poincaré Section Explorer
\(\mathcal{P}: \mathbf{x}\in\Sigma\mapsto \varphi_{\tau(\mathbf{x})}(\mathbf{x})\in\Sigma\)
High-resolution stroboscopic and hyperplane Poincaré sections for Duffing, the driven pendulum, Lorenz, Hénon–Heiles, the double pendulum, and the Chirikov standard map. Adaptive Tsit5 integrator in a Web Worker; multi-trajectory runs with per-orbit colouring.
Explore →
Potentials, Bifurcations & Discrete Maps
Tool 23 of 28
Particle in a Cubic Potential
\(V(x) = ax^3 + bx, \quad \ddot{x} = -V'(x)\)
Study motion and phase space trajectories in an asymmetric potential well. Explore bounded and unbounded orbits depending on initial energy.
Explore →
Tool 24 of 28
Particle in a Double Well Potential
\(V(x) = ax^4 - bx^2, \quad \ddot{x} = -V'(x)\)
Explore bistability in a symmetric double well. Observe how initial energy determines which well the particle occupies.
Explore →
Tool 25 of 28
Bifurcation Diagram
\(\text{plot }\{(r,\,x_n)\}\text{ as }r\text{ varies}\)
Watch period-doubling cascades, the Feigenbaum constant \(\delta\approx4.669\), and chaos emerge as the parameter sweeps. Dots coloured by Lyapunov exponent \(\lambda\). Lyapunov exponent plotted separately below. Custom map support.
Explore →
Tool 26 of 28
Bifurcation Family Explorer
\(\dot{x} = r + x^2,\ \ rx - x^2,\ \ rx - x^3,\ \ rx + x^3 - x^5\)
Coupled phase line and \((r,x^*)\) bifurcation diagram for the four canonical 1D bifurcations: saddle-node, transcritical, supercritical and subcritical pitchfork. Observe fixed-point creation, exchange of stability, symmetry breaking, and hysteresis.
Explore →
Tool 27 of 28
Logistic Map Explorer
\(x_{n+1} = r\,x_n(1-x_n)\)
Cobweb diagrams for the logistic map, tent map, and sine map. Explore the full period-doubling cascade — from stable fixed points through period-2, 4, 8 cycles to chaos. Eight quick-preset parameter values, real-time Lyapunov exponent, regime indicator, and step-by-step orbit animation.
Explore →
Tool 28 of 28
Cobweb Diagram
\(x_{n+1} = f(x_n)\)
Visualize cobweb diagrams for fixed-point iteration and discrete dynamical systems. Explore convergence, periodic orbits, and chaos with interactive map selection and parameter control.
Explore →