Ya vimos que la discretización de las ecuaciones de conservación en forma Lagrangiana para SPH quedaba: $latex frac{d}{dt}boldsymbol{v}_a = – sum_b m_b (frac{P_a}{rho_a^2} + frac{P_b}{rho_b^2}) nabla_a W{ab} = boldsymbol{F}_a$ $latex frac{d}{dt} e_a = frac{1}{2} sum_b m_b (frac{P_a}{rho_a^2} + frac{P_b}{rho_b^2}) boldsymbol{v}_{ab}cdot nabla_a W{ab} = E_a$ con $latex rho_a = sum_b m_b W_{ab}$ y $latex P_a$ de …