function [x1,x2]=myrootc( a,b,c)
diskr=b*b-4*a*c;
x1=(-b+sqrt(diskr))/(2.*a);
x2=(-b-sqrt(diskr))/(2.*a);
end

