function roots = nikola_petrov_quadr_eqn1 ( a, b, c )

% nikola_petrov_quadr_eqn1   A naive solver of a quadratic equation
%      
%      A code that finds the root(s) of the quadratic equation
%                a x^2 + b x + c = 0 
%      using the formula
%                x = ( -b +- sqrt(b^2-4ac) ) / (2a)
%

format long e

discr = b^2 - 4*a*c;
roots = [ (-b-sqrt(discr))/2/a ; (-b+sqrt(discr))/2/a ];
