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 = a4
newb4p = (24*b4 - 2*a1*a2*g + 12*a3*g + a1^2*g^2)/24
newc4p = (24*c4 - 2*a2*b1*g + 12*b3*g + 2*a2*g^2 + a1*b1*g^2)/24
newa5p = (720*a5 - a1^4*g)/720
newb5p = (720*b5 + 60*a1*a3*g + 360*a4*g - a1^3*b1*g + 2*a1^3*g^2)/720
newc5p = (360*c5 - 60*a1*a3*g - a1^3*b1*g - 3*a1^3*g^2)/360
newd5p = (720*d5 - 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
newe5p = (360*e5 + 30*a2^2*g - 30*a3*b1*g - a1^2*b1^2*g + 30*a1*b3*g + 
     180*b4*g - 15*a1*a2*g^2 + 30*a3*g^2 - a1^2*b1*g^2 + 3*a1^2*g^3)/360
newf5p = (720*f5 - a1*b1^3*g - 60*b1*b3*g + 360*c4*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 = a6
newb6p = (1440*b6 + 2*a1^3*a2*g + 120*a1*a4*g + 720*a5*g - a1^4*g^2)/1440
newc6p = (6*c6 - a1*a4*g)/6
newd6p = (2160*d6 - 180*a2*a3*g - 2*a1^2*a2*b1*g - 60*a1*b4*g - 360*c5*g - 
     6*a1^2*a2*g^2 + 60*a1*a3*g^2 + a1^3*b1*g^2 + 3*a1^3*g^3)/2160
newe6p = (720*e6 + 60*a2*a3*g + 2*a1^2*a2*b1*g - 60*a1*b4*g + 360*c5*g + 
     6*a1^2*a2*g^2 - 60*a1*a3*g^2 - a1^3*b1*g^2 - 3*a1^3*g^3)/720
newf6p = (864*f6 + 2*a1^2*a2*b1*g - 72*a4*b1*g + 24*a1*b4*g + 432*b5*g + 
     144*c5*g + 12*a1*a3*g^2 + 72*a4*g^2 - a1^3*b1*g^2)/864
newg6p = (-4*a1*a2*b1^2*g - 120*a2*b3*g + 120*b1*b4*g - 240*a1*c4*g - 
     720*e5*g - 60*a2^2*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 - 120*b4*g^2 + 12*a1*a2*g^3 + 
     2*a1^2*b1*g^3 - a1^2*g^4 + 1440*g6)/1440
newh6p = (2*a1*a2*b1^2*g + 24*a2*b3*g - 48*b1*b4*g + 24*a1*c4*g + 144*d5*g + 
     288*e5*g + 12*a2^2*g^2 + 4*a1*a2*b1*g^2 - 24*a3*b1*g^2 - a1^2*b1^2*g^2 + 
     48*b4*g^2 - 4*a1*a2*g^3 - 2*a1^2*b1*g^3 + 288*h6)/288
newi6p = (2*a2*b1^3*g - 120*b1*c4*g + 720*f5*g + 8*a2*b1^2*g^2 - 
     a1*b1^3*g^2 - 60*b1*b3*g^2 + 120*c4*g^2 - 8*a2*b1*g^3 - 4*a1*b1^2*g^3 - 
     2*a2*g^4 - a1*b1*g^4 + 1440*i6)/1440
newa7p = (30240*a7 + a1^6*g)/30240
newb7p = (30240*b7 - 42*a1^3*a3*g + 2520*a1*a5*g + 15120*a6*g + a1^5*b1*g - 
     3*a1^5*g^2)/30240
newc7p = (30240*c7 - 5040*a1*a5*g - a1^5*b1*g + 3*a1^5*g^2)/30240
newd7p = (30240*d7 + 84*a1^3*a3*g + 5*a1^5*b1*g + 6*a1^5*g^2)/30240
newe7p = (30240*e7 + 2520*a3^2*g - 2520*a2*a4*g - 2520*a5*b1*g - 
     a1^4*b1^2*g + 42*a1^3*b3*g + 15120*b6*g - 2520*a1*c5*g + 
     21*a1^3*a2*g^2 + 420*a1^2*a3*g^2 + 1260*a1*a4*g^2 + 2520*a5*g^2 + 
     6*a1^4*b1*g^2 + 9*a1^4*g^3)/30240
newf7p = (30240*f7 - 42*a1^2*a2^2*g - 5040*a3^2*g + 7560*a2*a4*g - 
     42*a1^2*a3*b1*g + 5*a1^4*b1^2*g - 84*a1^3*b3*g + 2520*a1*b5*g + 
     5040*a1*c5*g + 15120*c6*g - 756*a1^2*a3*g^2 - 2520*a1*a4*g^2 - 
     9*a1^4*b1*g^2 - 24*a1^4*g^3)/30240
newg7p = (84*a1^2*a2^2*g + 5040*a3^2*g - 15120*a2*a4*g + 84*a1^2*a3*b1*g + 
     5*a1^4*b1^2*g + 168*a1^3*b3*g - 15120*a1*b5*g - 5040*a1*c5*g + 
     252*a1^2*a3*g^2 + 33*a1^4*b1*g^2 + 18*a1^4*g^3 + 90720*g7)/90720
newh7p = (-42*a1^2*a2^2*g + 2520*a3^2*g + 252*a1^2*a3*b1*g + 4*a1^4*b1^2*g - 
     84*a1^3*b3*g - 5040*a1*c5*g + 756*a1^2*a3*g^2 + 18*a1^4*b1*g^2 + 
     27*a1^4*g^3 + 30240*h7)/30240
newi7p = (168*a1^2*a2^2*g - 5040*a3^2*g - 336*a1^2*a3*b1*g + a1^4*b1^2*g + 
     336*a1^3*b3*g + 5040*a1*c5*g - 1008*a1^2*a3*g^2 - 6*a1^4*b1*g^2 - 
     30*a1^4*g^3 + 90720*i7)/90720
newj7p = (-84*a1*a3*b1^2*g - a1^3*b1^3*g - 5040*a3*b3*g + 5040*a2*b4*g + 
     5040*a1*e5*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*j7)/30240
newk7p = (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 - 10080*a2*b4*g - 5040*a1*d5*g - 
     10080*a1*e5*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*k7
     )/30240
newl7p = (168*a1*a2^2*b1*g - 294*a1*a3*b1^2*g + a1^3*b1^3*g - 10080*a3*b3*g + 
     252*a1^2*b1*b3*g + 2520*a2*b4*g + 2520*b1*c5*g + 7560*a1*d5*g + 
     45360*d6*g - 336*a1*a2^2*g^2 - 3780*a2*a3*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 - 1260*a1*b4*g^2 - 
     2520*c5*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*l7)/90720
newm7p = (-126*a1*a2^2*b1*g + 168*a1*a3*b1^2*g + 4*a1^3*b1^3*g - 
     2520*a3*b3*g - 168*a1^2*b1*b3*g + 5040*a2*b4*g - 2520*b1*c5*g + 
     2520*a1*e5*g + 15120*e6*g - 168*a1*a2^2*g^2 + 1260*a2*a3*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 - 1260*a1*b4*g^2 + 2520*c5*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*m7)/30240
newn7p = (-42*a1*a2^2*b1*g - 84*a1*a3*b1^2*g + 5*a1^3*b1^3*g + 2520*a3*b3*g + 
     126*a1^2*b1*b3*g - 2520*a2*b4*g - 7560*b1*b5*g - 2520*b1*c5*g + 
     45360*f6*g + 84*a1*a2^2*g^2 + 105*a1^2*a2*b1*g^2 - 84*a1*a3*b1*g^2 - 
     3780*a4*b1*g^2 + 18*a1^3*b1^2*g^2 + 168*a1^2*b3*g^2 + 1260*a1*b4*g^2 + 
     7560*b5*g^2 + 2520*c5*g^2 + 252*a1*a3*g^3 - 9*a1^3*b1*g^3 + 3*a1^3*g^4 + 
     90720*n7)/90720
newo7p = (42*a2^2*b1^2*g + a1^2*b1^4*g + 84*a1*b1^2*b3*g + 2520*b3^2*g - 
     5040*a2*c4*g - 5040*a1*f5*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*o7)/30240
newp7p = (42*a2^2*b1^2*g - 42*a3*b1^3*g - a1^2*b1^4*g - 5040*b3^2*g + 
     7560*a2*c4*g + 2520*b1*e5*g + 2520*a1*f5*g - 168*a2^2*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 + 1260*b1*b4*g^2 - 2520*a1*c4*g^2 - 2520*e5*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 + 
     15120*g*g6 + 30240*p7)/30240
newq7p = (-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*c4*g - 2520*b1*d5*g - 5040*b1*e5*g - 
     84*a2^2*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 - 2520*b1*b4*g^2 + 
     1260*a1*c4*g^2 + 2520*d5*g^2 + 5040*e5*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 + 15120*g*h6 + 30240*q7)/30240
newr7p = (a1*b1^5*g + 42*b1^3*b3*g - 2520*b1*f5*g + 21*a2*b1^3*g^2 + 
     6*a1*b1^4*g^2 + 168*b1^2*b3*g^2 - 1260*b1*c4*g^2 + 2520*f5*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 + 15120*g*i6 + 30240*r7)
    /30240
