Euler:

dtemps(u;1)+2*(dtemps(a;1)/a)*u=((nabla)*delta*Phi)/a^2-((nabla)*delta*P)/(rho_0*a^2)

Conservation:

dtemps(lambda;1)=-pscalaire(nabla;u)