File: BESS.FT of Tape: Various/ETH/f2
(Source file text)
FUNCTION BESS(N,Z) X=Z*Z/4. A=1. B=1. K=0 10 K=K+1 B=-X*B/(K*(N+K)) A=A+B IF (ABS(0.0001*A)-ABS(B)) 10 P=1. DO 20 I=1,N 20 P=P*I BESS=A*(Z/2)**N/P END