OFFSET
1,4
FORMULA
a(n) = ceiling(sqrt(Pi - Sum_{i=1..n-1} 1/a(i)^2)).
EXAMPLE
a(4)=3 because Pi - 1/a(1)^2 - 1/a(2)^2 - 1/a(3)^2 = Pi - 1 - 1 - 1 = 0.1415926... and 3 is the smallest integer such that 1/3^2 = 0.1111111... < 0.1415926...
PROG
(PARI) f(x)=ceil(sqrt(1/x))
lista(n)={my(k=Pi, v=vector(n)); for(T=1, n, v[T]= f(k); k-=1/v[T]^2); v}
CROSSREFS
KEYWORD
nonn
AUTHOR
Hauke Worpel (hw1(AT)email.com), Nov 11 2006
STATUS
approved