std.float
exp, ln, pow for f64 without libm (core::pow needs C’s pow); ~1e-15 relative
pub exp :: fn(x: f64) -> f64e^x; 0 below ~-745, inf above ~709
pub ln :: fn(x: f64) -> f64natural log; -inf for x <= 0
pub pow :: fn(x: f64, y: f64) -> f64x^y; a negative x needs a whole y, else 0