create fractals in matlab

main.m
% Define the initial triangle
p1 = [0,0];
p2 = [1,0];
theta = 60 * pi / 180;
p3 = [0.5*sin(theta), 0.5*cos(theta)];

% Generate the Koch snowflake
n = 4; % Number of iterations
figure;
hold on;
koch_snowflake(p1, p2, n);
koch_snowflake(p2, p3, n);
koch_snowflake(p3, p1, n);
axis equal;
title(['Koch Snowflake, n=', num2str(n)]);

function koch_snowflake(p1, p2, n)
% Recursive function to generate a Koch snowflake
if n == 0
    plot([p1(1), p2(1)], [p1(2), p2(2)], 'k');
else
    delta = (p2 - p1) / 3;
    r = [cos(pi/3), -sin(pi/3); sin(pi/3), cos(pi/3)];
    q1 = p1 + delta;
    q2 = p2 - delta;
    q3 = q1 + r*(q2 - q1);
    koch_snowflake(p1, q1, n-1);
    koch_snowflake(q1, q3, n-1);
    koch_snowflake(q3, q2, n-1);
    koch_snowflake(q2, p2, n-1);
end
end
769 chars
33 lines