% quadform(a,b,c) returns the two roots of the quadratic equation
% ax^{2} + bx + c = 0

function [r1,r2] = quadform(a,b,c) 
  x=sqrt(b^2 - 4*a*c);
  r1=(-b-x)/(2*a);
  r2=(-b+x)/(2*a);
end
