Solution (source code)

= Solution

<Separation of variables> with $u=X(x)Y(y)$ and the three homogeneous boundary conditions gives $X_n=\sin(n\pi x)$ and $Y_n=\sinh(n\pi y)$. The boundary value at $y=1$ contains only the third <Fourier sine series> mode. Therefore
$$
\boxed{u(x,y)=2\frac{\sinh(3\pi y)}{\sinh(3\pi)}\sin(3\pi x).}
$$
This <harmonic function> has all four prescribed boundary values.