def make_factorial2_vals():
from mpmath import fac2
n = range(150)
fac2 = [fac2(val) for val in n]
return make_special_vals('factorial2_vals', ('n', n), ('fac2', fac2))
def make_factorial2_vals():
from mpmath import fac2
n = list(range(150))
fac2 = [fac2(val) for val in n]
return make_special_vals('factorial2_vals', ('n', n), ('fac2', fac2))
def sph_jn_power(n, z, terms=100):
zm = mpmathify(z)
s = sum((-zm**2/2)**k/(factorial(k) * fac2(2*n + 2*k + 1)) for k in xrange(terms))
return zm**n * s
def sph_jn_power(n, z, terms=100):
zm = mpmathify(z)
s = sum((-zm**2 / 2)**k / (factorial(k) * fac2(2 * n + 2 * k + 1))
for k in xrange(terms))
return zm**n * s
