[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] [bug #45926] error in generalized hypergeometric function
From: |
Patrick Alken |
Subject: |
[Bug-gsl] [bug #45926] error in generalized hypergeometric function |
Date: |
Fri, 11 Sep 2015 10:09:02 +0000 |
User-agent: |
Mozilla/5.0 (X11; Ubuntu; Linux x86_64; rv:40.0) Gecko/20100101 Firefox/40.0 |
URL:
<http://savannah.gnu.org/bugs/?45926>
Summary: error in generalized hypergeometric function
Project: GNU Scientific Library
Submitted by: psa
Submitted on: Fri 11 Sep 2015 10:09:01 AM GMT
Category: Runtime error
Severity: 3 - Normal
Operating System:
Status: None
Assigned to: None
Open/Closed: Open
Release:
Discussion Lock: Any
_______________________________________________________
Details:
from aim =at= stats =dot= uwo =dot= ca
The bug I am reporting arose when I tried to the gsl package in R:
gsl: wrapper for the Gnu Scientific Library
The authors of the gsl package have determined that the error is in the Gnu
Scientific Library.
The problem arises with the generalized hypergeometric function and is
illustrated
in the code below. Please excuse my odd syntax where I used 1+9 etc instead of
10
etc. Both Mathematica and MatLab get the result correct.
Browse[2]> hyperg_2F1(-0.2,-0.2 + 9,1 + 9,0.8)
[1] 0.7799897
Browse[2]> hyperg_2F1(-0.2,-0.2 + 10,1 + 10,0.8)
[1] NaN
MATHEMATICA GETS IT CORRECT
Hypergeometric2F1[-0.2, -0.2+9, 1+9, 0.8]
0.77999
Hypergeometric2F1[-0.2, -0.2+10, 1+10, 0.8]
0.775746
MATLAB GETS IT CORRECT
>> hypergeom([-0.2,-0.2+9],1+9,0.8)
ans =
0.7800
>> hypergeom([-0.2,-0.2+10],1+10,0.8)
ans =
0.7757
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?45926>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Bug-gsl] [bug #45926] error in generalized hypergeometric function,
Patrick Alken <=