Vi skriver om v'' + sin v = 0 som ett första ordningens system genom att införa funktionerna
y(x) = v(x) samt z(x) = v'(x) = y'(x). Systemet blir alltså
y' = z,         y(0) = v0
z' = -sin y,  z(0) = 0

Vi vill också jämföra med lösningen till v'' + v = 0. Jag löser därf två system och har definierat funktionerna
f, g och h som högerled i mina system. Funktionen f används i båda systemen.
Jag plottar enbart vinkeln (inte dess derivata) och vi ser, genom att variera glidaren, att överensstämmelsen
är god för små v0. Ett sätt att förklara detta är att studera Taylorutvecklingen, sin v ≈ v - v^3 / 6.
Avsaknaden av v^2-term gör att sinv ≈ v stämmer mycket bra för små v. Felet ökar med tiden.

f(x, y, z) = z
g(x, y, z) = (-sin(y))
h(x, y, z) = (-y)
v0 = 0.6 (slider from 0 to 1, step 0.1)
numeriskIntegral1, numeriskIntegral2 = NSolveODE({f, g}, 0, {v0, 0}, 60)
numeriskIntegral3, numeriskIntegral4 = NSolveODE({f, h}, 0, {v0, 0}, 60)