Function cubic_spline_kernel 

pub fn cubic_spline_kernel(r: f64, h: f64) -> f64

Cubic-spline SPH kernel W(r, h).

Returns the kernel value at distance r with smoothing length h.