macro
Beta B p q

#Copyright Mathews and Malnar, 12Nov01
#calculates complete beta function:
#B(p,q)=gamma(p+q)/gamma(p)/gamma(q) where gamma(n)=(n-1)!

#Example call:
#   mtb> %beta k1 4 12
#Answer is B = 5460

mconstant B p q i

If p>q
   let B=q
   let p=p-1
   do i=1:p
      let B=B*(q+i)/i
   enddo
   let p=p+1
else
   let B=p
   let q=q-1
   do i=1:q
      let B=B*(p+i)/i
   enddo
   let q=q+1
endif

print p q B
endmacro