We have $\Delta u=f$ in $D$, and $\dfrac{\partial u}{\partial n}+au=h$ on boundary of D, where $D$ is a domain in three dimension and $a$ is a positive constant. $\dfrac{\partial u}{\partial n}=\triangledown u\cdot n$ ($n$ is normal vector).
My thoughts: Suppose there are $u_1$ and $u_2$, satisfis the above equations. Let $w=u_1-u_2$, then we have $\Delta w=0$ in $D$, and $\dfrac{\partial w}{\partial n}=-aw$ on boundary of D. Maximum modulus principle may be useful but I don't know where to put it in. And energy method seems not helpful in this question.
Any help would be appreciated!