function koch
n = 5;
str = 'F';
alfa = 60;
for i = 1:n
    ns = length(str);
    spom = [];
    for j = 1:ns
        if(str(j) == 'F')
            spom = [spom,'F+F--F+F'];
        else
            spom = [spom,str(j)];
        end
    end
    str = spom;
end

% zelvi grafika
x = 0;
y = 0;
beta = 0;
for i = 1:length(str)
    switch(str(i))
        case 'F'
            x = [x,x(end)+cos(beta)];
            y = [y,y(end)+sin(beta)];
        case '+'
            beta = beta + alfa*pi/180;
        case '-'
            beta = beta - alfa*pi/180;
    end
end
figure
plot(x,y);
axis equal;