function [x1 x2] = quadform1(a,b,c)
% This function solves a quadratic equation

d = sqrt(b^2 - 4*a*c);
x1 = (-b + d) / (2*a);
x2 = (-b - d) / (2*a);