docs: def modular_exp( base, exponent, modulus ):
    "Raise base to exponent, reducing by modulus"
    if exponent < 0:
      raise NegativeExponentError( "Negative exponents (%d) not allowed" \
                                   % exponent )
    return pow( base, exponent, modulus )