newa1p = a1
newb1p = b1 + 2*g
newa3p = (6*a3 + a1^2*g)/6
newb3p = (6*b3 + a1*b1*g + a1*g^2)/6
newa5p = (180*a5 - 60*a1*a3*g - a1^3*b1*g - 6*a1^3*g^2)/180
newb5p = (360*b5 - a1^2*b1^2*g - 120*a1*b3*g - 12*a1^2*b1*g^2 - 8*a1^2*g^3)/
   360
newc5p = (180*c5 - 30*a3*b1*g - a1^2*b1^2*g + 30*a1*b3*g - 30*a3*g^2 - 
     2*a1^2*b1*g^2 - 3*a1^2*g^3)/180
newd5p = (360*d5 - a1*b1^3*g - 60*b1*b3*g - 8*a1*b1^2*g^2 - 60*b3*g^2 - 
     14*a1*b1*g^3 - 7*a1*g^4)/360
newa7p = (3780*a7 + 630*a3^2*g - 1260*a1*a5*g + 63*a1^2*a3*b1*g + 
     a1^4*b1^2*g - 21*a1^3*b3*g + 378*a1^2*a3*g^2 + 9*a1^4*b1*g^2 + 
     27*a1^4*g^3)/3780
newb7p = (45360*b7 - 5040*a3^2*g + 5040*a1*a5*g - 336*a1^2*a3*b1*g + 
     a1^4*b1^2*g + 336*a1^3*b3*g - 2016*a1^2*a3*g^2 - 12*a1^4*b1*g^2 - 
     120*a1^4*g^3)/45360
newc7p = (15120*c7 - 84*a1*a3*b1^2*g - a1^3*b1^3*g - 5040*a3*b3*g + 
     5040*a1*c5*g - 1008*a1*a3*b1*g^2 - 24*a1^3*b1^2*g^2 - 672*a1*a3*g^3 - 
     88*a1^3*b1*g^3 - 72*a1^3*g^4)/15120
newd7p = (15120*d7 + 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 + 2016*a1*a3*b1*g^2 + 
     78*a1^3*b1^2*g^2 + 1008*a1^2*b3*g^2 + 1344*a1*a3*g^3 + 272*a1^3*b1*g^3 + 
     192*a1^3*g^4)/15120
newe7p = (45360*e7 + 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 + 2520*a5*g^2 - 
     1848*a1*a3*b1*g^2 - 18*a1^3*b1^2*g^2 - 1848*a1^2*b3*g^2 - 
     1512*a1*a3*g^3 - 234*a1^3*b1*g^3 - 180*a1^3*g^4)/45360
newf7p = (7560*f7 - 1260*a5*b1*g + 84*a1*a3*b1^2*g + 2*a1^3*b1^3*g - 
     1260*a3*b3*g - 84*a1^2*b1*b3*g + 1260*a1*c5*g - 1260*a5*g^2 + 
     168*a1*a3*b1*g^2 + 6*a1^3*b1^2*g^2 - 84*a1^2*b3*g^2 + 252*a1*a3*g^3 + 
     15*a1^3*b1*g^3 + 18*a1^3*g^4)/7560
newg7p = (a1^2*b1^4*g + 84*a1*b1^2*b3*g + 2520*b3^2*g - 5040*a1*d5*g + 
     18*a1^2*b1^3*g^2 + 1008*a1*b1*b3*g^2 + 92*a1^2*b1^2*g^3 + 
     672*a1*b3*g^3 + 120*a1^2*b1*g^4 + 48*a1^2*g^5 + 15120*g7)/15120
newh7p = (-42*a3*b1^3*g - a1^2*b1^4*g - 5040*b3^2*g + 2520*b1*c5*g + 
     2520*a1*d5*g - 336*a3*b1^2*g^2 - 18*a1^2*b1^3*g^2 - 840*a1*b1*b3*g^2 + 
     2520*c5*g^2 - 588*a3*b1*g^3 - 106*a1^2*b1^2*g^3 - 420*a1*b3*g^3 - 
     294*a3*g^4 - 148*a1^2*b1*g^4 - 69*a1^2*g^5 + 15120*h7)/15120
newi7p = (84*a3*b1^3*g + 5*a1^2*b1^4*g + 42*a1*b1^2*b3*g + 2520*b3^2*g - 
     2520*b1*b5*g - 5040*b1*c5*g + 672*a3*b1^2*g^2 + 48*a1^2*b1^3*g^2 + 
     504*a1*b1*b3*g^2 - 2520*b5*g^2 - 5040*c5*g^2 + 1176*a3*b1*g^3 + 
     173*a1^2*b1^2*g^3 + 336*a1*b3*g^3 + 588*a3*g^4 + 236*a1^2*b1*g^4 + 
     114*a1^2*g^5 + 15120*i7)/15120
newj7p = (a1*b1^5*g + 42*b1^3*b3*g - 2520*b1*d5*g + 12*a1*b1^4*g^2 + 
     336*b1^2*b3*g^2 - 2520*d5*g^2 + 53*a1*b1^3*g^3 + 588*b1*b3*g^3 + 
     104*a1*b1^2*g^4 + 294*b3*g^4 + 93*a1*b1*g^5 + 31*a1*g^6 + 15120*j7)/15120
