Mercurial > hg > octave-nkf
view test/octave.test/arith/gamma-1.m @ 7217:c8abc18322b7
[project @ 2007-11-29 20:50:24 by jwe]
author | jwe |
---|---|
date | Thu, 29 Nov 2007 20:50:24 +0000 |
parents | 7050b5006ef9 |
children |
line wrap: on
line source
x = [.5, 1, 1.5, 2, 3, 4, 5]; v = [sqrt(pi), 1, .5*sqrt(pi), 1, 2, 6, 24]; all(abs(gamma(x) - v) < sqrt(eps)) && all(abs(lgamma(x) - log(v)) < sqrt(eps))