QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Mathematica check the incomplete gamma function results. 

Former-commit-id: 0eaac1b86051acdffd08f2861f4470a9030d2eed
This commit is contained in:
Marek Nečada 2018-08-20 16:59:39 +03:00
parent 7093ff3add
commit 2f6ab47d2e
2 changed files with 14 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
tests/gsl-intgamma/incgamma.mathematica Normal file
View File

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

4
tests/gsl-intgamma/quadratures.c
View File

@ -7,10 +7,10 @@
#define EPSABS 0 #define EPSABS 0
#endif #endif
#ifndef EPSREL #ifndef EPSREL
#define EPSREL 1e-8 #define EPSREL 1e-13
#endif #endif
#ifndef LIMIT #ifndef LIMIT
#define LIMIT 8192 //??? #define LIMIT 30000 //???
#endif #endif
#ifndef R0 #ifndef R0
#define R0 8e-6 #define R0 8e-6