Hi,
wie von pleindespoir angesprochen stimmt Dein homogener Teil schon nicht.
Es ist x^2+1 = 0 und damit
y_(h) = c*sin(x) + d*cos(x)
Deshalb funktioniert auch Dein partikulärer Ansatz nicht. Hast ja nen Resonanzfall. Nutze
y_(p) = (ax+b)*e^x + p*x*sin(x) +q*x*cos(x)
Wenn ich mich auf die Schnelle nicht vertan habe, kommst Du auf:
y = c*sin(x) + d*cos(x) + e^x*x - e^x + 1/2*x*sin(x)
Grüße