This problem came up when I was studying my final exam for my PDE class.
Find the solution to the first-order PDE $$u_x+xu_y=e^{x+y}$$ with the initial conditions $u(0, y)=y$.
I have only learnt two ways about find solutions to a linear PDE. The first method is by transforming the coordinates such that the equation becomes more simple, and the second method is using the method of characteristics(or this is what my professor calls it). It seems that there is some more general version of the "method of characteristics" that I am using like in this question. It seems that they introduce some dummy variable, and change it to a system of ODE's, but up to my knowledge, the method of characteristics is the method of finding the characteristic curves, in this case, which is the family of curves satisfying $\frac{dy}{dx}=\frac{x}{1}$ so that the function must be constant on these curves. We have $\frac{dy}{dx}=x$, $y-\frac{x^2}{2}=c$ in this case, so the homogeneous solution is given as $f(y-\frac{x^2}{2})$. All I need to find now is a particular solution to the non homogeneous equation, and I would add to the homogeneous solution. I have tried functions of the form $f(x)e^{x+y}, f(y)e^{x+y}, (ax+by)e^{x+y}$, but they did not work. I have also tried transforming the coordinates to $$\xi= y-\frac{1}{2}x^2, \quad \eta=x+y, x, y, f(x), f(y) (\text{and several variants})$$ but they did not work.
So, my question is, is there an elagant way to find either a particular solution to the non homogeneous solution, or a elegant transformation that makes the equation much simple?
Please tell me if some parts needs clarification, or edit. Thank you in advance.