在MATLAB中生成Barnsley Fern分形

import matlab.graphics.*;

figure;
axis off;
hold on;

function [x, y] = barnsley_fern(x, y)
    r = rand();
    if r < 0.01
        x_new = 0;
        y_new = 0.16 * y;
    elseif r < 0.86
        x_new = 0.85 * x + 0.04 * y;
        y_new = -0.04 * x + 0.85 * y + 1.6;
    elseif r < 0.93
        x_new = 0.2 * x - 0.26 * y;
        y_new = 0.23 * x + 0.22 * y + 1.6;
    else
        x_new = -0.15 * x + 0.28 * y;
        y_new = 0.26 * x + 0.24 * y + 0.44;
    end
    x = x_new;
    y = y_new;
end

x = 0;
y = 0;
iterations = 100000;

for i = 1:iterations
    [x, y] = barnsley_fern(x, y);
    plot(x, y, '.', 'MarkerSize', 1, 'Color', [0 0.5 0]);
end

set(gcf, 'Position', [100, 100, 800, 800]);
title('Barnsley Fern 分形');