read zeilWZP:
 
SUMMAND:=k*convert(binomial(n,k)^2*binomial(n+k,k)^2,factorial):
 
F:=normal((2*n^2+2*n+k+1)/(n+k+1)/(k+1)/(-n+k)+1/2/(k+1)-1/2/k):
G:=2*(n+1)/(n-k+1)/(n+k+1):
 
zeilWZP(SUMMAND,F,G,k,n,N):
 
quit: