Automatic discontinuity handling

Voorbeeld 2 stappen, op volgorde uitgevoerd in één sessie. Elke stap is tegen de engine uitgevoerd en de uitvoer hieronder is wat hij opleverde.

Stap 1

sol = NDSolve[{y'[t] == 15 / 4 - Floor[y[t]] + DiracDelta[t - 2], y[0] == .1}, y, {t, 0, 6}];

Geen uitvoer - deze stap zet iets klaar voor de volgende.

Stap 2

Plot[Evaluate[{y[t], y'[t]} /. sol], {t, 0, 6}, PlotLegends -> {"y(t)", "y'(t)"}]

Uitvoer

-Graphics-

Gebruikte functies