####
format long    % show all the digits
a = 1;  b = -(1e6+1e-6);  c = 1;
%% 
% First, we find the "good'' root using the quadratic forumla. 
%%
x1 = (-b + sqrt(b^2-4*a*c)) / (2*a);
%% 
% Then we use the better formula for computing the other root. 
%%
x2 = c/(a*x1)
##### SOURCE END #####
--></body></html>