fractal2.m

10000 puntos



1000 puntos



function X=fractal(N)

close allX=zeros(N,2);
X(1,:)=[0,0];
for K=1:N-1
r=rand;
if r<0.01
X(K+1,:)=T(X(K,:), 0, 0, 0, 0.16, 0, 0);
elseif r<0.86
X(K+1,:)=T(X(K,:), 0.85, 0.04, -0.04, 0.85, 0, 1.6);
elseif r<0.93
X(K+1,:)=T(X(K,:), 0.2, -0.26, 0.23, 0.22, 0, 1.6);
else
X(K+1,:)=T(X(K,:), -.15, 0.28, 0.26, 0.24, 0, 0.44);
end
end
figure(1)
plot(X(:,1), X(:,2),'.');

function U=T(X,a,b,c,d,e,f)
U=zeros(1,2);
U(1,1)=a*X(1)+b*X(2)+e;
U(1,2)=c*X(1)+d*X(2)+f;