To solve a non-homogeneous linear PDE $\displaystyle \frac{\partial^2 z}{\partial x^2}+\frac{\partial^2 z}{\partial x \, \partial y}+\frac{\partial z}{\partial y}-z=e^{-x}$
My Attempt: Putting $\displaystyle D=\frac{\partial z}{\partial x}$ and$\displaystyle D’=\frac{\partial z}{\partial y}$, we have $\displaystyle (D^2+DD’+D’-1)z=e^{-x}$
$$(D+1)(D+D’-1)=e^{-x}$$
C.F. $\displaystyle \implies z=e^{-x}\phi_1(y) +e^x\phi_2(y-x))$
$$ \mathrm{P.I.}=\frac{1}{f(D,D')}e^{-x}=\frac{1}{(1+1)(D+D'-1)}$$
$$ P.I.=\frac{1}{2(D+D'-1)}e^{-x}$$ (Putting $D=1$ into first factor)
Putting $D=1$ in next factor will get us zero, hence, we move to $$ \frac{1}{(D-mD')}F(x,y)=\int F(x,c-mx)dx=\int e^{-x}dx=-e^{-x}$$
So we finally get, C.F. + P.I. $=\displaystyle e^{-x}\phi_1(y) +e^x\phi_2(y-x))- \frac{e^{-x}}{2}$
The given answer: $$ e^{-x}\phi_1(y) +e^x\phi_2(y-x))- \frac{xe^{-x}}{2}$$
Where does the x come in here as per the given answer? Where am I going wrong?