a=1;
b=2;
c=3;
delta=b^2-4*a*c;
if delta>0
x1=(-b-sqrt(b^2-4*a*c))/(2*a)
x2=(-b+sqrt(b^2-4*a*c))/(2*a)
else
delta=sqrt(-b^2-4*a*c);
x1r=-b/(2*a);
x2r=-b/(2*a);
x1i=(b^2-4*a*c)/(2*a);
x2i=-(b^2-4*a*c)/(2*a);
x1r+i*x1i
x2r+i*x2i 
endif