view test/octave.test/arith/gammainc-1.m @ 4032:20d706a6740f

[project @ 2002-08-10 06:02:12 by jwe]
author jwe
date Sat, 10 Aug 2002 06:04:40 +0000
parents
children ac5ed096b496
line wrap: on
line source

a = [.5 .5 .5 .5 .5];
x = [0 1 2 3 4];
v1 = sqrt(pi)*erf(x)./gamma(a);
v3 = gammainc(x.*x,a);
all (abs(v1 - v2) < sqrt(eps))