To my surprise, the functions g0() and g1() are defined as (C1) g0(x); - x (D1) BESSEL_I(0, X) %E (C2) g1(x); - x (D2) BESSEL_I(1, X) %E I can't find any documentation on this.