At least for integer orders, Miller's backward recurrence algorithm evaluates bessel_j. This method would work for complex floats and bigfloats. For a description, see A & S. I can't think of reference to the Bessel functions with complex orders. Barton