File: gammai-1.m

package info (click to toggle)
octave 2.0.13.95-1
  • links: PTS
  • area: main
  • in suites: slink
  • size: 26,052 kB
  • ctags: 16,420
  • sloc: cpp: 67,184; fortran: 41,514; ansic: 26,607; sh: 7,291; makefile: 4,089; lex: 1,961; yacc: 1,852; perl: 1,676; lisp: 1,664; exp: 123
file content (6 lines) | stat: -rw-r--r-- 180 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
 
a = [.5 .5 .5 .5 .5];
x = [0 1 2 3 4];
v1 = sqrt(pi)*erf(x)./gamma(a);
v2 = gammai(a,x.*x);
v3 = gammainc(x.*x,a);
all (abs(v1 - v2) < sqrt(eps)) && all (abs(v2 - v3) < sqrt(eps))