# 
restart;
phi:=x->1+1/(1+x);
r0:=1;
r1:=phi(r0);
r2:=phi(r1);
a[0]:=1;
for n from 0 to 10 do
    a[n+1]:=phi(a[n]);
od;

for n from 0 to 11 do
    a[n]^2;
od;
a[0]:=1;
for n from 0 to 100 do
    a[n+1]:=phi(a[n]);
od:
a[101];
a[101]^2;
1/2*a[101]^2;

#