I want to solve a PDE coupled with an ODE:
$$ \nabla^2 T(x,y)=0 $$
defined over a rectangular region $x\in[0,L], y\in[0,l]$ subjected to:
$$\frac{\partial T(0,y)}{\partial x}=\frac{\partial T(L,y)}{\partial x}=0 \tag 1$$
$$\frac{\partial T(x,0)}{\partial y}=\gamma \tag 2$$
and on the top face $(y=l)$:
$$\frac{\partial T(x,l)}{\partial y}=\beta \Bigg[T(x,l)-t\Bigg] \tag 3$$
$t$ is not a constant and the equation governing it is:
$$\frac{\partial t}{\partial x}+\alpha(t-T(x,l))=0 \tag 4$$
where it is known that $t(x=0)=0$. $\gamma, \beta$ are constants.
Attempt Considering the b.c. (1), a solution of the following form can be assumed:
$$T(x,y)=\sum_{n=0}^{\infty}(c_3 e^{ay} + c_4 e^{-ay})\cos\bigg(\frac{n\pi x}{L}\bigg)$$
where $a=\frac{n\pi}{L}$ is the separation constant. Then we have two non-homogeneous b.c.(s) i.e. (2) and (3).
So probably these two conditions can be used one-by-one to form two simultaneous linear equations (after applying orthogonality) with two unknowns $c_3$ and $c_4$.