Betrachte stattdessen die Differentialgleichungen für
$$u = y_1 + y_2 \\ v = y_1 - y_2.$$
Sie lauten:
$$u' = y_1' + y_2' = y_1 + y_2 + 1 + \sin(x)\\ u' = u + 1 + sin(x)$$
und
$$v' = -v + 1 - sin(x)$$
So sind die Differentialgleichungen also entkoppelt und können wie zwei separate lineare, inhomogene Differentialgleichungen gelöst werden.
Anschließend erhält man die gesuchten Funktionen mittels
$$y_1 = \frac{u+v}{2} \\ y_2 = \frac{u-v}{2}$$