Solution (source code)

= Solution

If $u_1,u_2$ are two solutions with the same boundary values, put $w=u_1-u_2$. The <mean value theorem> gives
$$
\Delta w=q(x)w,
\qquad
q(x)=\int_0^1\partial_tg(x,u_2+s(u_1-u_2))\,ds\geq0.
$$
Therefore $\Delta w-qw=0$ with zero boundary data. Apply the <weak maximum principle for elliptic operators> to $w$ and $-w$ with zeroth-order coefficient $-q\leq0$. It follows that $w=0$, proving uniqueness.