2 DOUBLE PRECISION FUNCTION dlbeta (A, B)
28 DOUBLE PRECISION a, b, p, q, corr, sq2pil, d9lgmc, dgamma, dlngam,
32 DATA sq2pil / 0.9189385332 0467274178 0329736405 62 d0 /
37 IF (p .LE. 0.d0)
CALL xermsg (
'SLATEC',
'DLBETA',
38 +
'BOTH ARGUMENTS MUST BE GT ZERO', 1, 2)
40 IF (p.GE.10.d0)
GO TO 30
41 IF (q.GE.10.d0)
GO TO 20
45 dlbeta =
log(dgamma(p) * (dgamma(q)/dgamma(p+q)) )
50 20 corr = d9lgmc(q) - d9lgmc(p+q)
51 dlbeta = dlngam(p) + corr + p - p*
log(p+q)
52 1 + (q-0.5d0)*dlnrel(-p/(p+q))
57 30 corr = d9lgmc(p) + d9lgmc(q) - d9lgmc(p+q)
58 dlbeta = -0.5d0*
log(q) + sq2pil + corr + (p-0.5d0)*
log(p/(p+q))
59 1 + q*dlnrel(-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)