newa1 = a1
newb1 = b1 + 2*g
newa3 = a3 + (a1^2*g)/6
newb3 = b3 + (a1*b1*g)/6 + (a1*g^2)/6
newa5 = a5 - (a1^4*g)/360
newb5 = b5 + (a1*a3*g)/6 - (a1^3*b1*g)/360 + (a1^3*g^2)/90
newc5 = c5 - (a1*a3*g)/3 - (a1^3*b1*g)/180 - (a1^3*g^2)/30
newd5 = d5 - (a1^2*b1^2*g)/360 - (a1*b3*g)/3 - (a1^2*b1*g^2)/30 - 
   (a1^2*g^3)/45
newe5 = e5 - (a3*b1*g)/6 - (a1^2*b1^2*g)/180 + (a1*b3*g)/6 - (a3*g^2)/6 - 
   (a1^2*b1*g^2)/90 - (a1^2*g^3)/60
newf5 = f5 - (a1*b1^3*g)/360 - (b1*b3*g)/6 - (a1*b1^2*g^2)/45 - (b3*g^2)/6 - 
   (7*a1*b1*g^3)/180 - (7*a1*g^4)/360
newa7 = a7 + (a1^6*g)/15120
newb7 = b7 - (a1^3*a3*g)/360 + (a1*a5*g)/6 + (a1^5*b1*g)/15120 - 
   (a1^5*g^2)/2520
newc7 = c7 - (a1*a5*g)/3 - (a1^5*b1*g)/15120 + (a1^5*g^2)/2520
newd7 = d7 + (a1^3*a3*g)/180 + (a1^5*b1*g)/3024 + (a1^5*g^2)/1260
newe7 = e7 + (a3^2*g)/6 - (a5*b1*g)/6 - (a1^4*b1^2*g)/15120 + 
   (a1^3*b3*g)/360 - (a1*c5*g)/6 + (a1^2*a3*g^2)/18 - (a5*g^2)/6 + 
   (a1^4*b1*g^2)/1260 + (19*a1^4*g^3)/5040
newf7 = f7 - (a3^2*g)/3 - (a1^2*a3*b1*g)/360 + (a1^4*b1^2*g)/3024 - 
   (a1^3*b3*g)/180 + (a1*b5*g)/6 + (a1*c5*g)/3 - (a1^2*a3*g^2)/10 - 
   (a1^4*b1*g^2)/840 - (2*a1^4*g^3)/315
newg7 = (a3^2*g)/9 + (a1^2*a3*b1*g)/540 + (a1^4*b1^2*g)/9072 + 
   (a1^3*b3*g)/270 - (a1*b5*g)/3 - (a1*c5*g)/9 + (a1^2*a3*g^2)/90 + 
   (11*a1^4*b1*g^2)/7560 + (a1^4*g^3)/630 + g7
newh7 = (a3^2*g)/6 + (a1^2*a3*b1*g)/60 + (a1^4*b1^2*g)/3780 - 
   (a1^3*b3*g)/180 - (a1*c5*g)/3 + (a1^2*a3*g^2)/10 + (a1^4*b1*g^2)/420 + 
   (a1^4*g^3)/140 + h7
newi7 = -(a3^2*g)/9 - (a1^2*a3*b1*g)/135 + (a1^4*b1^2*g)/45360 + 
   (a1^3*b3*g)/135 + (a1*c5*g)/9 - (2*a1^2*a3*g^2)/45 - (a1^4*b1*g^2)/3780 - 
   (a1^4*g^3)/378 + i7
newj7 = -(a1*a3*b1^2*g)/180 - (a1^3*b1^3*g)/15120 - (a3*b3*g)/3 + 
   (a1*e5*g)/3 - (a1*a3*b1*g^2)/15 - (a1^3*b1^2*g^2)/630 - (2*a1*a3*g^3)/45 - 
   (11*a1^3*b1*g^3)/1890 - (a1^3*g^4)/210 + j7
newk7 = (a1*a3*b1^2*g)/90 + (a1^3*b1^3*g)/3024 + (2*a3*b3*g)/3 + 
   (a1^2*b1*b3*g)/90 - (a1*d5*g)/3 - (2*a1*e5*g)/3 + (2*a1*a3*b1*g^2)/15 + 
   (13*a1^3*b1^2*g^2)/2520 + (a1^2*b3*g^2)/15 + (4*a1*a3*g^3)/45 + 
   (17*a1^3*b1*g^3)/945 + (4*a1^3*g^4)/315 + k7
newl7 = (-7*a1*a3*b1^2*g)/1080 + (a1^3*b1^3*g)/45360 - (2*a3*b3*g)/9 + 
   (a1^2*b1*b3*g)/180 + (b1*c5*g)/18 + (a1*d5*g)/6 - (11*a1*a3*b1*g^2)/270 - 
   (a1^3*b1^2*g^2)/2520 - (11*a1^2*b3*g^2)/270 + (c5*g^2)/18 - 
   (a1*a3*g^3)/30 - (13*a1^3*b1*g^3)/2520 - (a1^3*g^4)/252 + l7
newm7 = (a1*a3*b1^2*g)/90 + (a1^3*b1^3*g)/3780 - (a3*b3*g)/6 - 
   (a1^2*b1*b3*g)/90 - (b1*c5*g)/6 + (a1*e5*g)/6 + (a1*a3*b1*g^2)/45 + 
   (a1^3*b1^2*g^2)/1260 - (a1^2*b3*g^2)/90 - (c5*g^2)/6 + (a1*a3*g^3)/30 + 
   (a1^3*b1*g^3)/504 + (a1^3*g^4)/420 + m7
newn7 = -(a1*a3*b1^2*g)/540 + (a1^3*b1^3*g)/9072 + (a3*b3*g)/18 + 
   (a1^2*b1*b3*g)/360 - (b1*b5*g)/6 - (b1*c5*g)/18 - (a1*a3*b1*g^2)/270 + 
   (a1^3*b1^2*g^2)/1260 + (a1^2*b3*g^2)/135 - (b5*g^2)/6 - (c5*g^2)/18 - 
   (a1*a3*g^3)/180 + (23*a1^3*b1*g^3)/15120 + (a1^3*g^4)/1890 + n7
newo7 = (a1^2*b1^4*g)/15120 + (a1*b1^2*b3*g)/180 + (b3^2*g)/6 - (a1*f5*g)/3 + 
   (a1^2*b1^3*g^2)/840 + (a1*b1*b3*g^2)/15 + (23*a1^2*b1^2*g^3)/3780 + 
   (2*a1*b3*g^3)/45 + (a1^2*b1*g^4)/126 + (a1^2*g^5)/315 + o7
newp7 = -(a3*b1^3*g)/360 - (a1^2*b1^4*g)/15120 - (b3^2*g)/3 + (b1*e5*g)/6 + 
   (a1*f5*g)/6 - (a3*b1^2*g^2)/45 - (a1^2*b1^3*g^2)/840 - (a1*b1*b3*g^2)/18 + 
   (e5*g^2)/6 - (7*a3*b1*g^3)/180 - (53*a1^2*b1^2*g^3)/7560 - 
   (a1*b3*g^3)/36 - (7*a3*g^4)/360 - (37*a1^2*b1*g^4)/3780 - 
   (23*a1^2*g^5)/5040 + p7
newq7 = (a3*b1^3*g)/180 + (a1^2*b1^4*g)/3024 + (a1*b1^2*b3*g)/360 + 
   (b3^2*g)/6 - (b1*d5*g)/6 - (b1*e5*g)/3 + (2*a3*b1^2*g^2)/45 + 
   (a1^2*b1^3*g^2)/315 + (a1*b1*b3*g^2)/30 - (d5*g^2)/6 - (e5*g^2)/3 + 
   (7*a3*b1*g^3)/90 + (173*a1^2*b1^2*g^3)/15120 + (a1*b3*g^3)/45 + 
   (7*a3*g^4)/180 + (59*a1^2*b1*g^4)/3780 + (19*a1^2*g^5)/2520 + q7
newr7 = (a1*b1^5*g)/15120 + (b1^3*b3*g)/360 - (b1*f5*g)/6 + 
   (a1*b1^4*g^2)/1260 + (b1^2*b3*g^2)/45 - (f5*g^2)/6 + 
   (53*a1*b1^3*g^3)/15120 + (7*b1*b3*g^3)/180 + (13*a1*b1^2*g^4)/1890 + 
   (7*b3*g^4)/360 + (31*a1*b1*g^5)/5040 + (31*a1*g^6)/15120 + r7
