a = 1.0
pred = 1.0
syms x
a = int(x/(1+x), 0, 1)
pred = vpa(a, 9)/3 + 1
for i = 2:4
pred = pred*(vpa(a,9)/3)+1
end