function c = mocnina_modulo(a,b,n)
%function c = mocnina_modulo(a,b,n)
%   mocnina a^b modulo n

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

end

