Harmonischer Oszillator
Bewgungsgleichung: a[x] = - k x

Verbesserter Algorithmus:

x1 = x0 + v0 dt + a[x0]/2 dt^2 ... dt^2-Term
v1 = v0 + a[ (x0+x1)/2 ] dt    ... Mittelwertbildung!

passt sagenhaft gut!

Beachte: a[ (x0+x1)/2 ] = (a[x0]+a[x1])/2 in diesem Fall!