12345678910111213141516171819202122232425 |
- import mpmath as mp
- def LeRu_cutoff(z):
- x = mp.fabs(z)
- return int(x + 11 * x**(1/3) + 1)
- def psi(n,z):
- return mp.sqrt( (mp.pi * z)/2 ) * mp.autoprec(mp.besselj)(n+1/2,z)
- def r(n,z):
- if n > 0:
- return psi(n-1,z)/psi(n,z)
- return mp.cos(z)/mp.sin(z)
- def D1(n,z):
- return r(n,z) - n/z
|