Function wcsph_pressure 

pub fn wcsph_pressure(rho: f64, rho0: f64, b: f64, gamma: f64) -> f64

Equation of state for WCSPH: p = B * ((rho/rho0)^gamma - 1).

Typical values: gamma = 7, B = rho0 * c_s^2 / gamma.