format long x1 = (sqrt(1900^2-4)+1900)/2; x2 = 1/x1; x1 x2 %x = 0; %xx = x - (x^2-1900*x+1)/(2*x-1900); %err = 10^(-31) %while abs(xx - x) < err % x = xx; % xx = x - (x^2-1900*x+1)/(2*x-1900); %end %xx %1/xx