Aus der PDG $$ u_{xy} + u_x = 1 - x - y $$ folgt mit \( v = u_x \) $$ v_y + v = 1 - x - y $$ und diese Gleichung hat die Lösung
$$ v = C(x) e^{-y} -x -y +2 $$ mit einer noch unbekannten Funktion \( C (x) \)
Also gilt $$ u(x,y) = e^{-y} \int C(x) dx - \frac{x^2}{2} - x y +2 x + G(y) $$ mit noch unbekannter Funktion \( G(y) \)
Aus den Randbedingungen folgt $$ 0 = u(x,0) = \int C(x) dx - \frac{x^2}{2} + 2x + G(0) $$ Also $$ \int C(x) dx = \frac{x^2}{2} - 2x - G(0) $$ und damit $$ u(x,y) = e^{-y} \left( \frac{x^2}{2} - 2x - G(0) \right) - \frac{x^2}{2} - x y +2 x + G(y) $$ Aus $$ 0 = u(0,y) = -G(0) e^{-y} + G(y) $$ folgt $$ G(y) = G(0) e^{-y} $$ Alles eingesetzt ergibt
$$ u(x,y) = \left( e^{-y} - 1 \right) \left( \frac{x^2}{2} - 2x \right) - xy $$