FUNCTION :  ETA[POWERqMODP] -  q-expansion of power of q-series mod p

CALLING SEQUENCE :  POWERqMODP(a,N,M,p)

PARAMETERS :
    a - q-series    
    M,p - positive integers
    N - integer

GLOBAL VARIABLES : 

SYNOPSIS :   
   Compute a^N to O(q^M) mod p (where N is ANY integer)

EXAMPLES :   

DISCUSSION :

SEE ALSO :