function newton(n)
% n %количество шагов
dx = 3/n; nits = 10; tol = 0.25; w1 = 1; w2 = exp(2*pi*i/3);
w3=w2*w2;
%Генерация матрицы...
[X,Y] = meshgrid(-1.5:dx:1.5); m = X + i*Y;
%Матрица размера int2str(prod(size(m))) создана
%выполняются int2str(nits) шагов
for k=1:nits
    m = m - (m.^3-1)./(3*m.^2);
end;
% mm записываем самый близкий корень
mm = 1*(abs(m-w1)<tol) + 2*(abs(m-w2)<tol) + 3*(abs(m-w3)<tol);
% the errs hold which entries are not converged yet
errs = find(mm==0); nerr = prod(size(errs));
xtra = 0;
%теперь мы заканчиваем вычислять те вхождения, которые должны все же сходиться
[jj,kk] = size(mm); for j = errs'
    x = m(j);
    xx = 1*(abs(x-w1)<tol)+2*(abs(x-w2)<tol)+3*(abs(x-w3)<tol);
    while xx==0
        x = x - (x*x*x-1) / (3*x*x);
        xtra = xtra+1;
        xx=1*(abs(x-w1)<tol)+2*(abs(x-w2)<tol)+3*(abs(x-w3)<tol);
    end; % end while
    m(j) = x;
    mm(j) = xx;
end;    % end for loops and if mm==0
% выводим рисунок на экран
image(mm); axis('image'); axis('off');
title('Множество Ньютона')
colormap(eye(3));   % for color monitors round
clear n tol m mm x xx l nits w1 w2 w3 b c dx pict errs j jj k kk nerr xtra Y