read AsyRec:


lu:=AsyBin(binomial(n,k)^2*binomial(n+k,k)^2,n,k,5,10000):

print(`The asymtotics for the Apery Sequence`):
print(Sum(binomial(n,k)^2*binomial(n+k,k)^2,k=0..n)):
print(`equals the following`):
print(lu[1]*lu[2]):
print():
print(`Note: The constant in front is a non-rigorous estimate.`):
print(`This took `, time(), `seconds.`):
quit: