%I #5 Dec 27 2013 14:10:00
%S 0,1,5,16,35,65,106,161,231,318,425,551,700,872,1069,1293,1545,1827,
%T 2140,2486,2866,3281,3734,4225,4755,5327,5942,6601,7306,8057,8856,
%U 9706,10606,11558,12565,13626,14744,15919,17154,18449,19805,21225,22709,24259,25876,27561
%N Floor(AGM(n^2, n^3)), where AGM denotes the arithmetic-geometric mean.
%o (Python)
%o import math
%o for n in range(67):
%o x = n*n
%o y = x*n
%o a = (x+y)*0.5
%o g = math.sqrt(x*y)
%o while abs(a-g)>0.00000000001:
%o x = a
%o y = g
%o a = (x+y)*0.5
%o g = math.sqrt(x*y)
%o print str(int(a))+',',
%Y Cf. A000290, A000578, A061979, A076145, A234363.
%K nonn
%O 0,3
%A _Alex Ratushnyak_, Dec 24 2013
|