gsl version: 1.7 computer system: PowerMac G4 operating system: Mac OS X 10.3.9 compiler: gcc 3.3
---------- Forwarded message ---------- From: Yoram Burak <[EMAIL PROTECTED]> Date: Dec 23, 2005 4:23 PM Subject: Error in evaluating scaled regular modified bessel functions of large arguments To: [email protected] Attempt to evaluate gsl_sf_bessel_In_scaled(n,x) results in an error for n>=2 and x larger than about 10^7. #include <stdio.h> #include <gsl/gsl_sf_bessel.h> int main() { double logx, x; for(logx = 0; logx < 20; logx++) { x = exp(logx); printf("%13le %13le\n", x, gsl_sf_bessel_In_scaled(2,x)); } } -- Yoram Burak Kavli Institute for Theoretical Physics UC Santa Barbara, California 93106 -- Yoram Burak Kavli Institute for Theoretical Physics UC Santa Barbara, California 93106 Voice: +1 805 893-7310 Fax: +1 805 893-2431 E-mail: [EMAIL PROTECTED] _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
