Mercurial > hg > octave-lyh

comparison scripts/specfun/beta.m @ 2312:204cc7db6f4a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 1996-07-11 21:20:36 by jwe]
author jwe
date Thu, 11 Jul 1996 21:20:36 +0000
parents 2b5788792cad
children 5ca126254d15
comparison
equal deleted inserted replaced
2311:2b5788792cad 2312:204cc7db6f4a
20 ## usage: beta (a, b) 20 ## usage: beta (a, b)
21 ## 21 ##
22 ## Returns the beta function beta(a,b) = gamma(a) * gamma(b) / gamma(a+b) 22 ## Returns the beta function beta(a,b) = gamma(a) * gamma(b) / gamma(a+b)
23 ## of a and b. 23 ## of a and b.
24 24
25 ## Author: KH <Kurt.Hornik@ci.tuwien.ac.at>
26 ## Created: 13 June 1993
27 ## Adapted-By: jwe
28
25 function retval = beta (a, b) 29 function retval = beta (a, b)
26 30
27 ## Written by KH (Kurt.Hornik@ci.tuwien.ac.at) on Jun 13, 1993
28
29 if (nargin != 2) 31 if (nargin != 2)
30 usage ("beta (a, b)"); 32 usage ("beta (a, b)");
31 endif 33 endif
32 34
33 retval = exp (lgamma (a) + lgamma (b) - lgamma (a+b)); 35 retval = exp (lgamma (a) + lgamma (b) - lgamma (a+b));