Constant RHO_SCALE 

const RHO_SCALE: f64 = 0.9999;

ρ = RHO_SCALE·tanh(x) keeps ρ inside the formulas’ bound.