newa1p = a1
newb1p = b1 + g
newa2p = (2*a2 - a1*g)/2
newa3p = (12*a3 + a1^2*g)/12
newb3p = (12*b3 + 6*a2*g + a1*b1*g - a1*g^2)/12
newa4p = (24*a4 - 2*a1*a2*g + 12*a3*g + a1^2*g^2)/24
newb4p = (24*b4 - 2*a2*b1*g + 12*b3*g + 2*a2*g^2 + a1*b1*g^2)/24
newa5p = (360*a5 - 60*a1*a3*g - a1^3*b1*g - 3*a1^3*g^2)/360
newb5p = (720*b5 - 60*a2^2*g - a1^2*b1^2*g - 120*a1*b3*g - 6*a1^2*b1*g^2 - 
     2*a1^2*g^3)/720
newc5p = (360*c5 + 30*a2^2*g + 180*a4*g - 30*a3*b1*g - a1^2*b1^2*g + 
     30*a1*b3*g - 15*a1*a2*g^2 + 30*a3*g^2 - a1^2*b1*g^2 + 3*a1^2*g^3)/360
newd5p = (720*d5 - a1*b1^3*g - 60*b1*b3*g + 360*b4*g - 30*a2*b1*g^2 - 
     4*a1*b1^2*g^2 + 60*b3*g^2 + 4*a1*b1*g^3 + a1*g^4)/720
newa6p = (2160*a6 - 180*a2*a3*g - 60*a1*a4*g - 360*a5*g - 2*a1^2*a2*b1*g - 
     6*a1^2*a2*g^2 + 60*a1*a3*g^2 + a1^3*b1*g^2 + 3*a1^3*g^3)/2160
newb6p = (720*b6 + 60*a2*a3*g - 60*a1*a4*g + 360*a5*g + 2*a1^2*a2*b1*g + 
     6*a1^2*a2*g^2 - 60*a1*a3*g^2 - a1^3*b1*g^2 - 3*a1^3*g^3)/720
newc6p = (1440*c6 + 120*a4*b1*g - 4*a1*a2*b1^2*g - 120*a2*b3*g - 
     240*a1*b4*g - 720*c5*g - 60*a2^2*g^2 - 120*a4*g^2 - 4*a1*a2*b1*g^2 + 
     60*a3*b1*g^2 + 2*a1^2*b1^2*g^2 - 60*a1*b3*g^2 + 12*a1*a2*g^3 + 
     2*a1^2*b1*g^3 - a1^2*g^4)/1440
newd6p = (288*d6 - 48*a4*b1*g + 2*a1*a2*b1^2*g + 24*a2*b3*g + 24*a1*b4*g + 
     144*b5*g + 288*c5*g + 12*a2^2*g^2 + 48*a4*g^2 + 4*a1*a2*b1*g^2 - 
     24*a3*b1*g^2 - a1^2*b1^2*g^2 - 4*a1*a2*g^3 - 2*a1^2*b1*g^3)/288
newe6p = (1440*e6 + 2*a2*b1^3*g - 120*b1*b4*g + 720*d5*g + 8*a2*b1^2*g^2 - 
     a1*b1^3*g^2 - 60*b1*b3*g^2 + 120*b4*g^2 - 8*a2*b1*g^3 - 4*a1*b1^2*g^3 - 
     2*a2*g^4 - a1*b1*g^4)/1440
newa7p = (30240*a7 - 42*a1^2*a2^2*g + 2520*a3^2*g - 5040*a1*a5*g + 
     252*a1^2*a3*b1*g + 4*a1^4*b1^2*g - 84*a1^3*b3*g + 756*a1^2*a3*g^2 + 
     18*a1^4*b1*g^2 + 27*a1^4*g^3)/30240
newb7p = (90720*b7 + 168*a1^2*a2^2*g - 5040*a3^2*g + 5040*a1*a5*g - 
     336*a1^2*a3*b1*g + a1^4*b1^2*g + 336*a1^3*b3*g - 1008*a1^2*a3*g^2 - 
     6*a1^4*b1*g^2 - 30*a1^4*g^3)/90720
newc7p = (30240*c7 + 5040*a2*a4*g - 84*a1*a3*b1^2*g - a1^3*b1^3*g - 
     5040*a3*b3*g + 5040*a1*c5*g - 504*a1*a3*b1*g^2 - 12*a1^3*b1^2*g^2 - 
     168*a1*a3*g^3 - 22*a1^3*b1*g^3 - 9*a1^3*g^4)/30240
newd7p = (30240*d7 - 10080*a2*a4*g + 84*a1*a2^2*b1*g + 168*a1*a3*b1^2*g + 
     5*a1^3*b1^3*g + 10080*a3*b3*g + 168*a1^2*b1*b3*g - 5040*a1*b5*g - 
     10080*a1*c5*g + 252*a1*a2^2*g^2 + 1008*a1*a3*b1*g^2 + 39*a1^3*b1^2*g^2 + 
     504*a1^2*b3*g^2 + 336*a1*a3*g^3 + 68*a1^3*b1*g^3 + 24*a1^3*g^4)/30240
newe7p = (90720*e7 + 2520*a2*a4*g + 45360*a6*g + 168*a1*a2^2*b1*g + 
     2520*a5*b1*g - 294*a1*a3*b1^2*g + a1^3*b1^3*g - 10080*a3*b3*g + 
     252*a1^2*b1*b3*g + 7560*a1*b5*g - 336*a1*a2^2*g^2 - 3780*a2*a3*g^2 - 
     1260*a1*a4*g^2 - 2520*a5*g^2 - 42*a1^2*a2*b1*g^2 - 924*a1*a3*b1*g^2 - 
     9*a1^3*b1^2*g^2 - 924*a1^2*b3*g^2 - 126*a1^2*a2*g^3 + 252*a1*a3*g^3 - 
     48*a1^3*b1*g^3 + 9*a1^3*g^4)/90720
newf7p = (30240*f7 + 5040*a2*a4*g - 126*a1*a2^2*b1*g - 2520*a5*b1*g + 
     168*a1*a3*b1^2*g + 4*a1^3*b1^3*g - 2520*a3*b3*g - 168*a1^2*b1*b3*g + 
     15120*b6*g + 2520*a1*c5*g - 168*a1*a2^2*g^2 + 1260*a2*a3*g^2 - 
     1260*a1*a4*g^2 + 2520*a5*g^2 + 42*a1^2*a2*b1*g^2 + 168*a1*a3*b1*g^2 + 
     6*a1^3*b1^2*g^2 - 84*a1^2*b3*g^2 + 126*a1^2*a2*g^3 - 504*a1*a3*g^3 - 
     3*a1^3*b1*g^3 - 27*a1^3*g^4)/30240
newg7p = (42*a2^2*b1^2*g + a1^2*b1^4*g + 84*a1*b1^2*b3*g + 2520*b3^2*g - 
     5040*a2*b4*g - 5040*a1*d5*g + 252*a2^2*b1*g^2 + 9*a1^2*b1^3*g^2 + 
     504*a1*b1*b3*g^2 + 84*a2^2*g^3 + 23*a1^2*b1^2*g^3 + 168*a1*b3*g^3 + 
     15*a1^2*b1*g^4 + 3*a1^2*g^5 + 30240*g7)/30240
newh7p = (42*a2^2*b1^2*g - 42*a3*b1^3*g - a1^2*b1^4*g - 5040*b3^2*g + 
     7560*a2*b4*g + 2520*b1*c5*g + 15120*c6*g + 2520*a1*d5*g - 
     168*a2^2*b1*g^2 + 1260*a4*b1*g^2 - 42*a1*a2*b1^2*g^2 - 168*a3*b1^2*g^2 - 
     9*a1^2*b1^3*g^2 - 1260*a2*b3*g^2 - 420*a1*b1*b3*g^2 - 2520*a1*b4*g^2 - 
     2520*c5*g^2 - 336*a2^2*g^3 - 42*a1*a2*b1*g^3 + 168*a3*b1*g^3 - 
     16*a1^2*b1^2*g^3 - 420*a1*b3*g^3 + 21*a1*a2*g^4 + 42*a3*g^4 - 
     8*a1^2*b1*g^4 - 3*a1^2*g^5 + 30240*h7)/30240
newi7p = (-84*a2^2*b1^2*g + 84*a3*b1^3*g + 5*a1^2*b1^4*g + 42*a1*b1^2*b3*g + 
     2520*b3^2*g - 2520*a2*b4*g - 2520*b1*b5*g - 5040*b1*c5*g + 15120*d6*g - 
     84*a2^2*b1*g^2 - 2520*a4*b1*g^2 + 105*a1*a2*b1^2*g^2 + 336*a3*b1^2*g^2 + 
     24*a1^2*b1^3*g^2 + 1260*a2*b3*g^2 + 252*a1*b1*b3*g^2 + 1260*a1*b4*g^2 + 
     2520*b5*g^2 + 5040*c5*g^2 + 252*a2^2*g^3 + 210*a1*a2*b1*g^3 - 
     336*a3*b1*g^3 + 17*a1^2*b1^2*g^3 + 84*a1*b3*g^3 - 84*a3*g^4 - 
     23*a1^2*b1*g^4 - 6*a1^2*g^5 + 30240*i7)/30240
newj7p = (a1*b1^5*g + 42*b1^3*b3*g - 2520*b1*d5*g + 15120*e6*g + 
     21*a2*b1^3*g^2 + 6*a1*b1^4*g^2 + 168*b1^2*b3*g^2 - 1260*b1*b4*g^2 + 
     2520*d5*g^2 + 84*a2*b1^2*g^3 + 8*a1*b1^3*g^3 - 168*b1*b3*g^3 + 
     21*a2*b1*g^4 - 8*a1*b1^2*g^4 - 42*b3*g^4 - 6*a1*b1*g^5 - a1*g^6 + 
     30240*j7)/30240
