newa1p = a1 + d
newb1p = b1
newa2p = (2*a2 + b1*d)/2
newa3p = (12*a3 + 6*a2*d - a1*b1*d + b1*d^2)/12
newb3p = (12*b3 - b1^2*d)/12
newa4p = (24*a4 - 2*a1*a2*d + 12*a3*d + 2*a2*d^2 - a1*b1*d^2)/24
newb4p = (24*b4 - 2*a2*b1*d + 12*b3*d - b1^2*d^2)/24
newc4p = c4
newa5p = (720*a5 - 60*a1*a3*d + 360*a4*d + a1^3*b1*d - 30*a1*a2*d^2 + 
     60*a3*d^2 + 4*a1^2*b1*d^2 - 4*a1*b1*d^3 - b1*d^4)/720
newb5p = (720*b5 + 60*a2^2*d - 120*a3*b1*d + a1^2*b1^2*d + 6*a1*b1^2*d^2 + 
     2*b1^2*d^3)/720
newc5p = (360*c5 - 30*a2^2*d + 30*a3*b1*d + a1^2*b1^2*d - 30*a1*b3*d + 
     180*b4*d - 15*a2*b1*d^2 + a1*b1^2*d^2 + 30*b3*d^2 - 3*b1^2*d^3)/360
newd5p = (a1*b1^3*d + 60*b1*b3*d + 360*c4*d - 2*b1^3*d^2 + 720*d5)/720
newe5p = (a1*b1^3*d - 60*b1*b3*d + 3*b1^3*d^2 + 360*e5)/360
newf5p = (b1^4*d + 720*f5)/720
newa6p = (1440*a6 + 2*a1^3*a2*d - 120*a1*a4*d + 720*a5*d + 8*a1^2*a2*d^2 - 
     60*a1*a3*d^2 + 120*a4*d^2 + a1^3*b1*d^2 - 8*a1*a2*d^3 + 4*a1^2*b1*d^3 - 
     2*a2*d^4 + a1*b1*d^4)/1440
newb6p = (1440*b6 + 120*a2*a3*d - 4*a1^2*a2*b1*d - 240*a4*b1*d + 
     120*a1*b4*d - 720*c5*d + 60*a2^2*d^2 - 4*a1*a2*b1*d^2 - 60*a3*b1*d^2 - 
     2*a1^2*b1^2*d^2 + 60*a1*b3*d^2 - 120*b4*d^2 + 12*a2*b1*d^3 - 
     2*a1*b1^2*d^3 + b1^2*d^4)/1440
newc6p = (288*c6 - 24*a2*a3*d + 2*a1^2*a2*b1*d + 24*a4*b1*d - 48*a1*b4*d + 
     144*b5*d + 288*c5*d - 12*a2^2*d^2 + 4*a1*a2*b1*d^2 + a1^2*b1^2*d^2 - 
     24*a1*b3*d^2 + 48*b4*d^2 - 4*a2*b1*d^3 + 2*a1*b1^2*d^3)/288
newd6p = (2*a1*a2*b1^2*d + 40*a2*b3*d + 80*b1*b4*d - 120*a1*c4*d - 
     4*a2*b1^2*d^2 + a1*b1^3*d^2 + 60*b1*b3*d^2 + 120*c4*d^2 - 2*b1^3*d^3 + 
     720*d*d5 + 1440*d6)/1440
newe6p = (2*a1*a2*b1^2*d - 60*a2*b3*d - 60*b1*b4*d + 6*a2*b1^2*d^2 + 
     a1*b1^3*d^2 - 60*b1*b3*d^2 + 3*b1^3*d^3 + 360*d*e5 + 720*e6)/720
newf6p = (a2*b3*d - b1*b4*d + 18*f6)/18
newg6p = (2*a2*b1^3*d + 120*b1*c4*d + b1^4*d^2 + 720*d*f5 + 1440*g6)/1440
newh6p = (-(b1*c4*d) + 6*h6)/6
newi6p = i6
newa7p = (30240*a7 + 42*a1^3*a3*d - 2520*a1*a5*d + 15120*a6*d - a1^5*b1*d + 
     21*a1^3*a2*d^2 + 168*a1^2*a3*d^2 - 1260*a1*a4*d^2 + 2520*a5*d^2 - 
     6*a1^4*b1*d^2 + 84*a1^2*a2*d^3 - 168*a1*a3*d^3 - 8*a1^3*b1*d^3 + 
     21*a1*a2*d^4 - 42*a3*d^4 + 8*a1^2*b1*d^4 + 6*a1*b1*d^5 + b1*d^6)/30240
newb7p = (30240*b7 - 42*a1^2*a2^2*d - 2520*a3^2*d + 5040*a2*a4*d + 
     84*a1^2*a3*b1*d - 5040*a5*b1*d - a1^4*b1^2*d - 252*a1*a2^2*d^2 + 
     504*a1*a3*b1*d^2 - 9*a1^3*b1^2*d^2 - 84*a2^2*d^3 + 168*a3*b1*d^3 - 
     23*a1^2*b1^2*d^3 - 15*a1*b1^2*d^4 - 3*b1^2*d^5)/30240
newc7p = (30240*c7 - 42*a1^2*a2^2*d + 5040*a3^2*d - 7560*a2*a4*d + 
     2520*a5*b1*d + a1^4*b1^2*d - 42*a1^3*b3*d + 15120*b6*d + 2520*a1*c5*d + 
     168*a1*a2^2*d^2 + 1260*a2*a3*d^2 - 42*a1^2*a2*b1*d^2 - 
     420*a1*a3*b1*d^2 - 2520*a4*b1*d^2 + 9*a1^3*b1^2*d^2 - 168*a1^2*b3*d^2 + 
     1260*a1*b4*d^2 - 2520*c5*d^2 + 336*a2^2*d^3 - 42*a1*a2*b1*d^3 - 
     420*a3*b1*d^3 + 16*a1^2*b1^2*d^3 + 168*a1*b3*d^3 + 21*a2*b1*d^4 + 
     8*a1*b1^2*d^4 + 42*b3*d^4 + 3*b1^2*d^5)/30240
newd7p = (84*a1^2*a2^2*d - 2520*a3^2*d + 2520*a2*a4*d + 42*a1^2*a3*b1*d - 
     5*a1^4*b1^2*d + 84*a1^3*b3*d - 2520*a1*b5*d - 5040*a1*c5*d + 
     15120*c6*d + 84*a1*a2^2*d^2 - 1260*a2*a3*d^2 + 105*a1^2*a2*b1*d^2 + 
     252*a1*a3*b1*d^2 + 1260*a4*b1*d^2 - 24*a1^3*b1^2*d^2 + 336*a1^2*b3*d^2 - 
     2520*a1*b4*d^2 + 2520*b5*d^2 + 5040*c5*d^2 - 252*a2^2*d^3 + 
     210*a1*a2*b1*d^3 + 84*a3*b1*d^3 - 17*a1^2*b1^2*d^3 - 336*a1*b3*d^3 + 
     23*a1*b1^2*d^4 - 84*b3*d^4 + 6*b1^2*d^5 + 30240*d7)/30240
newe7p = (a1^3*b1^3*d + 5040*a3*b3*d - 84*a1^2*b1*b3*d - 5040*a2*b4*d + 
     5040*b1*c5*d + 12*a1^2*b1^3*d^2 - 504*a1*b1*b3*d^2 + 22*a1*b1^3*d^3 - 
     168*b1*b3*d^3 + 9*b1^3*d^4 + 30240*e7)/30240
newf7p = (-84*a1*a2^2*b1*d + 168*a1*a3*b1^2*d - 5*a1^3*b1^3*d - 
     10080*a3*b3*d + 168*a1^2*b1*b3*d + 10080*a2*b4*d - 5040*b1*b5*d - 
     10080*b1*c5*d - 252*a2^2*b1*d^2 + 504*a3*b1^2*d^2 - 39*a1^2*b1^3*d^2 + 
     1008*a1*b1*b3*d^2 - 68*a1*b1^3*d^3 + 336*b1*b3*d^3 - 24*b1^3*d^4 + 
     30240*f7)/30240
newg7p = (-42*a1*a2^2*b1*d + 84*a1*a3*b1^2*d - 5*a1^3*b1^3*d + 
     12600*a3*b3*d - 126*a1^2*b1*b3*d - 7560*a2*b4*d + 7560*b1*b5*d + 
     2520*b1*c5*d + 504*a2^2*b1*d^2 - 1008*a3*b1^2*d^2 + 3*a1^2*b1^3*d^2 + 
     2520*a2*b3*d^2 - 756*a1*b1*b3*d^2 - 2520*b1*b4*d^2 + 51*a1*b1^3*d^3 - 
     252*b1*b3*d^3 + 18*b1^3*d^4 + 45360*d*f6 + 90720*g7)/90720
newh7p = (126*a1*a2^2*b1*d - 168*a1*a3*b1^2*d - 4*a1^3*b1^3*d + 
     2520*a3*b3*d + 168*a1^2*b1*b3*d - 5040*a2*b4*d + 2520*b1*c5*d + 
     168*a2^2*b1*d^2 + 42*a1*a2*b1^2*d^2 - 84*a3*b1^2*d^2 - 6*a1^2*b1^3*d^2 - 
     1260*a2*b3*d^2 + 168*a1*b1*b3*d^2 - 1260*b1*b4*d^2 + 126*a2*b1^2*d^3 + 
     3*a1*b1^3*d^3 - 504*b1*b3*d^3 + 27*b1^3*d^4 - 2520*a1*d*e5 + 
     2520*d^2*e5 + 15120*d*e6 + 30240*h7)/30240
newi7p = (-84*a1*a2^2*b1*d + 294*a1*a3*b1^2*d - a1^3*b1^3*d - 5040*a3*b3*d - 
     252*a1^2*b1*b3*d + 7560*a2*b4*d - 2520*b1*c5*d - 252*a2^2*b1*d^2 + 
     63*a1*a2*b1^2*d^2 + 252*a3*b1^2*d^2 - 12*a1^2*b1^3*d^2 + 
     1260*a2*b3*d^2 - 252*a1*b1*b3*d^2 + 2520*b1*b4*d^2 - 3780*a1*c4*d^2 - 
     126*a2*b1^2*d^3 + 6*a1*b1^3*d^3 + 756*b1*b3*d^3 - 30*b1^3*d^4 - 
     7560*a1*d*d5 + 7560*d^2*d5 + 45360*d*d6 + 90720*i7)/90720
newj7p = (42*a3*b1^3*d + a1^2*b1^4*d - 2520*b3^2*d + 2520*a2*c4*d + 
     21*a2*b1^3*d^2 - 6*a1*b1^4*d^2 + 420*b1^2*b3*d^2 + 1260*b1*c4*d^2 - 
     9*b1^4*d^3 - 2520*b1*d*e5 - 2520*a1*d*f5 + 2520*d^2*f5 + 15120*d*g6 + 
     30240*j7)/30240
newk7p = (42*a2^2*b1^2*d - 84*a3*b1^3*d - 5*a1^2*b1^4*d - 42*a1*b1^2*b3*d + 
     5040*b3^2*d - 7560*a2*c4*d + 9*a1*b1^4*d^2 - 756*b1^2*b3*d^2 - 
     2520*b1*c4*d^2 + 24*b1^4*d^3 + 2520*b1*d*d5 + 5040*b1*d*e5 + 
     15120*d*h6 + 30240*k7)/30240
newl7p = (-126*a2^2*b1^2*d + 252*a3*b1^3*d - a1^2*b1^4*d - 168*a1*b1^2*b3*d - 
     2520*b3^2*d + 15120*a2*c4*d - 15*a1*b1^4*d^2 - 504*b1^2*b3*d^2 + 
     9*b1^4*d^3 - 15120*b1*d*d5 + 90720*l7)/90720
newm7p = (42*a2^2*b1^2*d - 84*a3*b1^3*d - 4*a1^2*b1^4*d + 252*a1*b1^2*b3*d - 
     2520*b3^2*d - 18*a1*b1^4*d^2 + 756*b1^2*b3*d^2 - 27*b1^4*d^3 - 
     5040*b1*d*e5 + 30240*m7)/30240
newn7p = (-126*a2^2*b1^2*d + 252*a3*b1^3*d - 5*a1^2*b1^4*d - 
     84*a1*b1^2*b3*d + 2520*b3^2*d - 12*a1*b1^4*d^2 - 252*b1^2*b3*d^2 + 
     3*b1^4*d^3 + 90720*n7)/90720
newo7p = (-(a1*b1^5*d) - 42*b1^3*b3*d + 3*b1^5*d^2 + 2520*b1*d*f5 + 
     15120*d*i6 + 30240*o7)/30240
newp7p = (a1*b1^5*d - 3*b1^5*d^2 - 5040*b1*d*f5 + 30240*p7)/30240
newq7p = (-5*a1*b1^5*d + 84*b1^3*b3*d - 6*b1^5*d^2 + 30240*q7)/30240
newr7p = (-(b1^6*d) + 30240*r7)/30240
