x = input('enter your guess')
for ii = [1:1000]
    x=x-((x*x-2)/(2*x));
    if mod(ii,100)==0
        fprintf('%e\n',x)
    end
end