re and im N razy plain

% script m-file: ex1.m
clc; clear; close all;

N=6;
for t=1:N

num=rand()+rand()*1i;
%num = 7+9i;
n = 10;
k = 0:1:n;
rt = abs(num).^(1/n)*exp(1i*(2*pi*k+angle(num))/n);
%rr = roots([1, zeros(1,n-1), -num])

hm = abs(max(real([rt,num]))-min(real([rt,num]))).*0.5; %marginesy
vm = abs(max(imag([rt,num]))-min(imag([rt,num]))).*0.5;
hrange = [min(real([rt,num]))-hm, max(real([rt,num]))+hm];
vrange = [min(imag([rt,num]))-vm, max(imag([rt,num]))+vm];
figure(1);
plot([0, 0], vrange, 'k'); hold on; grid on;
plot(hrange, [0, 0], 'k');
axis([hrange, vrange]);

plot(real(num), imag(num), 'ro');
plot(real(rt), imag(rt), 'bo'); hold off;
pause;
end