function s=sierp_koberec(n)
%sierp_koberec - Siešrpenskeho
%n...pocet iteraci (zaroven urcuje rad matice: 3^n), n<8
%s..koberec (matice 3^n x 3^n)
%priklad: imagesc(sierp_koberec(5));

if nargin<1
    n=3;
end

if ~isscalar(n) || ~isreal(n) || n<=0 || n>8 || n~=fix(n)
   error('Zadejte kladne cele cislo mezi 1 a 8!') 
end

s=1;%priprava promenne kterou budeme nafukovat
for i=1:n
    m=size(s,1); %pocet radku matice 's'
    s=[s s s;s zeros(m) s; s s s];    %skladame z malych kousku
end