Notes

%n = size of matrix
%l = number of fractal levels
function Final_Lab(n,l)
mold=zeros(n);
mnew=zeros(n);
mold(:,1)=1;
mold(1,:)=1;
mold(n,:)=1;
mold(:,n)=1;
%mold
a = 0.5;
b = 0;
c = 0;
d = 0.5;
e = [1, 1, n/2];
f = [1, n/2, n/2];
%loop over each row and coloum.
for z=1:l
for x=1:n
for y=1:n
if mold(x,y)==1 %determine the cooridinate = 1
x1=floor(a*x+b*y+e(1));
x2=floor(a*x+b*y+e(2));
x3=floor(a*x+b*y+e(3));
y1=floor(c*x+d*y+f(1));
y2=floor(c*x+d*y+f(2));
y3=floor(c*x+d*y+f(3));
mnew(x1,y1)=1;
mnew(x2,y2)=1;
mnew(x3,y3)=1;
end
end
end
mold=mnew; %save the coordinates into mnew
mnew=zeros(n);
end
figure
hold on
for x=1:n
for y=1:n-1
if mold(x,y)==1 && mold(x,y+1)==1 %see if verticle adjacents are both equal to 1
plot([x,x],[y,y+1], 'b')%plot those lines
end
end
end
for x=1:n-1
for y=1:n
if mold(x,y)== 1 && mold(x+1,y)==1 %see if horizontal adjacents are 1
plot([x,x+1],[y,y], 'b') %plot those lines
end
end
end


end