Function sinh_ratio 

fn sinh_ratio(sigma: f64, x: f64, h: f64) -> f64

sinh(σx) / sinh(σh) for 0 ≤ x ≤ h, written so a large σh cannot overflow: both exponentials are of non-positive argument.