x=input("x =")
eps=input("precision = ")
a=0;b=1;m=1;M=%e
while M-m>eps
    if x<(a+b)/2 then
        b=(a+b)/2
        M=sqrt(m*M)
    else 
        a=(a+b)/2
        m=sqrt(m*M)
    end
end
afficher(m)
afficher(M)