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
newb4 = b4 - (a1*a2*g)/12 + (a3*g)/2 + (a1^2*g^2)/24
newc4 = c4 - (a2*b1*g)/12 + (b3*g)/2 + (a2*g^2)/12 + (a1*b1*g^2)/24
newa5 = a5 - (a1^4*g)/720
newb5 = b5 + (a1*a3*g)/12 + (a4*g)/2 - (a1^3*b1*g)/720 + (a1^3*g^2)/360
newc5 = c5 - (a1*a3*g)/6 - (a1^3*b1*g)/360 - (a1^3*g^2)/120
newd5 = d5 - (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
newe5 = e5 + (a2^2*g)/12 - (a3*b1*g)/12 - (a1^2*b1^2*g)/360 + (a1*b3*g)/12 + 
   (b4*g)/2 - (a1*a2*g^2)/24 + (a3*g^2)/12 - (a1^2*b1*g^2)/360 + 
   (a1^2*g^3)/120
newf5 = f5 - (a1*b1^3*g)/720 - (b1*b3*g)/12 + (c4*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
newb6 = b6 + (a1^3*a2*g)/720 + (a1*a4*g)/12 + (a5*g)/2 - (a1^4*g^2)/1440
newc6 = c6 - (a1*a4*g)/6
newd6 = d6 - (a2*a3*g)/12 - (a1^2*a2*b1*g)/1080 - (a1*b4*g)/36 - (c5*g)/6 - 
   (a1^2*a2*g^2)/360 + (a1*a3*g^2)/36 + (a1^3*b1*g^2)/2160 + (a1^3*g^3)/720
newe6 = e6 + (a2*a3*g)/12 + (a1^2*a2*b1*g)/360 - (a1*b4*g)/12 + (c5*g)/2 + 
   (a1^2*a2*g^2)/120 - (a1*a3*g^2)/12 - (a1^3*b1*g^2)/720 - (a1^3*g^3)/240
newf6 = f6 + (a1^2*a2*b1*g)/432 - (a4*b1*g)/12 + (a1*b4*g)/36 + (b5*g)/2 + 
   (c5*g)/6 + (a1*a3*g^2)/72 + (a4*g^2)/12 - (a1^3*b1*g^2)/864
newg6 = -(a1*a2*b1^2*g)/360 - (a2*b3*g)/12 + (b1*b4*g)/12 - (a1*c4*g)/6 - 
   (e5*g)/2 - (a2^2*g^2)/24 - (a1*a2*b1*g^2)/360 + (a3*b1*g^2)/24 + 
   (a1^2*b1^2*g^2)/720 - (a1*b3*g^2)/24 - (b4*g^2)/12 + (a1*a2*g^3)/120 + 
   (a1^2*b1*g^3)/720 - (a1^2*g^4)/1440 + g6
newh6 = (a1*a2*b1^2*g)/144 + (a2*b3*g)/12 - (b1*b4*g)/6 + (a1*c4*g)/12 + 
   (d5*g)/2 + e5*g + (a2^2*g^2)/24 + (a1*a2*b1*g^2)/72 - (a3*b1*g^2)/12 - 
   (a1^2*b1^2*g^2)/288 + (b4*g^2)/6 - (a1*a2*g^3)/72 - (a1^2*b1*g^3)/144 + h6
newi6 = (a2*b1^3*g)/720 - (b1*c4*g)/12 + (f5*g)/2 + (a2*b1^2*g^2)/180 - 
   (a1*b1^3*g^2)/1440 - (b1*b3*g^2)/24 + (c4*g^2)/12 - (a2*b1*g^3)/180 - 
   (a1*b1^2*g^3)/360 - (a2*g^4)/720 - (a1*b1*g^4)/1440 + i6
newa7 = a7 + (a1^6*g)/30240
newb7 = b7 - (a1^3*a3*g)/720 + (a1*a5*g)/12 + (a6*g)/2 + (a1^5*b1*g)/30240 - 
   (a1^5*g^2)/10080
newc7 = c7 - (a1*a5*g)/6 - (a1^5*b1*g)/30240 + (a1^5*g^2)/10080
newd7 = d7 + (a1^3*a3*g)/360 + (a1^5*b1*g)/6048 + (a1^5*g^2)/5040
newe7 = e7 + (a3^2*g)/12 - (a2*a4*g)/12 - (a5*b1*g)/12 - 
   (a1^4*b1^2*g)/30240 + (a1^3*b3*g)/720 + (b6*g)/2 - (a1*c5*g)/12 + 
   (a1^3*a2*g^2)/1440 + (a1^2*a3*g^2)/72 + (a1*a4*g^2)/24 + (a5*g^2)/12 + 
   (a1^4*b1*g^2)/5040 + (a1^4*g^3)/3360
newf7 = f7 - (a1^2*a2^2*g)/720 - (a3^2*g)/6 + (a2*a4*g)/4 - 
   (a1^2*a3*b1*g)/720 + (a1^4*b1^2*g)/6048 - (a1^3*b3*g)/360 + (a1*b5*g)/12 + 
   (a1*c5*g)/6 + (c6*g)/2 - (a1^2*a3*g^2)/40 - (a1*a4*g^2)/12 - 
   (a1^4*b1*g^2)/3360 - (a1^4*g^3)/1260
newg7 = (a1^2*a2^2*g)/1080 + (a3^2*g)/18 - (a2*a4*g)/6 + 
   (a1^2*a3*b1*g)/1080 + (a1^4*b1^2*g)/18144 + (a1^3*b3*g)/540 - 
   (a1*b5*g)/6 - (a1*c5*g)/18 + (a1^2*a3*g^2)/360 + (11*a1^4*b1*g^2)/30240 + 
   (a1^4*g^3)/5040 + g7
newh7 = -(a1^2*a2^2*g)/720 + (a3^2*g)/12 + (a1^2*a3*b1*g)/120 + 
   (a1^4*b1^2*g)/7560 - (a1^3*b3*g)/360 - (a1*c5*g)/6 + (a1^2*a3*g^2)/40 + 
   (a1^4*b1*g^2)/1680 + (a1^4*g^3)/1120 + h7
newi7 = (a1^2*a2^2*g)/540 - (a3^2*g)/18 - (a1^2*a3*b1*g)/270 + 
   (a1^4*b1^2*g)/90720 + (a1^3*b3*g)/270 + (a1*c5*g)/18 - (a1^2*a3*g^2)/90 - 
   (a1^4*b1*g^2)/15120 - (a1^4*g^3)/3024 + i7
newj7 = -(a1*a3*b1^2*g)/360 - (a1^3*b1^3*g)/30240 - (a3*b3*g)/6 + 
   (a2*b4*g)/6 + (a1*e5*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 + j7
newk7 = (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 - (a2*b4*g)/3 - (a1*d5*g)/6 - 
   (a1*e5*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 + k7
newl7 = (a1*a2^2*b1*g)/540 - (7*a1*a3*b1^2*g)/2160 + (a1^3*b1^3*g)/90720 - 
   (a3*b3*g)/9 + (a1^2*b1*b3*g)/360 + (a2*b4*g)/36 + (b1*c5*g)/36 + 
   (a1*d5*g)/12 + (d6*g)/2 - (a1*a2^2*g^2)/270 - (a2*a3*g^2)/24 - 
   (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*b4*g^2)/72 - (c5*g^2)/36 - (a1^2*a2*g^3)/720 + 
   (a1*a3*g^3)/360 - (a1^3*b1*g^3)/1890 + (a1^3*g^4)/10080 + l7
newm7 = -(a1*a2^2*b1*g)/240 + (a1*a3*b1^2*g)/180 + (a1^3*b1^3*g)/7560 - 
   (a3*b3*g)/12 - (a1^2*b1*b3*g)/180 + (a2*b4*g)/6 - (b1*c5*g)/12 + 
   (a1*e5*g)/12 + (e6*g)/2 - (a1*a2^2*g^2)/180 + (a2*a3*g^2)/24 + 
   (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*b4*g^2)/24 + (c5*g^2)/12 + (a1^2*a2*g^3)/240 - 
   (a1*a3*g^3)/60 - (a1^3*b1*g^3)/10080 - (a1^3*g^4)/1120 + m7
newn7 = -(a1*a2^2*b1*g)/2160 - (a1*a3*b1^2*g)/1080 + (a1^3*b1^3*g)/18144 + 
   (a3*b3*g)/36 + (a1^2*b1*b3*g)/720 - (a2*b4*g)/36 - (b1*b5*g)/12 - 
   (b1*c5*g)/36 + (f6*g)/2 + (a1*a2^2*g^2)/1080 + (a1^2*a2*b1*g^2)/864 - 
   (a1*a3*b1*g^2)/1080 - (a4*b1*g^2)/24 + (a1^3*b1^2*g^2)/5040 + 
   (a1^2*b3*g^2)/540 + (a1*b4*g^2)/72 + (b5*g^2)/12 + (c5*g^2)/36 + 
   (a1*a3*g^3)/360 - (a1^3*b1*g^3)/10080 + (a1^3*g^4)/30240 + n7
newo7 = (a2^2*b1^2*g)/720 + (a1^2*b1^4*g)/30240 + (a1*b1^2*b3*g)/360 + 
   (b3^2*g)/12 - (a2*c4*g)/6 - (a1*f5*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 + o7
newp7 = (a2^2*b1^2*g)/720 - (a3*b1^3*g)/720 - (a1^2*b1^4*g)/30240 - 
   (b3^2*g)/6 + (a2*c4*g)/4 + (b1*e5*g)/12 + (a1*f5*g)/12 - 
   (a2^2*b1*g^2)/180 - (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 + 
   (b1*b4*g^2)/24 - (a1*c4*g^2)/12 - (e5*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 + (g*g6)/2 + p7
newq7 = -(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*c4*g)/12 - (b1*d5*g)/12 - 
   (b1*e5*g)/6 - (a2^2*b1*g^2)/360 + (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 - (b1*b4*g^2)/12 + (a1*c4*g^2)/24 + (d5*g^2)/12 + 
   (e5*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 + (g*h6)/2 + q7
newr7 = (a1*b1^5*g)/30240 + (b1^3*b3*g)/720 - (b1*f5*g)/12 + 
   (a2*b1^3*g^2)/1440 + (a1*b1^4*g^2)/5040 + (b1^2*b3*g^2)/180 - 
   (b1*c4*g^2)/24 + (f5*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 + (g*i6)/2 + r7
