newa1 = a1
newb1 = b1 + g
newa2 = a2 - (a1*g)/2
newa3 = a3 + (a1^2*g)/12
newb3 = b3 + (a2*g)/2 + (a1*b1*g)/12 - (a1*g^2)/12
newa4 = a4 - (a1*a2*g)/12 + (a3*g)/2 + (a1^2*g^2)/24
newb4 = b4 - (a2*b1*g)/12 + (b3*g)/2 + (a2*g^2)/12 + (a1*b1*g^2)/24
newa5 = a5 - (a1*a3*g)/6 - (a1^3*b1*g)/360 - (a1^3*g^2)/120
newb5 = b5 - (a2^2*g)/12 - (a1^2*b1^2*g)/720 - (a1*b3*g)/6 - 
   (a1^2*b1*g^2)/120 - (a1^2*g^3)/360
newc5 = c5 + (a2^2*g)/12 + (a4*g)/2 - (a3*b1*g)/12 - (a1^2*b1^2*g)/360 + 
   (a1*b3*g)/12 - (a1*a2*g^2)/24 + (a3*g^2)/12 - (a1^2*b1*g^2)/360 + 
   (a1^2*g^3)/120
newd5 = d5 - (a1*b1^3*g)/720 - (b1*b3*g)/12 + (b4*g)/2 - (a2*b1*g^2)/24 - 
   (a1*b1^2*g^2)/180 + (b3*g^2)/12 + (a1*b1*g^3)/180 + (a1*g^4)/720
newa6 = a6 - (a2*a3*g)/12 - (a1*a4*g)/36 - (a5*g)/6 - (a1^2*a2*b1*g)/1080 - 
   (a1^2*a2*g^2)/360 + (a1*a3*g^2)/36 + (a1^3*b1*g^2)/2160 + (a1^3*g^3)/720
newb6 = b6 + (a2*a3*g)/12 - (a1*a4*g)/12 + (a5*g)/2 + (a1^2*a2*b1*g)/360 + 
   (a1^2*a2*g^2)/120 - (a1*a3*g^2)/12 - (a1^3*b1*g^2)/720 - (a1^3*g^3)/240
newc6 = c6 + (a4*b1*g)/12 - (a1*a2*b1^2*g)/360 - (a2*b3*g)/12 - (a1*b4*g)/6 - 
   (c5*g)/2 - (a2^2*g^2)/24 - (a4*g^2)/12 - (a1*a2*b1*g^2)/360 + 
   (a3*b1*g^2)/24 + (a1^2*b1^2*g^2)/720 - (a1*b3*g^2)/24 + (a1*a2*g^3)/120 + 
   (a1^2*b1*g^3)/720 - (a1^2*g^4)/1440
newd6 = d6 - (a4*b1*g)/6 + (a1*a2*b1^2*g)/144 + (a2*b3*g)/12 + (a1*b4*g)/12 + 
   (b5*g)/2 + c5*g + (a2^2*g^2)/24 + (a4*g^2)/6 + (a1*a2*b1*g^2)/72 - 
   (a3*b1*g^2)/12 - (a1^2*b1^2*g^2)/288 - (a1*a2*g^3)/72 - (a1^2*b1*g^3)/144
newe6 = e6 + (a2*b1^3*g)/720 - (b1*b4*g)/12 + (d5*g)/2 + (a2*b1^2*g^2)/180 - 
   (a1*b1^3*g^2)/1440 - (b1*b3*g^2)/24 + (b4*g^2)/12 - (a2*b1*g^3)/180 - 
   (a1*b1^2*g^3)/360 - (a2*g^4)/720 - (a1*b1*g^4)/1440
newa7 = a7 - (a1^2*a2^2*g)/720 + (a3^2*g)/12 - (a1*a5*g)/6 + 
   (a1^2*a3*b1*g)/120 + (a1^4*b1^2*g)/7560 - (a1^3*b3*g)/360 + 
   (a1^2*a3*g^2)/40 + (a1^4*b1*g^2)/1680 + (a1^4*g^3)/1120
newb7 = b7 + (a1^2*a2^2*g)/540 - (a3^2*g)/18 + (a1*a5*g)/18 - 
   (a1^2*a3*b1*g)/270 + (a1^4*b1^2*g)/90720 + (a1^3*b3*g)/270 - 
   (a1^2*a3*g^2)/90 - (a1^4*b1*g^2)/15120 - (a1^4*g^3)/3024
newc7 = c7 + (a2*a4*g)/6 - (a1*a3*b1^2*g)/360 - (a1^3*b1^3*g)/30240 - 
   (a3*b3*g)/6 + (a1*c5*g)/6 - (a1*a3*b1*g^2)/60 - (a1^3*b1^2*g^2)/2520 - 
   (a1*a3*g^3)/180 - (11*a1^3*b1*g^3)/15120 - (a1^3*g^4)/3360
newd7 = d7 - (a2*a4*g)/3 + (a1*a2^2*b1*g)/360 + (a1*a3*b1^2*g)/180 + 
   (a1^3*b1^3*g)/6048 + (a3*b3*g)/3 + (a1^2*b1*b3*g)/180 - (a1*b5*g)/6 - 
   (a1*c5*g)/3 + (a1*a2^2*g^2)/120 + (a1*a3*b1*g^2)/30 + 
   (13*a1^3*b1^2*g^2)/10080 + (a1^2*b3*g^2)/60 + (a1*a3*g^3)/90 + 
   (17*a1^3*b1*g^3)/7560 + (a1^3*g^4)/1260
newe7 = e7 + (a2*a4*g)/36 + (a6*g)/2 + (a1*a2^2*b1*g)/540 + (a5*b1*g)/36 - 
   (7*a1*a3*b1^2*g)/2160 + (a1^3*b1^3*g)/90720 - (a3*b3*g)/9 + 
   (a1^2*b1*b3*g)/360 + (a1*b5*g)/12 - (a1*a2^2*g^2)/270 - (a2*a3*g^2)/24 - 
   (a1*a4*g^2)/72 - (a5*g^2)/36 - (a1^2*a2*b1*g^2)/2160 - 
   (11*a1*a3*b1*g^2)/1080 - (a1^3*b1^2*g^2)/10080 - (11*a1^2*b3*g^2)/1080 - 
   (a1^2*a2*g^3)/720 + (a1*a3*g^3)/360 - (a1^3*b1*g^3)/1890 + (a1^3*g^4)/10080
newf7 = f7 + (a2*a4*g)/6 - (a1*a2^2*b1*g)/240 - (a5*b1*g)/12 + 
   (a1*a3*b1^2*g)/180 + (a1^3*b1^3*g)/7560 - (a3*b3*g)/12 - 
   (a1^2*b1*b3*g)/180 + (b6*g)/2 + (a1*c5*g)/12 - (a1*a2^2*g^2)/180 + 
   (a2*a3*g^2)/24 - (a1*a4*g^2)/24 + (a5*g^2)/12 + (a1^2*a2*b1*g^2)/720 + 
   (a1*a3*b1*g^2)/180 + (a1^3*b1^2*g^2)/5040 - (a1^2*b3*g^2)/360 + 
   (a1^2*a2*g^3)/240 - (a1*a3*g^3)/60 - (a1^3*b1*g^3)/10080 - (a1^3*g^4)/1120
newg7 = (a2^2*b1^2*g)/720 + (a1^2*b1^4*g)/30240 + (a1*b1^2*b3*g)/360 + 
   (b3^2*g)/12 - (a2*b4*g)/6 - (a1*d5*g)/6 + (a2^2*b1*g^2)/120 + 
   (a1^2*b1^3*g^2)/3360 + (a1*b1*b3*g^2)/60 + (a2^2*g^3)/360 + 
   (23*a1^2*b1^2*g^3)/30240 + (a1*b3*g^3)/180 + (a1^2*b1*g^4)/2016 + 
   (a1^2*g^5)/10080 + g7
newh7 = (a2^2*b1^2*g)/720 - (a3*b1^3*g)/720 - (a1^2*b1^4*g)/30240 - 
   (b3^2*g)/6 + (a2*b4*g)/4 + (b1*c5*g)/12 + (c6*g)/2 + (a1*d5*g)/12 - 
   (a2^2*b1*g^2)/180 + (a4*b1*g^2)/24 - (a1*a2*b1^2*g^2)/720 - 
   (a3*b1^2*g^2)/180 - (a1^2*b1^3*g^2)/3360 - (a2*b3*g^2)/24 - 
   (a1*b1*b3*g^2)/72 - (a1*b4*g^2)/12 - (c5*g^2)/12 - (a2^2*g^3)/90 - 
   (a1*a2*b1*g^3)/720 + (a3*b1*g^3)/180 - (a1^2*b1^2*g^3)/1890 - 
   (a1*b3*g^3)/72 + (a1*a2*g^4)/1440 + (a3*g^4)/720 - (a1^2*b1*g^4)/3780 - 
   (a1^2*g^5)/10080 + h7
newi7 = -(a2^2*b1^2*g)/360 + (a3*b1^3*g)/360 + (a1^2*b1^4*g)/6048 + 
   (a1*b1^2*b3*g)/720 + (b3^2*g)/12 - (a2*b4*g)/12 - (b1*b5*g)/12 - 
   (b1*c5*g)/6 + (d6*g)/2 - (a2^2*b1*g^2)/360 - (a4*b1*g^2)/12 + 
   (a1*a2*b1^2*g^2)/288 + (a3*b1^2*g^2)/90 + (a1^2*b1^3*g^2)/1260 + 
   (a2*b3*g^2)/24 + (a1*b1*b3*g^2)/120 + (a1*b4*g^2)/24 + (b5*g^2)/12 + 
   (c5*g^2)/6 + (a2^2*g^3)/120 + (a1*a2*b1*g^3)/144 - (a3*b1*g^3)/90 + 
   (17*a1^2*b1^2*g^3)/30240 + (a1*b3*g^3)/360 - (a3*g^4)/360 - 
   (23*a1^2*b1*g^4)/30240 - (a1^2*g^5)/5040 + i7
newj7 = (a1*b1^5*g)/30240 + (b1^3*b3*g)/720 - (b1*d5*g)/12 + (e6*g)/2 + 
   (a2*b1^3*g^2)/1440 + (a1*b1^4*g^2)/5040 + (b1^2*b3*g^2)/180 - 
   (b1*b4*g^2)/24 + (d5*g^2)/12 + (a2*b1^2*g^3)/360 + (a1*b1^3*g^3)/3780 - 
   (b1*b3*g^3)/180 + (a2*b1*g^4)/1440 - (a1*b1^2*g^4)/3780 - (b3*g^4)/720 - 
   (a1*b1*g^5)/5040 - (a1*g^6)/30240 + j7
