Die Lösung zu
\(\ddot x(t) + w^2 x(t) = 0\)
ist bekannt:
\(x_{h}(t)= A sin(\omega t) + Bcos(\omega t)\)
Und jetzt machst du Var. der Konstanten, der Einfachheit halber erstmal nur für A,daher
\(x_p (t)= A(t) sin(\omega t)\). Das setzt du in die DGL ein, es kommt eine neue DGL in
\( \ddot A(t), \dot A(t)\) heraus, welche man mit einer weiteren Variation der Konstanten lösen kann. Die Rechnung ist nicht so schwer, aber ist mir zuviel Schreibarbeit im Moment ;).