function [y]=Resoudre(x)

if x(1) ~= 0 then

  delta = (x(2))^2 - 4*x(1)*x(3);
  if delta > 0 then
    y(1) = (-x(2)+sqrt(delta))/(2*x(1));
    y(2) = (-x(2)-sqrt(delta))/(2*x(1));
  elseif delta == 0 then
    y(1) = (-x(2))/(2*x(1));
    y(2) = (-x(2))/(2*x(1));
  else
    y(1) = (-x(2)+%i*sqrt(-delta))/(2*x(1));
    y(2) = (-x(2)-%i*sqrt(-delta))/(2*x(1));
  end;
 else
	y(1)=0;
	y(2)=0;
 end;
endfunction 

b=[0,0]
a=[1,-10^5,-1];
b=Resoudre(a)


