34 DATA sq2pil / 0.9189385332 0467274 e0 /
39 IF (p .LE. 0.0) CALL
xermsg(
'SLATEC',
'ALBETA',
40 +
'BOTH ARGUMENTS MUST BE GT ZERO', 1, 2)
41 IF (p.GE.10.0) go
to 30
42 IF (q.GE.10.0) go
to 20
46 albeta =
log(gamma(p) * (gamma(q)/gamma(p+q)) )
59 albeta = -0.5*
log(q) + sq2pil + corr + (p-0.5)*
log(p/(p+q))
OCTAVE_EXPORT octave_value_list or N dimensional array whose elements are all equal to the base of natural logarithms The constant ex $e satisfies the equation log(e)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
charNDArray max(char d, const charNDArray &m)
subroutine xermsg(LIBRAR, SUBROU, MESSG, NERR, LEVEL)
charNDArray min(char d, const charNDArray &m)