format long a=1; b=-1900; c=1; f=@(x) x^2-1900*x+1; x1=(-b+sqrt(b^2-4*a*c))/(2*a); x2=1/x1; x1 x2 f(x1) f(x2)