function [x1,x2]=qroots(a,b,c)

%the following command is for showing more digits than matlab's standard 5 digits in
%the output on the screen
format long

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