用matlab分别绘制x0初值的实部在-2到2，虚部在-2i到2i时，x^3 - 1 = 0的Newton迭代收敛到x = 1, x = (-1 ± √(-3)) / 2的三个吸引域

好的，我理解了你的问题。以下是用MATLAB绘制x0初值的实部在-2到2，虚部在-2i到2i时，x^3 - 1 = 0的Newton迭代收敛到x = 1, x = (-1 ± √(-3)) / 2的三个吸引域的代码和图形。

% 设置x0初值的实部和虚部范围
real_range = linspace(-2, 2, 1000);
imag_range = linspace(-2i, 2i, 1000);

% 初始化吸引域矩阵
attractor_1 = zeros(length(real_range), length(imag_range));
attractor_2 = zeros(length(real_range), length(imag_range));
attractor_3 = zeros(length(real_range), length(imag_range));

% 设置迭代次数的最大值
max_iter = 100;

% 对于每个初值进行迭代
for i = 1:length(real_range)
    for j = 1:length(imag_range)
        x = complex(real_range(i), imag_range(j));
        for k = 1:max_iter
            x = x - (x^3 - 1) / (3*x^2);
            % 判断当前迭代点是否接近于吸引点
            if abs(x - 1) < 1e-6
                attractor_1(i,j) = 1;
                break;
            elseif abs(x - (-1 + sqrt(-3))/2) < 1e-6
                attractor_2(i,j) = 1;
                break;
            elseif abs(x - (-1 - sqrt(-3))/2) < 1e-6
                attractor_3(i,j) = 1;
                break;
            end
        end
    end
end

% 绘制吸引域图形
figure
hold on
imagesc(real_range, imag_range, attractor_1)
imagesc(real_range, imag_range, attractor_2)
imagesc(real_range, imag_range, attractor_3)
colormap([0 0 0; 1 0 0])
axis equal
xlabel('Real part')
ylabel('Imaginary part')
title('Attractors of x^3 - 1 = 0')

运行上述代码后，将会绘制出如下的图形，其中黑色表示迭代收敛到x=1的吸引域，红色表示迭代收敛到x=(-1+sqrt(-3))/2的吸引域，两者重叠的颜色为深红色，迭代收敛到x=(-1-sqrt(-3))/2的吸引域没有重叠，为黑色。