Using the substitution $u(x, y)=v(x, y) e^{-x^{2}}$, find $u(x, y)$ that satisfies

$$u_{x}+x u_{y}+2 x u=e^{-x^{2}}$$

with boundary data $u(0, y)=y e^{-y^{2}}$.