f[x_]:=x+1/2(a-x^2);

odm = N[(1/2)^(1/2),20]

x1=0;a=1/2;n=0;

x1=f[x1];
 N[x1,20]
odm
n=n+1