berechne erst die homogene LÖsung, indem Du die rechte Seite 0 setzt.
Du kommst auf die Eigenwerte λ_(1) = -5 und λ_(2) = 1.
Da die rechte Seite keine reellen e-Funktion beinhaltet, haben wir schonmal keinen Resonanzfall vorliegen.
Der zu wählende Ansatz ist y_(p) = a*sin(x) + b*cos(x).
Dieser ist zweimal abzuleiten und in die DGL einzusetzen. Ein Koeffizientenvergleich liefert a = -2 und b = -3.
Gesamtlösung ist folglich:
y = c_(1)*e^{-5x} + c_(2)*e^x - 2*sin(x) - 3*cos(x)
Grüße