def Approx(seuil): n=0, u=0 while 1/4^n>seuil: u=1/(2+u) n=n+1 return u