function [r1,r2]=quad_rootfinder(a,b,c)
% 
%function [r1,r2]=root_finder(a,b,c)
%
%This function computes the roots of 
% the quadratic equation 
% ax^2 + bx + c =0
%
% To run type 
% [r1,r2]=quad_rootfinder(1,-4,4)
%
sqrt_d = sqrt(b^2 - 4*a*c) ;
r1 = (-b - sqrt_d)/(2*a) ;
r2 = (-b + sqrt_d)/(2*a) ;