/* B(p,q) */
new; cls;
print B(1/2,1/2);



proc B(p,q);
   retp(gamma(p).*gamma(q)./gamma(p+q));
endp;