Trouble with bessik


statmobile
11-19-2003, 11:38 PM
Hey Pplz,
I'm trying to use bessik.c, but it keeps telling me that I have bad arguments in bessik. I'm using the sample code xbessik.c, but that isn't working either (nr.h has a previous declaration of 'select'). Can someone please explain how one can get at a Modified Bessel Function of the third kind with this code? I've been using the code:

bessik(2.0,1.0,&xri,&xrk,&xrip,&xrkp);

Saul Teukolsky
11-20-2003, 08:19 AM
Dear Stat,

1) The "select" problem is discussed at www.nr.com, click on "Problems", click on "common problems". The simplest fix is to comment out the select declaration in nr.h, if you're not using that recipe.

2) To call bessik, make sure all your arguments are declared float. Declare float variables instead of using 1.0 and 2.0 as arguments, depending on your compiler options.

Saul Teukolsky

statmobile
11-20-2003, 09:53 AM
That did the trick, thanks for your help