function Kohova_pahulja(Iter)

if nargin ~= 1
	error('Nema dovoljno ulaznih podataka');
end

xl = zeros(10,1);
xr = xl;
yl = xl;
yr = yl;
xr(Iter) = 1;
r = sqrt(1/3^2 - 1/6^2);
clf;
% podesavanje prikaza ...
set(gca,'FontSize',14);
set(gcf,'Color',[1,1,1]);
hold on;
pahuljica(xl,xr,yl,yr,Iter,r);
title('Kohova pahulja');
text(0.5,-0.05,(['Broj iteracija: ' num2str(Iter)]), ...
		 'HorizontalAlign','center','FontSize',12); 
hold off;
axis equal; axis tight; axis off;

%------------------------------------------------------------
function pahuljica(xl,xr,yl,yr,Iter,r)

if (Iter < 2)
    plot([xl(1) xr(1)],[-yl(1) -yr(1)],'b-') 
    return
end

Iter = Iter-1;

xl(Iter) = xl(Iter+1);
yl(Iter) = yl(Iter+1);
xr(Iter) = 1/3 * xr(Iter+1) + 2/3 * xl(Iter+1);
yr(Iter) = 1/3 * yr(Iter+1) + 2/3 * yl(Iter+1);
pahuljica(xl,xr,yl,yr,Iter,r);

xl(Iter) = xr(Iter);
yl(Iter) = yr(Iter);
xr(Iter) = .5 * xr(Iter+1) + .5 * xl(Iter+1) - r * (yl(Iter+1) - yr(Iter + 1));
yr(Iter) = .5 * yr(Iter+1)+.5*yl(Iter+1)+r*(xl(Iter+1)-xr(Iter+1));
pahuljica(xl,xr,yl,yr,Iter,r);

xl(Iter) = xr(Iter);
yl(Iter) = yr(Iter);
xr(Iter) = 2/3 * xr(Iter+1) + 1/3 * xl(Iter+1);
yr(Iter) = 2/3 * yr(Iter+1) + 1/3 * yl(Iter+1);
pahuljica(xl,xr,yl,yr,Iter,r);

xl(Iter) = xr(Iter);
yl(Iter) = yr(Iter);
xr(Iter) = xr(Iter+1);
yr(Iter) = yr(Iter+1);
pahuljica(xl,xr,yl,yr,Iter,r);
Iter = Iter + 1;
return;