magix2

To reproduce the problem: ````julia using OrdinaryDiffEq ω = 1 x₀ = [0.0] dx₀ = [π/2] tspan = (0.0, 2π) function harmonic_oscillator(ddu,du,u,ω,t) ddu .= -ω^2 * u end prob =...