exponentmod.m

function a=exponentmod(b,e,m)
% This fucntion calculates 
% [ b^e (mod m)]
if m==1
    a=0;
    return
end
i=1;
a=1;
while i<=e
    a = mod(a*b,m);
    i = i+1;
end
end