帮我找一段绘制三维mandelbrot集的Matlab代码

以下是一个绘制三维Mandelbrot集的Matlab代码：

% 设置参数
N = 100;
x = linspace(-2, 2, N);
y = linspace(-2, 2, N);
z = linspace(-2, 2, N);
[X,Y,Z] = meshgrid(x,y,z);
C = X + 1i*Y + Z*1i;

% 计算Mandelbrot集
Z0 = zeros(N,N,N);
M = zeros(N,N,N);
for k = 1:20
    Z0 = Z0.^2 + C;
    M(abs(Z0) > 2 & M == 0) = k;
end

% 绘制三维图像
figure
isosurface(X,Y,Z,M,10)
axis equal
xlabel('Real axis')
ylabel('Imaginary axis')
zlabel('Z axis')

这段代码将会生成一个三维的Mandelbrot集图像。你可以根据需要更改参数N，以获得更高的分辨率图像。