$$2\frac{\partial^2 z}{\partial x^2}-5\frac{\partial^2 z}{\partial x\partial y}+\frac{\partial^2 z}{\partial y^2} = 5\sin(2x+y)$$
is an inhomogeneous PDE.
$$\left( \frac{\partial}{\partial x}-2\frac{\partial}{\partial y}\right)\left( 2\frac{\partial}{\partial x}-\frac{\partial}{\partial y}\right)z(x,y)=5\sin(2x+y)$$
Solving the associated homogeneous PDE : $\left( \frac{\partial}{\partial x}-2\frac{\partial}{\partial y}\right)\left( 2\frac{\partial}{\partial x}-\frac{\partial}{\partial y}\right)Z(x,y)=0$
$$Z(x,y)=f(2x+y)+g(x+2y)$$
with any differentiable functions $f$ and $g$.
Search for a particular solution of the inhomogeneous PDE:
We could use the method of variation of parameter in searching a solution of the form $z=h(x,y)Z(x,y)$. This would be tiresome. By luck the form of the term $5\sin(2x+y)$ is the same as $f(2x+y)$.
So, it is much simpler to search a solution on the form: $c_1\:x\sin(2x+y)+c_2\:x\cos(2x+y)$. Putting it into the PDE allows to identify $c_1=0$ and $c_2=-\frac{5}{3}$. Thus, the general solution of the PDE is :
$$z(x,y)=f(2x+y)+g(x+2y)-\frac{5}{3}\cos(2x+y)$$
with any differentiable functions $f$ and $g$.