de Moivres Nth Root Function on returns 2 out 5 roots? Spot the fault?

Question

FortuitousMonkey 2018-3-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/388376-de-moivres-nth-root-function-on-returns-2-out-5-roots-spot-the-fault

评论： FortuitousMonkey 2018-3-15

Hi,

I wonder if anyone could help me find the fault in the code? When I call the function using the below it only returns 2 values, the first and last root, but misses 3 in the middle, subsequently my plot is also incorrect (it should return 5 equally spaced points). I can't find any error in the code.

>> complex_root1(2+3j,5)
D =
   1.2675 + 0.2524i   0.6317 - 1.1275i

My function is

function res=complex_root(z,N)
%z is the complex number
%n is the number of roots, must be positive and larger than o
k=[0,N-1];
%Logic statement if N above 0 proceed with calculation
if N>0
  a=complex(z);
  magnitude=abs(a);
  r=(abs(a));
  t=angle(a);
  %de Moivre's forumla 
  D=r^(1/N)*((cos((t+2*pi*k)/N))+((1i*sin((t+2*pi*k)/N))))
  polarplot (D)
%Logical statement, if N below zero return error message
elseif N<=0;
  disp 'N must be a positive integer above 0'
end
end