離散事象を含む微分方程式 

正方形の箱の側面にぶつかると方向を変えるボールの動きをモデル化する．

In[1]:=

X sol = NDSolve[{x'[t] == a[t], y'[t] == b[t], x[0] == 0, y[0] == 0, a[0] == 1, b[0] == Sqrt[2], WhenEvent[x[t]^2 - 1, a[t] -> -a[t]], WhenEvent[y[t]^2 - 1, b[t] -> -b[t]]}, {x, y}, {t, 0, 100}, DiscreteVariables -> {a, b}]; Map[ParametricPlot[{x[t], y[t]} /. sol, {t, 0, #}, Frame -> True, FrameTicks -> None, PlotRange -> 1, Axes -> False] &, {5, 20, 40, 80}]

Out[1]=