By the way: It seems that nothing depends on T[0,0], so the call to limit in limit(Z0[j+1]/((j+1)!), T[0,0],0) seems unnecessary. Also, maybe you need to change one pi to %pi in the definition of C[0]. --Barton ________________________________