function kapradina

n = 100000;
X = zeros(n,1);
Y = zeros(n,1);

for i = 2:n
    r = rand;
    
    if(r < 0.01)
        X(i) = 0;
        Y(i) = 0.16*Y(i-1);
    elseif(r > 0.01 && r < 0.87)
        X(i) = 0.85*X(i-1) + 0.04*Y(i-1);
        Y(i) = -0.04*X(i-1)+ 0.85*Y(i-1) + 1.6;
    elseif(r > 0.87 && r < 0.94)
        X(i) = 0.2*X(i-1) -0.26*Y(i-1);
        Y(i) = 0.23*X(i-1)+ 0.22*Y(i-1) + 1.6;
    else
        X(i) = -0.15*X(i-1) + 0.28*Y(i-1);
        Y(i) = 0.26*X(i-1) + 0.24*Y(i-1) + 0.44;
    end
end
figure
plot(X,Y,'.k');
axis equal;