Die Lösung zu
x¨(t)+w2x(t)=0
ist bekannt:
xh(t)=Asin(ωt)+Bcos(ωt)
Und jetzt machst du Var. der Konstanten, der Einfachheit halber erstmal nur für A,daher
xp(t)=A(t)sin(ωt). Das setzt du in die DGL ein, es kommt eine neue DGL in
A¨(t),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 ;).