This is the mail archive of the fortran@gcc.gnu.org mailing list for the GNU Fortran project.

Index Nav:	[Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav:	[Date Prev] [Date Next]	[Thread Prev] [Thread Next]
Other format:	[Raw text]

Re: Problems with GAMMA functions?

From: Jerry DeLisle <jvdelisle at verizon dot net>
To: Thomas Koenig <tkoenig at netcologne dot de>
Cc: Fortran List <fortran at gcc dot gnu dot org>
Date: 2007年10月08日 23:18:54 -0700
Subject: Re: Problems with GAMMA functions?
References: <Pine.CYG.4.58.0710062310490.724@homepc> <1191709485.3259.6.camel@meiner.onlinehome.de> <Pine.OSF.4.58.0710081324410.25589@ax0rm1.roma1.infn.it> <1191908978.3292.3.camel@meiner.onlinehome.de> <470B1A42.104@verizon.net> <470B1AAD.4010006@verizon.net>

Heres some test data. Its failing in the do loop for two values at n=18 and 19:

n xs ts abs(gamma(xs)-ts)/ts

 1 1.5000000 0.88622695 0.0000000
 2 2.5000000 1.3293403 0.0000000
 3 3.5000000 3.3233509 7.17404163E-08
 4 4.5000000 11.631728 2.45967129E-07
 5 5.5000000 52.342777 2.18637453E-07
 6 6.5000000 287.88528 2.12012068E-07
 7 7.5000000 1871.2543 0.0000000
 8 8.5000000 14034.407 4.17500701E-07
 9 9.5000000 119292.46 4.58432169E-07
10 10.500000 1133278.4 8.82395739E-07
11 11.500000 11899423. 1.26056534E-06
12 12.500000 1.36843360E+08 8.18454055E-07
13 13.500000 1.71054208E+09 3.74150403E-07
14 14.500000 2.30923182E+10 1.50768756E-06
15 15.500000 3.34838596E+11 1.56579324E-06
16 16.500000 5.18999822E+12 1.31324589E-06
17 17.500000 8.56349706E+13 8.81619656E-07
18 18.500000 1.49861207E+15 2.14947249E-06
19 19.500000 2.77243231E+16 2.55612940E-06
20 20.500000 5.40624301E+17 0.0000000

Regards,

Jerry

Follow-Ups:
- Re: Problems with GAMMA functions?
  - From: Thomas Koenig

References:
- Problems with GAMMA functions?
  - From: Angelo Graziosi
- Re: Problems with GAMMA functions?
  - From: Thomas Koenig
- Re: Problems with GAMMA functions?
  - From: Angelo Graziosi
- Re: Problems with GAMMA functions?
  - From: Thomas Koenig
- Re: Problems with GAMMA functions?
  - From: Jerry DeLisle
- Re: Problems with GAMMA functions?
  - From: Jerry DeLisle

Index Nav:	[Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav:	[Date Prev] [Date Next]	[Thread Prev] [Thread Next]