newa1 = a1 + 2*d
newa3 = a3 + 2*d^3
newa5 = a5 + 2*d^5
newb5 = b5 - (a1*a3*d)/6 - (a3*d^2)/6 + (a1^2*d^3)/6 + (a1*d^4)/6
newa7 = a7 + 2*d^7
newb7 = b7 - (a1*a5*d)/6 - (a5*d^2)/6 + (a1^2*d^5)/6 + (a1*d^6)/6
newc7 = c7 + (a1^3*a3*d)/360 - (a1*b5*d)/6 + (a1^2*a3*d^2)/45 - (b5*d^2)/6 - 
   (a1^4*d^3)/360 + (7*a1*a3*d^3)/180 - (a1^3*d^4)/45 + (7*a3*d^4)/360 - 
   (7*a1^2*d^5)/180 - (7*a1*d^6)/360
newd7 = -(a3^2*d)/6 + (a1*a3*d^3)/6 - (a3*d^4)/6 + (a1*d^6)/6 + d7
newa9 = a9 + 2*d^9
newb9 = b9 - (a1*a7*d)/6 - (a7*d^2)/6 + (a1^2*d^7)/6 + (a1*d^8)/6
newc9 = c9 + (a1^3*a5*d)/360 - (a1*b7*d)/6 + (a1^2*a5*d^2)/45 - (b7*d^2)/6 + 
   (7*a1*a5*d^3)/180 + (7*a5*d^4)/360 - (a1^4*d^5)/360 - (a1^3*d^6)/45 - 
   (7*a1^2*d^7)/180 - (7*a1*d^8)/360
newd9 = (a3*a5*d)/6 - (a1*a5*d^3)/3 - (a5*d^4)/6 + (a1*a3*d^5)/6 + 
   (a3*d^6)/3 - (a1*d^8)/6 + d9
newe9 = -(a3*a5*d)/3 + (a1*a5*d^3)/6 - (a5*d^4)/6 + (a1*a3*d^5)/6 - 
   (a3*d^6)/6 + (a1*d^8)/3 + e9
newf9 = (a1^2*a3^2*d)/180 + (a3*b5*d)/6 + (a1*a3^2*d^2)/90 - 
   (a1^3*a3*d^3)/180 + (a3^2*d^3)/60 - (a1*b5*d^3)/3 + (a1^2*a3*d^4)/45 - 
   (b5*d^4)/6 + (2*a1*a3*d^5)/45 - (a1^3*d^6)/30 + (7*a3*d^6)/180 - 
   (11*a1^2*d^7)/180 - (7*a1*d^8)/180 - (a1*d*d7)/6 - (d^2*d7)/6 + f9
newg9 = (a1^2*a3^2*d)/360 - (a3*b5*d)/3 + (a1*a3^2*d^2)/30 - 
   (a1^3*a3*d^3)/360 + (a3^2*d^3)/45 + (a1*b5*d^3)/6 - (2*a1^2*a3*d^4)/45 - 
   (b5*d^4)/6 - (a1*a3*d^5)/180 + (a1^3*d^6)/90 + (7*a3*d^6)/360 - 
   (a1^2*d^7)/60 - (7*a1*d^8)/360 + g9
newh9 = -(a1^5*a3*d)/15120 + (a1^3*b5*d)/360 - (a1*c7*d)/6 - 
   (a1^4*a3*d^2)/1260 + (a1^2*b5*d^2)/45 - (c7*d^2)/6 + (a1^6*d^3)/15120 - 
   (53*a1^3*a3*d^3)/15120 + (7*a1*b5*d^3)/180 + (a1^5*d^4)/1260 - 
   (13*a1^2*a3*d^4)/1890 + (7*b5*d^4)/360 + (53*a1^4*d^5)/15120 - 
   (31*a1*a3*d^5)/5040 + (13*a1^3*d^6)/1890 - (31*a3*d^6)/15120 + 
   (31*a1^2*d^7)/5040 + (31*a1*d^8)/15120 + h9
newa11 = a11 + 2*d^11
newb11 = b11 - (a1*a9*d)/6 - (a9*d^2)/6 + (a1^2*d^9)/6 + (a1*d^10)/6
newc11 = c11 + (a1^3*a7*d)/360 - (a1*b9*d)/6 + (a1^2*a7*d^2)/45 - 
   (b9*d^2)/6 + (7*a1*a7*d^3)/180 + (7*a7*d^4)/360 - (a1^4*d^7)/360 - 
   (a1^3*d^8)/45 - (7*a1^2*d^9)/180 - (7*a1*d^10)/360
newd11 = -(a1^5*a5*d)/15120 + (a1^3*b7*d)/360 - (a1*c9*d)/6 - 
   (a1^4*a5*d^2)/1260 + (a1^2*b7*d^2)/45 - (c9*d^2)/6 - 
   (53*a1^3*a5*d^3)/15120 + (7*a1*b7*d^3)/180 - (13*a1^2*a5*d^4)/1890 + 
   (7*b7*d^4)/360 + (a1^6*d^5)/15120 - (31*a1*a5*d^5)/5040 + 
   (a1^5*d^6)/1260 - (31*a5*d^6)/15120 + (53*a1^4*d^7)/15120 + 
   (13*a1^3*d^8)/1890 + (31*a1^2*d^9)/5040 + (31*a1*d^10)/15120 + d11
newe11 = (a1^7*a3*d)/604800 - (a1^5*b5*d)/15120 + (a1^3*c7*d)/360 + 
   (a1^6*a3*d^2)/37800 - (a1^4*b5*d^2)/1260 + (a1^2*c7*d^2)/45 - 
   (a1^8*d^3)/604800 + (79*a1^5*a3*d^3)/453600 - (53*a1^3*b5*d^3)/15120 + 
   (7*a1*c7*d^3)/180 - (a1^7*d^4)/37800 + (23*a1^4*a3*d^4)/37800 - 
   (13*a1^2*b5*d^4)/1890 + (7*c7*d^4)/360 - (79*a1^6*d^5)/453600 + 
   (367*a1^3*a3*d^5)/302400 - (31*a1*b5*d^5)/5040 - (23*a1^5*d^6)/37800 + 
   (157*a1^2*a3*d^6)/113400 - (31*b5*d^6)/15120 - (367*a1^4*d^7)/302400 + 
   (127*a1*a3*d^7)/151200 - (157*a1^3*d^8)/113400 + (127*a3*d^8)/604800 - 
   (127*a1^2*d^9)/151200 - (127*a1*d^10)/604800 + e11 - (a1*d*h9)/6 - 
   (d^2*h9)/6
newf11 = (a3*a7*d)/6 - (a1*a7*d^3)/3 - (a7*d^4)/6 + (a1*a3*d^7)/6 + 
   (a3*d^8)/3 - (a1*d^10)/6 + f11
newg11 = -(a3*a7*d)/3 + (a1*a7*d^3)/6 - (a7*d^4)/6 + (a1*a3*d^7)/6 - 
   (a3*d^8)/6 + (a1*d^10)/3 + g11
newh11 = -(a5^2*d)/6 + (a1*a5*d^5)/6 - (a5*d^6)/6 + (a1*d^10)/6 + h11
newi11 = -(a3*a5*d^3)/6 + (a3^2*d^5)/6 - (a5*d^6)/6 + (a3*d^8)/6 + i11
newj11 = -(a1^2*a3*a5*d)/120 + (a5*b5*d)/6 - (2*a1*a3*a5*d^2)/45 + 
   (a1^3*a5*d^3)/90 - (7*a3*a5*d^3)/180 + (a1^2*a5*d^4)/15 - 
   (a1^3*a3*d^5)/360 + (7*a1*a5*d^5)/90 + (a1*b5*d^5)/6 - 
   (2*a1^2*a3*d^6)/45 + (7*a5*d^6)/360 + (b5*d^6)/3 - (7*a1*a3*d^7)/60 + 
   (a1^3*d^8)/45 - (7*a3*d^8)/90 + (7*a1^2*d^9)/90 + (7*a1*d^10)/120 - 
   (a1*d*d9)/6 - (d^2*d9)/6 + j11
newk11 = (a1^2*a3*a5*d)/180 - (a5*b5*d)/3 - (a3*b7*d)/3 + (a1*a3*a5*d^2)/15 - 
   (a1^3*a5*d^3)/360 + (2*a3*a5*d^3)/45 + (a1*b7*d^3)/6 - 
   (2*a1^2*a5*d^4)/45 - (b7*d^4)/6 - (a1^3*a3*d^5)/360 - (a1*a5*d^5)/180 + 
   (a1*b5*d^5)/6 - (2*a1^2*a3*d^6)/45 + (7*a5*d^6)/360 - (b5*d^6)/6 - 
   (a1*a3*d^7)/180 + (a1^3*d^8)/45 + (7*a3*d^8)/360 - (a1^2*d^9)/30 - 
   (7*a1*d^10)/180 + k11
newl11 = -(a1^2*a3*a5*d)/120 + (5*a5*b5*d)/6 + (a3*b7*d)/6 - 
   (a1*a3*a5*d^2)/10 + (a1^3*a5*d^3)/90 - (a3*a5*d^3)/15 - (a1*b7*d^3)/3 + 
   (11*a1^2*a5*d^4)/90 - (b7*d^4)/6 - (a1^3*a3*d^5)/360 + 
   (19*a1*a5*d^5)/180 - (a1*b5*d^5)/6 + (a1^2*a3*d^6)/90 + (7*a5*d^6)/360 + 
   (2*b5*d^6)/3 - (4*a1*a3*d^7)/45 - (a1^3*d^8)/30 - (7*a3*d^8)/90 + 
   (a1^2*d^9)/20 + (7*a1*d^10)/120 + l11
newm11 = (7*a1^2*a3*a5*d)/360 - (2*a5*b5*d)/3 + (11*a1*a3*a5*d^2)/90 - 
   (a1^3*a5*d^3)/60 + (a3*a5*d^3)/10 - (a1^2*a5*d^4)/10 - (a1^3*a3*d^5)/360 - 
   (11*a1*a5*d^5)/180 - (a1*b5*d^5)/6 + (a1^2*a3*d^6)/90 + (7*a5*d^6)/360 - 
   (5*b5*d^6)/6 + (2*a1*a3*d^7)/15 - (a1^3*d^8)/30 + (7*a3*d^8)/60 - 
   (31*a1^2*d^9)/180 - (49*a1*d^10)/360 - (a1*d*e9)/6 - (d^2*e9)/6 + m11
newn11 = (a1*a3^3*d)/360 - (a3^3*d^2)/90 - (a1^2*a3^2*d^3)/360 + 
   (2*a1*a3^2*d^4)/45 + (a3^2*d^5)/60 - (a1^2*a3*d^6)/30 + (a1*a3*d^7)/180 + 
   (7*a3*d^8)/360 - (a1^2*d^9)/45 - (7*a1*d^10)/360 + (a3*d*d7)/6 - 
   (a1*d^3*d7)/3 - (d^4*d7)/6 + n11
newo11 = (a1*a3^3*d)/180 + (a3^3*d^2)/30 - (a1^2*a3^2*d^3)/180 - 
   (a3*b5*d^3)/6 - (a1*a3^2*d^4)/45 + (11*a3^2*d^5)/180 - (a1^2*a3*d^6)/90 - 
   (b5*d^6)/6 - (2*a1*a3*d^7)/45 + (7*a3*d^8)/180 - (a1^2*d^9)/60 - 
   (7*a1*d^10)/180 - (a3*d*d7)/3 + (a1*d^3*d7)/6 - (d^4*d7)/6 + o11
newp11 = -(a1^4*a3^2*d)/15120 + (a1^2*a3*b5*d)/180 - (b5^2*d)/6 - 
   (a3*c7*d)/3 - (a1^3*a3^2*d^2)/840 + (a1*a3*b5*d^2)/15 + 
   (a1^5*a3*d^3)/15120 - (23*a1^2*a3^2*d^3)/3780 - (a1^3*b5*d^3)/360 + 
   (2*a3*b5*d^3)/45 + (a1*c7*d^3)/6 + (a1^4*a3*d^4)/630 - (a1*a3^2*d^4)/126 - 
   (2*a1^2*b5*d^4)/45 - (c7*d^4)/6 + (131*a1^3*a3*d^5)/15120 - 
   (a3^2*d^5)/315 - (a1*b5*d^5)/180 - (a1^5*d^6)/2520 + 
   (17*a1^2*a3*d^6)/1890 + (7*b5*d^6)/360 - (13*a1^4*d^7)/5040 + 
   (a1*a3*d^7)/5040 - (a1^3*d^8)/945 - (31*a3*d^8)/15120 + (a1^2*d^9)/336 + 
   (31*a1*d^10)/15120 + p11
newq11 = (a1^4*a3^2*d)/15120 + (b5^2*d)/3 + (a3*c7*d)/6 + 
   (a1^3*a3^2*d^2)/840 - (a1*a3*b5*d^2)/18 - (a1^5*a3*d^3)/15120 + 
   (53*a1^2*a3^2*d^3)/7560 - (a3*b5*d^3)/36 - (a1*c7*d^3)/3 - 
   (a1^4*a3*d^4)/630 + (37*a1*a3^2*d^4)/3780 + (a1^2*b5*d^4)/18 - 
   (c7*d^4)/6 - (53*a1^3*a3*d^5)/5040 + (23*a3^2*d^5)/5040 + (a1*b5*d^5)/36 + 
   (a1^5*d^6)/2520 - (4*a1^2*a3*d^6)/315 + (53*a1^4*d^7)/15120 - 
   (a1*a3*d^7)/336 + (11*a1^3*d^8)/3780 + (31*a3*d^8)/15120 - 
   (a1^2*d^9)/630 - (31*a1*d^10)/15120 - (a1^3*d*d7)/360 - (a1^2*d^2*d7)/45 - 
   (7*a1*d^3*d7)/180 - (7*d^4*d7)/360 + (a1*d*f9)/6 + (d^2*f9)/6 + q11
newr11 = -(a1^4*a3^2*d)/3024 + (a1^2*a3*b5*d)/360 - (b5^2*d)/6 - 
   (a1^3*a3^2*d^2)/315 + (a1*a3*b5*d^2)/30 + (a1^5*a3*d^3)/3024 - 
   (173*a1^2*a3^2*d^3)/15120 + (a1^3*b5*d^3)/180 + (a3*b5*d^3)/45 + 
   (a1^4*a3*d^4)/420 - (59*a1*a3^2*d^4)/3780 + (a1^2*b5*d^4)/30 + 
   (3*a1^3*a3*d^5)/560 - (19*a3^2*d^5)/2520 + (17*a1*b5*d^5)/180 + 
   (a1^5*d^6)/1260 - (a1^2*a3*d^6)/315 + (7*b5*d^6)/120 + 
   (23*a1^4*d^7)/3780 - (79*a1*a3*d^7)/5040 + (71*a1^3*d^8)/3780 - 
   (31*a3*d^8)/3024 + (13*a1^2*d^9)/560 + (31*a1*d^10)/3024 + 
   (a1^3*d*d7)/180 + (2*a1^2*d^2*d7)/45 + (7*a1*d^3*d7)/90 + (7*d^4*d7)/180 - 
   (a1*d*f9)/3 - (d^2*f9)/3 - (a1*d*g9)/6 - (d^2*g9)/6 + r11
