function roots = quadratic_equation(a,b,c)

%Calculate the roots
root1 = (-b + sqrt(b^2 - 4*a*c))/(2*a);
root2 = (-b - sqrt(b^2 - 4*a*c))/(2*a);

roots = [root1;
         root2];