h=1/(1-z); ha = PadeApproximant[h^a,{z,0,{3,3}}]