[x1,x2,diskr]=myroot(3,3,1)

function [x1,x2,diskr]=myroot( a,b,c)
diskr=b*b-4*a*c;
if(diskr>0)
x1=(-b+sqrt(diskr))/(2.*a);
x2=(-b-sqrt(diskr))/(2.*a);
return
end
if(diskr==0)
x1=-b/(2.*a);
return
end
x1=-b/(2.*a);
x2=(sqrt(-diskr))/(2.*a);
return
end

