I am attempting to solve a 4th order pde in two variables on a square
$$u_{xxxx}+2u_{xxyy}+u_{yyyy}=0$$
with some Neumann boundary condition on the bottom $u_y=\sin x$ and $u_x=0$ on the top we may assume $u_{yy}=u_{xx}$.
What I have tried was to use separation of variables, $u=XY$
The problem is that when I do that I get the equation
$$X^{(4)}Y+2X^{(2)}Y^{(2)}+XY^{(4)}=0$$
I decided to divide by $XY$ but that leaves the middle term to be problematic. Can someone give me some advice for possible approaches here?