Yes, the problem is solved at the moment as I only needed Bessel functions for orders ranging from 0 to about 350 and Fortran 2008 does the job. I should still look into the problem you mentioned as at some point of time I might have to come back to SLATEC for other mathematical functions. I...
Problem solved! Fortran 2008 has bessel_jn(n, x) intrinsic.
http://fortranwiki.org/fortran/show/bessel_jn
I compiled my .f95 using this function and of course setting the standard as 2008, i.e.,
gfortran -c -std=f2008 .\testf2008.f95
and then
gfortran -o a .\testf2008.o
Result...
Hello all,
I am trying to compute the Bessel function of the second kind (Bessel_y) using the SLATEC's Amos library available on Netlib. Here is the SLATEC code I use Link. Below I have pasted my test program that calls SLATEC routine CBESY...
This site uses cookies to help personalise content, tailor your experience and to keep you logged in if you register.
By continuing to use this site, you are consenting to our use of cookies.