The old bessel(arg, order) function has been deprecated for quite some time now. All functionality that was in bessel was always in bessel_j. so there's no loss. (Except the args for bessel were reverse of those for bessel_j and friends.) I think it's time for bessel to go. Any objections? Ray