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*a2*b1*d + 12*b3*d - b1^2*d^2)/24
newb4p = b4
newa5p = (360*a5 - 30*a2^2*d + 180*a4*d + 30*a3*b1*d + a1^2*b1^2*d - 
     30*a1*b3*d - 15*a2*b1*d^2 + a1*b1^2*d^2 + 30*b3*d^2 - 3*b1^2*d^3)/360
newb5p = (720*b5 + a1*b1^3*d + 60*b1*b3*d + 360*b4*d - 2*b1^3*d^2)/720
newc5p = (360*c5 + a1*b1^3*d - 60*b1*b3*d + 3*b1^3*d^2)/360
newd5p = (b1^4*d + 720*d5)/720
newa6p = (1440*a6 + 80*a4*b1*d + 2*a1*a2*b1^2*d + 40*a2*b3*d - 120*a1*b4*d + 
     720*b5*d - 4*a2*b1^2*d^2 + a1*b1^3*d^2 + 60*b1*b3*d^2 + 120*b4*d^2 - 
     2*b1^3*d^3)/1440
newb6p = (720*b6 - 60*a4*b1*d + 2*a1*a2*b1^2*d - 60*a2*b3*d + 360*c5*d + 
     6*a2*b1^2*d^2 + a1*b1^3*d^2 - 60*b1*b3*d^2 + 3*b1^3*d^3)/720
newc6p = (1440*c6 + 2*a2*b1^3*d + 120*b1*b4*d + b1^4*d^2 + 720*d*d5)/1440
newd6p = (-(b1*b4*d) + 6*d6)/6
newe6p = e6
newa7p = (30240*a7 - 5040*a2*a4*d + 126*a1*a2^2*b1*d + 2520*a5*b1*d - 
     168*a1*a3*b1^2*d - 4*a1^3*b1^3*d + 2520*a3*b3*d + 168*a1^2*b1*b3*d + 
     15120*b6*d - 2520*a1*c5*d + 168*a2^2*b1*d^2 - 1260*a4*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 + 2520*c5*d^2 + 126*a2*b1^2*d^3 + 3*a1*b1^3*d^3 - 
     504*b1*b3*d^3 + 27*b1^3*d^4)/30240
newb7p = (90720*b7 + 7560*a2*a4*d + 45360*a6*d - 84*a1*a2^2*b1*d - 
     2520*a5*b1*d + 294*a1*a3*b1^2*d - a1^3*b1^3*d - 5040*a3*b3*d - 
     252*a1^2*b1*b3*d - 7560*a1*b5*d - 252*a2^2*b1*d^2 + 2520*a4*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 - 3780*a1*b4*d^2 + 7560*b5*d^2 - 
     126*a2*b1^2*d^3 + 6*a1*b1^3*d^3 + 756*b1*b3*d^3 - 30*b1^3*d^4)/90720
newc7p = (30240*c7 + 42*a3*b1^3*d + a1^2*b1^4*d - 2520*b3^2*d + 
     2520*a2*b4*d - 2520*b1*c5*d + 15120*c6*d + 21*a2*b1^3*d^2 - 
     6*a1*b1^4*d^2 + 420*b1^2*b3*d^2 + 1260*b1*b4*d^2 - 9*b1^4*d^3 - 
     2520*a1*d*d5 + 2520*d^2*d5)/30240
newd7p = (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*b4*d + 2520*b1*b5*d + 5040*b1*c5*d + 
     9*a1*b1^4*d^2 - 756*b1^2*b3*d^2 - 2520*b1*b4*d^2 + 24*b1^4*d^3 + 
     15120*d*d6 + 30240*d7)/30240
newe7p = (-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*b4*d - 15120*b1*b5*d - 15*a1*b1^4*d^2 - 
     504*b1^2*b3*d^2 + 9*b1^4*d^3 + 90720*e7)/90720
newf7p = (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 - 5040*b1*c5*d - 18*a1*b1^4*d^2 + 756*b1^2*b3*d^2 - 
     27*b1^4*d^3 + 30240*f7)/30240
newg7p = (-(a1*b1^5*d) - 42*b1^3*b3*d + 3*b1^5*d^2 + 2520*b1*d*d5 + 
     15120*d*e6 + 30240*g7)/30240
newh7p = (a1*b1^5*d - 3*b1^5*d^2 - 5040*b1*d*d5 + 30240*h7)/30240
newi7p = (-5*a1*b1^5*d + 84*b1^3*b3*d - 6*b1^5*d^2 + 30240*i7)/30240
newj7p = (-(b1^6*d) + 30240*j7)/30240
