oldtests/gsl-intgamma/incgamma.mathematica

   1 l:= ReadList["incgamma.in", Table[Number, {2}]];
   2 For[i = 1, i <= Length[l], i++,
   3   p := l[[i]];
   4   j := p[[1]];
   5   x := p[[2]];
   6   res := NIntegrate[SigmaIntegrand[n, k, r, x], {x,eta,Infinity}];
   7   Print[j, " ", CForm[x], " | ", CForm[N[Gamma[1/2-j,x],16]],
   8                           " | ", CForm[N[Gamma[1/2-j,x*I],16]],
   9                           " | ", CForm[N[Gamma[1/2-j,-x],16]],
  10                           " | ", CForm[N[Gamma[1/2-j,-x*I],16]]]
  11 ];
  12