function [wynik]=root(a,n)
  r=abs(a);
  if (imag(a) >= 0)
    fi = acos(real(a)/abs(a));
  else
    fi = -acos(real(a)/abs(a));
  end
  for i=0:1:n-1
    psi=(fi+2*%pi*i)/n;
    pierwiastki(i+1)=(r^(1/n))*(cos(psi)+%i*sin(psi)); //r^(1/n) to z0
    wynik(i+1)=pierwiastki(i+1);
  end
endfunction