How to solve $u_{xy}+u_x+u_y=2$?

$$u_{xy}+u_x+u_y=2$$ Let $\quad u(x,y)=x+y+v(x,y)\quad$ which leads to : $$v_{xy}+v_x+v_y=0$$ Search for particular solutions of the form $v=X(x)Y(y)$

$X'Y'+X'Y+XY'=0\quad;\quad\frac{X'}{X}=-\frac{Y'}{Y+Y'}=\lambda \qquad\begin{cases} X=e^{\lambda x} \\ Y=e^{-\frac{\lambda}{1+\lambda}y}\end{cases}$

Particular solutions : $\quad v_\lambda(x,y)=e^{\lambda x-\frac{\lambda}{1+\lambda}y}$

The general solution is any linear combination of the particular solutions.

General solution expressed on the form of integral : $$v(x,y)=\int f(\lambda)e^{\lambda x-\frac{\lambda}{1+\lambda}y}d\lambda$$ $f(\lambda)$ is an arbitrary function, in so far the integral be convergent. $$u(x,y)=x+y+\int f(\lambda)e^{\lambda x-\frac{\lambda}{1+\lambda}y}d\lambda$$ Or, equivalently : $$u(x,y)=x+y+\int g(\mu)e^{-\frac{\mu}{1+\mu}x+\mu y}d\mu$$ The function $f(\lambda)$ , or $g(\mu)$ , as well as the bounds of the integral, have to be determined according to the boundary conditions. This generally draw to solve an integral equation. Without well defined boundary conditions no further calculus is possible, even not to say if it is possible to analytically solve it and if there is a closed form for the solution.