Проблема с уравнением Пуассона в Mathematica

[Deneira]

Есть уравнение Пуассона $\triangle u =1$ со следующими краевыми условиями: $u(0,y)=u(4,y)=C (\operatorname{const}), \frac{\partial u}{\partial y}(x,0)=\frac{\partial u}{\partial y}(x,2)=0; x\in[0,4], y\in[0,2]$.
При попытке решить данное уравнение в Mathematica 10 возникли проблемы с этим условием: $\frac{\partial u}{\partial y}(x,0)=0$. Подскажите, пожалуйста, как исправить.

NDSolve[{D[u[x, y], x, x] + D[u[x, y], y, y] == 1, u[0, y] == u[4, y] == C1, (D[u[x, y], y] /. y -> 0) == 0, (D[u[x, y], y] /. y -> 2) == 0}, u[x, y], {x, 0, 4}, {y, 0, 2}]

[Vince Diesel]

Для численного решения вместо постоянной $C$ должно быть какое-то конкретное число. Во-вторых, в 10й версии граничные условия задаются с помощью DirichletCondition и NeumannValue, примеры есть в справке по этим командам.