Euler:

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

Conservation:

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

Poisson:

(nabla)^2 *Phi_0+(nabla)^2 *delta*Phi=4*pi*G*rho_0+4*pi*G*delta*rho