function [x1,x2] = quadform(a,b,c)
%QUADFORM Uses the quadratic formula to 
% find the two roots of a*x^2 + b*x + c

d=sqrt(b^2 - 4*a*c);  %discriminant

x1 = (-b + d)/(2*a); 
x2 = (-b - d)/(2*a);

end

