function c = mocnina(a,b)
%function c = mocnina(a,b)
%   mocnina a^b

if b<0, error('Zaporny exponent'); end
if round(b)~=b, error('Necelociselny exponent'); end
c=1;
while b>0
    if mod(b,2)
        c=c*a;
    end
    b=floor(b/2);
    a=a*a;
end

end

