Numerische Lösung DGL des harmonischen Oszillator