Computes the b-th power of a in ℤ/pℤ. This function is here to also deal with large p.
i1 : modularPower(8,10003,101) o1 = 7