bessel_j[spherical](n, z) computes \(j_n(z)\) where the plain call computes \(J_n(z)\). Accepted by the four Bessel families.