Cornish-Fisher expansion of w

We apply cfexpw to cumulantw.

The following zformula is defined as .

-v0 + w + o*(-cr[0] - td[al1, au1] + 
   w^2*(-cr[2] + tp3[9, 9, 9]/6) - tp3[9, 9, 9]/6 - 
   (v0^2*tp3[9, 9, 9])/3 + (v0*w*tp3[9, 9, 9])/6) + 
 o^2*(v0^3*(tp3[9, 9, 9]^2/18 + (tp3[9, 9, al1]*tp3[9, 9, au1])/8 - 
     tp4[9, 9, 9, 9]/8) + v0*(td[al1, au2]*td[al2, au1] - 
     (cr[0]*tp3[9, 9, 9])/6 - (td[al1, au1]*tp3[9, 9, 9])/6 + 
     (5*tp3[9, 9, 9]^2)/72 + (tp3[9, 9, al1]*tp3[9, 9, au1])/8 - 
     tp4[9, 9, 9, 9]/24) + v0*w^2*(-(cr[2]*tp3[9, 9, 9])/6 - 
     tp3[9, 9, 9]^2/24 + tp4[9, 9, 9, 9]/24) + 
   w^3*(-cr[3] - (cr[2]*tp3[9, 9, 9])/3 - tp3[9, 9, 9]^2/72 + 
     tp4[9, 9, 9, 9]/24) + w*(-cr[1] + td[al1, au2]*td[al2, au1] - 
     (cr[0]*tp3[9, 9, 9])/3 + (td[al1, au1]*tp3[9, 9, 9])/6 + 
     (13*tp3[9, 9, 9]^2)/72 + (tp3[9, 9, al1]*tp3[9, 9, au1])/2 - 
     td[au1, au2]*tp3[9, al1, al2] + 
     (tp3[9, al1, au2]*tp3[9, al2, au1])/2 + 
     v0^2*(-(tp3[9, 9, al1]*tp3[9, 9, au1])/8 + 
       tp4[9, 9, 9, 9]/24) - tp4[9, 9, 9, 9]/8 - 
     tp4[9, 9, al1, au1]/4))

Get the coefficients of for zformula.