%
%
 Square - equation
%
% ax^2+bx+c=0
%
a=input('a ? ');
b=input('b ? ');
c=input('c ? ');

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

if(discr<0)
   disp([' Two complex solutions']);
   disp(['x1:' num2str(x1) '   x2:' num2str(x2)]);
else
     if(discr>0)
      disp([' Two real solutions']);
      disp(['x1:' num2str(x1) '   x2:' num2str(x2)]);
     else
      disp([' One real solutions']);
      disp(['x:' num2str(x1) ]);
   end
end