format long e

a=1890
b=1910
x2=a+(b-a)/2;
f = @(x) x^2 - 1900*x + 1;
while (abs(f(x2)) > 10^(-14))
  x2=a+(b-a)/2;
  if(f(a)*f(x2)>0)a=x2;
  else b=x2;
  end
endwhile
x2

a=0
b=1
x1=a+(b-a)/2;
f = @(x) x^2 - 1900*x + 1;
while (abs(f(x1)) > 10^(-14))
  x1=a+(b-a)/2;
  if(f(a)*f(x1)>0)a=x1;
  else b=x1;
  end
endwhile
x1