Physic Labs

Frontier physics

Nonlinear dynamical systems, fixed points, and bifurcations

Sweep the parameter μ of the normal form x˙=μ−x2\dot{x} = \mu - x^2 to find the bifurcation where an equilibrium pair is born or annihilated, then probe stability with different initial conditions.

Advanced

Equipment

  • “Parameter μ” slider (−100…+100)
  • “Initial condition x₀” slider
  • Vector-field plot of $\dot{x}$ versus $x$ with evolution arrows
  • Plot of the solution x(t) integrated by Euler's method

Procedure

  1. Find the equilibria on the vector field

    Set μ > 0 and watch the curve x˙=μ−x2\dot{x} = \mu - x^2 cross the horizontal axis at two points: x∗=±μx^* = \pm\sqrt{\mu}. Increase μ and watch the pair move apart; lower μ toward 0 and see them approach each other.

  2. Probe stability with initial conditions

    Keep μ > 0 and drag “Initial condition x₀” across the upper equilibrium, then watch x(t). Trajectories starting from x0>−μx_0 > -\sqrt{\mu} converge to +μ+\sqrt{\mu} (stable); those starting below drift to −∞-\infty — confirming the upper fixed point is unstable, consistent with the sign of f′(x)=−2xf'(x) = -2x.

  3. Cross the saddle-node bifurcation

    Lower μ through 0: the two equilibria merge at x∗=0x^* = 0 and annihilate — a saddle-node bifurcation. For μ < 0 every trajectory escapes; as a check, record the μ value where the pair appears and compare with the algebraic prediction μc=0\mu_c = 0.

Simulation

Experiment history

The qualitative approach to dynamical systems began with Henri Poincaré: in the three-body problem (1889) and the « Méthodes nouvelles de la mécanique céleste » (1892–99) he studied solution curves, equilibria and stability directly instead of seeking analytic solutions — the seed of modern bifurcation theory. In the 1920s Balthasar van der Pol described self-oscillations in electronic circuits; in 1963 Edward Lorenz found chaotic motion in a three-equation atmosphere model, and in 1975–78 Mitchell Feigenbaum proved the universality of period-doubling cascades. The normal form x˙=μ−x2\dot{x} = \mu - x^2 in this simulation is the simplest template for a local bifurcation.

Related physicists

Related library topics