マンデルブロ集合を書くプログラム

% パラメータ設定
maxIter = 5000; % 最大反復回数
xlim = [-2, 1]; % x範囲
ylim = [-1.5, 1.5]; % y範囲
resolution = 1000; % 解像度

% 複素数平面のメッシュグリッドを作成
[x, y] = meshgrid(linspace(xlim(1), xlim(2), resolution), linspace(ylim(1), ylim(2), resolution));
c = x + 1i*y;
z = zeros(size(c));
mandelbrotSet = zeros(size(c));

% マンデルブロ集合の計算
for k = 1:maxIter
    z = z.^2 + c;
    mandelbrotSet = mandelbrotSet + (abs(z) < 2);
end

% 結果をプロット
imagesc(xlim, ylim, mandelbrotSet);
colormap(hot);
axis off;
title('Mandelbrot Set');