∇²u = -ρ/ε