format long 
x1 = (1900+sqrt(1900^2-4))/2
x2 = 1/x1