def sph_inkn(n,z): """Compute the spherical Bessel functions, in(z) and kn(z) and their derivatives for all orders up to and including n. """ if not (isscalar(n) and isscalar(z)): raise ValueError("arguments must be scalars.") if (n!= floor(n)) or (n<0): raise ValueError("n must be a non-negative integer.") if iscomplex(z) or less(z,0): nm,In,Inp,kn,knp = specfun.csphik(n,z) else: nm,In,Inp = specfun.sphi(n,z) nm,kn,knp = specfun.sphk(n,z) return In,Inp,kn,knp
def sph_inkn(n, z): """Compute the spherical Bessel functions, in(z) and kn(z) and their derivatives for all orders up to and including n. """ if not (isscalar(n) and isscalar(z)): raise ValueError, "arguments must be scalars." if (n != floor(n)) or (n < 0): raise ValueError, "n must be a non-negative integer." if iscomplex(z) or less(z, 0): nm, In, Inp, kn, knp = specfun.csphik(n, z) else: nm, In, Inp = specfun.sphi(n, z) nm, kn, knp = specfun.sphk(n, z) return In, Inp, kn, knp
def sph_in(n,z): """Compute the spherical Bessel function in(z) and its derivative for all orders up to and including n. """ if not (isscalar(n) and isscalar(z)): raise ValueError("arguments must be scalars.") if (n!= floor(n)) or (n<0): raise ValueError("n must be a non-negative integer.") if (n < 1): n1 = 1 else: n1 = n if iscomplex(z): nm,In,Inp,kn,knp = specfun.csphik(n1,z) else: nm,In,Inp = specfun.sphi(n1,z) return In[:(n+1)], Inp[:(n+1)]
def sph_in(n, z): """Compute the spherical Bessel function in(z) and its derivative for all orders up to and including n. """ if not (isscalar(n) and isscalar(z)): raise ValueError, "arguments must be scalars." if (n != floor(n)) or (n < 0): raise ValueError, "n must be a non-negative integer." if (n < 1): n1 = 1 else: n1 = n if iscomplex(z): nm, In, Inp, kn, knp = specfun.csphik(n1, z) else: nm, In, Inp = specfun.sphi(n1, z) return In[:(n + 1)], Inp[:(n + 1)]