Function regular_bessel 

fn regular_bessel(n: i32, x: Complex) -> Complex

x^-n J_n(x), entire in x^2, or NaN where AMOS fails.