Euler :

drond*delta*v/drondt+(pscalaire(v_0;nabla))*delta*v+pscalaire(delta*v;nabla)*v_0=-(nabladelta*P)/rho_0-nabladelta*Phi

Conservation:

drond*delta/drondt+(pscalaire(v_0;nabla))*lambda=-pscalaire(nabla;delta*v)

Poisson:

(nabla)^2 *delta*Phi=4*pi*G*rho_0*lambda