Thanks Pauli

But I have already very shortly built  for bessel function, but the code you gave me is in Fortran.. I also used f2py but I could not manage to read fortran codes..that is why I have asked in Python what is wrong??


Пятница, 21 декабря 2012, 12:46 UTC от Pauli Virtanen <pav@iki.fi>:
Happyman <bahtiyor_zohidov <at> mail.ru> writes:
[clip]
> IF I GIVE ( it is necessary value for my program ):
> a , b = sph_jn ( 536 , 2513.2741228718346 + 201.0619298974676j )

The implementation of the spherical Bessel functions is through
this Fortran code:

https://github.com/scipy/scipy/blob/master/scipy/special/specfun/specfun.f#L1091

It does not have asymptotic expansions for dealing with parts of
the complex plane where the computation via the recurrence does not
work.

--
Pauli Virtanen


_______________________________________________
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion