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)) )
51 20 corr = r9lgmc(q) - r9lgmc(p+q)
52 albeta = alngam(p) + corr + p - p*
log(p+q) +
53 1 (q-0.5)*alnrel(-p/(p+q))
58 30 corr = r9lgmc(p) + r9lgmc(q) - r9lgmc(p+q)
59 albeta = -0.5*
log(q) + sq2pil + corr + (p-0.5)*
log(p/(p+q))
60 1 + q*alnrel(-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)
charNDArray max(char d, const charNDArray &m)
charNDArray min(char d, const charNDArray &m)