function [wynik]=root(z , n)

  modul=sqrt((real(z)*real(z))+(imag(z)*imag(z)));
  kat=acos(real(z)/modul);

  i=0;
  while i<n
    wynik(i+1)=(modul^(1/n))*(cos((kat+2*i*%pi)/n)+%i*sin((kat+2*i*%pi)/n));
    i=i+1;
  end

endfunction