Re: hypergeometric function

["Dan, Ho-Jin" <hjdan@sys713.kaist.ac.kr>]

[-- Attachment #1: Type: text/plain, Size: 422 bytes --]

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


[-- Attachment #2: gsl_sf_h1.c --]
[-- Type: text/plain, Size: 240 bytes --]

#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);
  }
}



[-- Attachment #3: hyperg.gif --]
[-- Type: image/gif, Size: 14044 bytes --]