newa1p = a1 + 2*d
newa3p = a3 + 2*d^3
newa5p = a5 + 2*d^5
newb5p = (6*b5 - a1*a3*d - a3*d^2 + a1^2*d^3 + a1*d^4)/6
newa7p = a7 + 2*d^7
newb7p = (6*b7 - a1*a5*d - a5*d^2 + a1^2*d^5 + a1*d^6)/6
newc7p = (360*c7 + a1^3*a3*d - 60*a1*b5*d + 8*a1^2*a3*d^2 - 60*b5*d^2 - 
     a1^4*d^3 + 14*a1*a3*d^3 - 8*a1^3*d^4 + 7*a3*d^4 - 14*a1^2*d^5 - 7*a1*d^6)
    /360
newd7p = (-(a3^2*d) + a1*a3*d^3 - a3*d^4 + a1*d^6 + 6*d7)/6
newa9p = a9 + 2*d^9
newb9p = (6*b9 - a1*a7*d - a7*d^2 + a1^2*d^7 + a1*d^8)/6
newc9p = (360*c9 + a1^3*a5*d - 60*a1*b7*d + 8*a1^2*a5*d^2 - 60*b7*d^2 + 
     14*a1*a5*d^3 + 7*a5*d^4 - a1^4*d^5 - 8*a1^3*d^6 - 14*a1^2*d^7 - 7*a1*d^8)
    /360
newd9p = (a3*a5*d - 2*a1*a5*d^3 - a5*d^4 + a1*a3*d^5 + 2*a3*d^6 - a1*d^8 + 
     6*d9)/6
newe9p = (-2*a3*a5*d + a1*a5*d^3 - a5*d^4 + a1*a3*d^5 - a3*d^6 + 2*a1*d^8 + 
     6*e9)/6
newf9p = (a1^2*a3^2*d + 30*a3*b5*d + 2*a1*a3^2*d^2 - a1^3*a3*d^3 + 
     3*a3^2*d^3 - 60*a1*b5*d^3 + 4*a1^2*a3*d^4 - 30*b5*d^4 + 8*a1*a3*d^5 - 
     6*a1^3*d^6 + 7*a3*d^6 - 11*a1^2*d^7 - 7*a1*d^8 - 30*a1*d*d7 - 
     30*d^2*d7 + 180*f9)/180
newg9p = (a1^2*a3^2*d - 120*a3*b5*d + 12*a1*a3^2*d^2 - a1^3*a3*d^3 + 
     8*a3^2*d^3 + 60*a1*b5*d^3 - 16*a1^2*a3*d^4 - 60*b5*d^4 - 2*a1*a3*d^5 + 
     4*a1^3*d^6 + 7*a3*d^6 - 6*a1^2*d^7 - 7*a1*d^8 + 360*g9)/360
newh9p = (-(a1^5*a3*d) + 42*a1^3*b5*d - 2520*a1*c7*d - 12*a1^4*a3*d^2 + 
     336*a1^2*b5*d^2 - 2520*c7*d^2 + a1^6*d^3 - 53*a1^3*a3*d^3 + 
     588*a1*b5*d^3 + 12*a1^5*d^4 - 104*a1^2*a3*d^4 + 294*b5*d^4 + 
     53*a1^4*d^5 - 93*a1*a3*d^5 + 104*a1^3*d^6 - 31*a3*d^6 + 93*a1^2*d^7 + 
     31*a1*d^8 + 15120*h9)/15120
newa11p = a11 + 2*d^11
newb11p = (6*b11 - a1*a9*d - a9*d^2 + a1^2*d^9 + a1*d^10)/6
newc11p = (360*c11 + a1^3*a7*d - 60*a1*b9*d + 8*a1^2*a7*d^2 - 60*b9*d^2 + 
     14*a1*a7*d^3 + 7*a7*d^4 - a1^4*d^7 - 8*a1^3*d^8 - 14*a1^2*d^9 - 7*a1*d^10
     )/360
newd11p = (-(a1^5*a5*d) + 42*a1^3*b7*d - 2520*a1*c9*d - 12*a1^4*a5*d^2 + 
     336*a1^2*b7*d^2 - 2520*c9*d^2 - 53*a1^3*a5*d^3 + 588*a1*b7*d^3 - 
     104*a1^2*a5*d^4 + 294*b7*d^4 + a1^6*d^5 - 93*a1*a5*d^5 + 12*a1^5*d^6 - 
     31*a5*d^6 + 53*a1^4*d^7 + 104*a1^3*d^8 + 93*a1^2*d^9 + 31*a1*d^10 + 
     15120*d11)/15120
newe11p = (3*a1^7*a3*d - 120*a1^5*b5*d + 5040*a1^3*c7*d + 48*a1^6*a3*d^2 - 
     1440*a1^4*b5*d^2 + 40320*a1^2*c7*d^2 - 3*a1^8*d^3 + 316*a1^5*a3*d^3 - 
     6360*a1^3*b5*d^3 + 70560*a1*c7*d^3 - 48*a1^7*d^4 + 1104*a1^4*a3*d^4 - 
     12480*a1^2*b5*d^4 + 35280*c7*d^4 - 316*a1^6*d^5 + 2202*a1^3*a3*d^5 - 
     11160*a1*b5*d^5 - 1104*a1^5*d^6 + 2512*a1^2*a3*d^6 - 3720*b5*d^6 - 
     2202*a1^4*d^7 + 1524*a1*a3*d^7 - 2512*a1^3*d^8 + 381*a3*d^8 - 
     1524*a1^2*d^9 - 381*a1*d^10 + 1814400*e11 - 302400*a1*d*h9 - 
     302400*d^2*h9)/1814400
newf11p = (a3*a7*d - 2*a1*a7*d^3 - a7*d^4 + a1*a3*d^7 + 2*a3*d^8 - a1*d^10 + 
     6*f11)/6
newg11p = (-2*a3*a7*d + a1*a7*d^3 - a7*d^4 + a1*a3*d^7 - a3*d^8 + 2*a1*d^10 + 
     6*g11)/6
newh11p = (-(a5^2*d) + a1*a5*d^5 - a5*d^6 + a1*d^10 + 6*h11)/6
newi11p = (-(a3*a5*d^3) + a3^2*d^5 - a5*d^6 + a3*d^8 + 6*i11)/6
newj11p = (-3*a1^2*a3*a5*d + 60*a5*b5*d - 16*a1*a3*a5*d^2 + 4*a1^3*a5*d^3 - 
     14*a3*a5*d^3 + 24*a1^2*a5*d^4 - a1^3*a3*d^5 + 28*a1*a5*d^5 + 
     60*a1*b5*d^5 - 16*a1^2*a3*d^6 + 7*a5*d^6 + 120*b5*d^6 - 42*a1*a3*d^7 + 
     8*a1^3*d^8 - 28*a3*d^8 + 28*a1^2*d^9 + 21*a1*d^10 - 60*a1*d*d9 - 
     60*d^2*d9 + 360*j11)/360
newk11p = (2*a1^2*a3*a5*d - 120*a5*b5*d - 120*a3*b7*d + 24*a1*a3*a5*d^2 - 
     a1^3*a5*d^3 + 16*a3*a5*d^3 + 60*a1*b7*d^3 - 16*a1^2*a5*d^4 - 60*b7*d^4 - 
     a1^3*a3*d^5 - 2*a1*a5*d^5 + 60*a1*b5*d^5 - 16*a1^2*a3*d^6 + 7*a5*d^6 - 
     60*b5*d^6 - 2*a1*a3*d^7 + 8*a1^3*d^8 + 7*a3*d^8 - 12*a1^2*d^9 - 
     14*a1*d^10 + 360*k11)/360
newl11p = (-3*a1^2*a3*a5*d + 300*a5*b5*d + 60*a3*b7*d - 36*a1*a3*a5*d^2 + 
     4*a1^3*a5*d^3 - 24*a3*a5*d^3 - 120*a1*b7*d^3 + 44*a1^2*a5*d^4 - 
     60*b7*d^4 - a1^3*a3*d^5 + 38*a1*a5*d^5 - 60*a1*b5*d^5 + 4*a1^2*a3*d^6 + 
     7*a5*d^6 + 240*b5*d^6 - 32*a1*a3*d^7 - 12*a1^3*d^8 - 28*a3*d^8 + 
     18*a1^2*d^9 + 21*a1*d^10 + 360*l11)/360
newm11p = (7*a1^2*a3*a5*d - 240*a5*b5*d + 44*a1*a3*a5*d^2 - 6*a1^3*a5*d^3 + 
     36*a3*a5*d^3 - 36*a1^2*a5*d^4 - a1^3*a3*d^5 - 22*a1*a5*d^5 - 
     60*a1*b5*d^5 + 4*a1^2*a3*d^6 + 7*a5*d^6 - 300*b5*d^6 + 48*a1*a3*d^7 - 
     12*a1^3*d^8 + 42*a3*d^8 - 62*a1^2*d^9 - 49*a1*d^10 - 60*a1*d*e9 - 
     60*d^2*e9 + 360*m11)/360
newn11p = (a1*a3^3*d - 4*a3^3*d^2 - a1^2*a3^2*d^3 + 16*a1*a3^2*d^4 + 
     6*a3^2*d^5 - 12*a1^2*a3*d^6 + 2*a1*a3*d^7 + 7*a3*d^8 - 8*a1^2*d^9 - 
     7*a1*d^10 + 60*a3*d*d7 - 120*a1*d^3*d7 - 60*d^4*d7 + 360*n11)/360
newo11p = (a1*a3^3*d + 6*a3^3*d^2 - a1^2*a3^2*d^3 - 30*a3*b5*d^3 - 
     4*a1*a3^2*d^4 + 11*a3^2*d^5 - 2*a1^2*a3*d^6 - 30*b5*d^6 - 8*a1*a3*d^7 + 
     7*a3*d^8 - 3*a1^2*d^9 - 7*a1*d^10 - 60*a3*d*d7 + 30*a1*d^3*d7 - 
     30*d^4*d7 + 180*o11)/180
newp11p = (-(a1^4*a3^2*d) + 84*a1^2*a3*b5*d - 2520*b5^2*d - 5040*a3*c7*d - 
     18*a1^3*a3^2*d^2 + 1008*a1*a3*b5*d^2 + a1^5*a3*d^3 - 92*a1^2*a3^2*d^3 - 
     42*a1^3*b5*d^3 + 672*a3*b5*d^3 + 2520*a1*c7*d^3 + 24*a1^4*a3*d^4 - 
     120*a1*a3^2*d^4 - 672*a1^2*b5*d^4 - 2520*c7*d^4 + 131*a1^3*a3*d^5 - 
     48*a3^2*d^5 - 84*a1*b5*d^5 - 6*a1^5*d^6 + 136*a1^2*a3*d^6 + 294*b5*d^6 - 
     39*a1^4*d^7 + 3*a1*a3*d^7 - 16*a1^3*d^8 - 31*a3*d^8 + 45*a1^2*d^9 + 
     31*a1*d^10 + 15120*p11)/15120
newq11p = (a1^4*a3^2*d + 5040*b5^2*d + 2520*a3*c7*d + 18*a1^3*a3^2*d^2 - 
     840*a1*a3*b5*d^2 - a1^5*a3*d^3 + 106*a1^2*a3^2*d^3 - 420*a3*b5*d^3 - 
     5040*a1*c7*d^3 - 24*a1^4*a3*d^4 + 148*a1*a3^2*d^4 + 840*a1^2*b5*d^4 - 
     2520*c7*d^4 - 159*a1^3*a3*d^5 + 69*a3^2*d^5 + 420*a1*b5*d^5 + 
     6*a1^5*d^6 - 192*a1^2*a3*d^6 + 53*a1^4*d^7 - 45*a1*a3*d^7 + 
     44*a1^3*d^8 + 31*a3*d^8 - 24*a1^2*d^9 - 31*a1*d^10 - 42*a1^3*d*d7 - 
     336*a1^2*d^2*d7 - 588*a1*d^3*d7 - 294*d^4*d7 + 2520*a1*d*f9 + 
     2520*d^2*f9 + 15120*q11)/15120
newr11p = (-5*a1^4*a3^2*d + 42*a1^2*a3*b5*d - 2520*b5^2*d - 
     48*a1^3*a3^2*d^2 + 504*a1*a3*b5*d^2 + 5*a1^5*a3*d^3 - 
     173*a1^2*a3^2*d^3 + 84*a1^3*b5*d^3 + 336*a3*b5*d^3 + 36*a1^4*a3*d^4 - 
     236*a1*a3^2*d^4 + 504*a1^2*b5*d^4 + 81*a1^3*a3*d^5 - 114*a3^2*d^5 + 
     1428*a1*b5*d^5 + 12*a1^5*d^6 - 48*a1^2*a3*d^6 + 882*b5*d^6 + 
     92*a1^4*d^7 - 237*a1*a3*d^7 + 284*a1^3*d^8 - 155*a3*d^8 + 351*a1^2*d^9 + 
     155*a1*d^10 + 84*a1^3*d*d7 + 672*a1^2*d^2*d7 + 1176*a1*d^3*d7 + 
     588*d^4*d7 - 5040*a1*d*f9 - 5040*d^2*f9 - 2520*a1*d*g9 - 2520*d^2*g9 + 
     15120*r11)/15120
