a=1;
b=3;
c=1;
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)
endif