PDE with Integral constraint

I am trying to solve the Non-linear Schrodinger equation

$-\Delta \psi(r) + \psi(r) - |\psi(r)|^2\psi(r) = 0$ where $r \in \Omega$

In a square domain ($(x,y) \in \Omega$ where $\Omega=[0,1]\times [0,1]$) with the Dirichlet condition

$\psi(r) = 0$ for $r \in \partial \Omega$

$\int_\Omega d^2r \, |\psi(r)|^2 = 1$

Does somebody know how to include the integral condition into NDSolve?

• Hola Alvaro, welcome to Mathematica.SE, please consider taking the tour so you learn the basic rules of the site. Here its considered helpful to share your code attempts, in a well formatted easy to copy&paste form. Don't erase your nicely formatted $\LaTeX$, just add Mathematica code in addition. – rhermans Oct 12 '14 at 16:20
• Hello, actually my first attempt was using a Lagrange multiplier, but NDSolve accepts only algebraic constraints, as for example: $f(r)+g(r)=1$. To put it differently it seems to me that it only accepts local constraints and not non-local as $\int_\Omega d^2 r \, |\psi(r)|^2 = 1$. In fact my question is precisely how to circumvent this problem. Thanks – Alvaro Rojo Bravo Oct 13 '14 at 13:54