pub fn hyperg_U(a: f64, b: f64, x: f64) -> f64
This routine computes the confluent hypergeometric function U(a,b,x).