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#L...
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