This is the mail archive of the gsl-discuss@sources.redhat.com mailing list for the GSL 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: hypergeometric function


Hi,
I couldn't find the discontinuty of the hypergeometry graph.
please, check the file hyperg.gif,

Yours,
Dan, Ho-Jin

Stefan Koch wrote:

>Hello,
>
>I'm a new user of the gsl-library. Everything works well. But when I tried
>to look at the hypergeometric function, I recognized a discontinuity (jump)
>in the graphic-drawing of the function. I think this should not be so. I
>used the following code segment:
>
--snip

#include <stdio.h>
#include <gsl/gsl_sf.h>

int main()
{
  int i, m=1000;
  int n=3, p=8;
  double result, x;

  for (i=0; i<m; i++) {
    x = i*0.1;
    result = gsl_sf_hyperg_1F1_int(-n, p+1, x);
    printf("%f %f\n", x, result);
  }
}


GIF image


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