newa1p = a1 + 2*d
newb1p = b1
newa3p = (6*a3 - a1*b1*d - b1*d^2)/6
newb3p = (6*b3 - b1^2*d)/6
newa5p = (360*a5 - 60*a1*a3*d + a1^3*b1*d - 60*a3*d^2 + 8*a1^2*b1*d^2 + 
     14*a1*b1*d^3 + 7*b1*d^4)/360
newb5p = (360*b5 - 120*a3*b1*d + a1^2*b1^2*d + 12*a1*b1^2*d^2 + 8*b1^2*d^3)/
   360
newc5p = (180*c5 + 30*a3*b1*d + a1^2*b1^2*d - 30*a1*b3*d + 2*a1*b1^2*d^2 - 
     30*b3*d^2 + 3*b1^2*d^3)/180
newd5p = (a1*b1^3*d + 60*b1*b3*d - 4*b1^3*d^2 + 360*d5)/360
newe5p = (a1*b1^3*d - 60*b1*b3*d + 6*b1^3*d^2 + 180*e5)/180
newf5p = (b1^4*d + 360*f5)/360
newa7p = (15120*a7 + 42*a1^3*a3*d - 2520*a1*a5*d - a1^5*b1*d + 
     336*a1^2*a3*d^2 - 2520*a5*d^2 - 12*a1^4*b1*d^2 + 588*a1*a3*d^3 - 
     53*a1^3*b1*d^3 + 294*a3*d^4 - 104*a1^2*b1*d^4 - 93*a1*b1*d^5 - 31*b1*d^6)
    /15120
newb7p = (15120*b7 - 2520*a3^2*d + 84*a1^2*a3*b1*d - 5040*a5*b1*d - 
     a1^4*b1^2*d + 1008*a1*a3*b1*d^2 - 18*a1^3*b1^2*d^2 + 672*a3*b1*d^3 - 
     92*a1^2*b1^2*d^3 - 120*a1*b1^2*d^4 - 48*b1^2*d^5)/15120
newc7p = (15120*c7 + 5040*a3^2*d + 2520*a5*b1*d + a1^4*b1^2*d - 
     42*a1^3*b3*d + 2520*a1*c5*d - 840*a1*a3*b1*d^2 + 18*a1^3*b1^2*d^2 - 
     336*a1^2*b3*d^2 + 2520*c5*d^2 - 420*a3*b1*d^3 + 106*a1^2*b1^2*d^3 - 
     588*a1*b3*d^3 + 148*a1*b1^2*d^4 - 294*b3*d^4 + 69*b1^2*d^5)/15120
newd7p = (-2520*a3^2*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 + 504*a1*a3*b1*d^2 - 48*a1^3*b1^2*d^2 + 
     672*a1^2*b3*d^2 - 2520*b5*d^2 - 5040*c5*d^2 + 336*a3*b1*d^3 - 
     173*a1^2*b1^2*d^3 + 1176*a1*b3*d^3 - 236*a1*b1^2*d^4 + 588*b3*d^4 - 
     114*b1^2*d^5 + 15120*d7)/15120
newe7p = (a1^3*b1^3*d + 5040*a3*b3*d - 84*a1^2*b1*b3*d + 5040*b1*c5*d + 
     24*a1^2*b1^3*d^2 - 1008*a1*b1*b3*d^2 + 88*a1*b1^3*d^3 - 672*b1*b3*d^3 + 
     72*b1^3*d^4 + 15120*e7)/15120
newf7p = (168*a1*a3*b1^2*d - 5*a1^3*b1^3*d - 10080*a3*b3*d + 
     168*a1^2*b1*b3*d - 5040*b1*b5*d - 10080*b1*c5*d + 1008*a3*b1^2*d^2 - 
     78*a1^2*b1^3*d^2 + 2016*a1*b1*b3*d^2 - 272*a1*b1^3*d^3 + 
     1344*b1*b3*d^3 - 192*b1^3*d^4 + 15120*f7)/15120
newg7p = (84*a1*a3*b1^2*d - 5*a1^3*b1^3*d + 12600*a3*b3*d - 
     126*a1^2*b1*b3*d + 7560*b1*b5*d + 2520*b1*c5*d - 2016*a3*b1^2*d^2 + 
     6*a1^2*b1^3*d^2 - 1512*a1*b1*b3*d^2 + 204*a1*b1^3*d^3 - 1008*b1*b3*d^3 + 
     144*b1^3*d^4 + 45360*g7)/45360
newh7p = (-84*a1*a3*b1^2*d - 2*a1^3*b1^3*d + 1260*a3*b3*d + 84*a1^2*b1*b3*d + 
     1260*b1*c5*d - 84*a3*b1^2*d^2 - 6*a1^2*b1^3*d^2 + 168*a1*b1*b3*d^2 - 
     15*a1*b1^3*d^3 + 252*b1*b3*d^3 - 18*b1^3*d^4 - 1260*a1*d*e5 - 
     1260*d^2*e5 + 7560*h7)/7560
newi7p = (294*a1*a3*b1^2*d - a1^3*b1^3*d - 5040*a3*b3*d - 252*a1^2*b1*b3*d - 
     2520*b1*c5*d + 504*a3*b1^2*d^2 - 24*a1^2*b1^3*d^2 - 504*a1*b1*b3*d^2 - 
     39*a1*b1^3*d^3 - 756*b1*b3*d^3 + 12*b1^3*d^4 - 7560*a1*d*d5 - 
     7560*d^2*d5 + 45360*i7)/45360
newj7p = (42*a3*b1^3*d + a1^2*b1^4*d - 2520*b3^2*d - 12*a1*b1^4*d^2 + 
     840*b1^2*b3*d^2 - 57*b1^4*d^3 - 2520*b1*d*e5 - 2520*a1*d*f5 - 
     2520*d^2*f5 + 15120*j7)/15120
newk7p = (-84*a3*b1^3*d - 5*a1^2*b1^4*d - 42*a1*b1^2*b3*d + 5040*b3^2*d + 
     18*a1*b1^4*d^2 - 1512*b1^2*b3*d^2 + 96*b1^4*d^3 + 2520*b1*d*d5 + 
     5040*b1*d*e5 + 15120*k7)/15120
newl7p = (252*a3*b1^3*d - a1^2*b1^4*d - 168*a1*b1^2*b3*d - 2520*b3^2*d - 
     30*a1*b1^4*d^2 - 1008*b1^2*b3*d^2 + 36*b1^4*d^3 - 15120*b1*d*d5 + 
     45360*l7)/45360
newm7p = (-21*a3*b1^3*d - a1^2*b1^4*d + 63*a1*b1^2*b3*d - 630*b3^2*d - 
     9*a1*b1^4*d^2 + 378*b1^2*b3*d^2 - 27*b1^4*d^3 - 1260*b1*d*e5 + 3780*m7)/
   3780
newn7p = (252*a3*b1^3*d - 5*a1^2*b1^4*d - 84*a1*b1^2*b3*d + 2520*b3^2*d - 
     24*a1*b1^4*d^2 - 504*b1^2*b3*d^2 + 12*b1^4*d^3 + 45360*n7)/45360
newo7p = (-(a1*b1^5*d) - 42*b1^3*b3*d + 6*b1^5*d^2 + 2520*b1*d*f5 + 15120*o7)/
   15120
newp7p = (a1*b1^5*d - 6*b1^5*d^2 - 5040*b1*d*f5 + 15120*p7)/15120
newq7p = (-5*a1*b1^5*d + 84*b1^3*b3*d - 12*b1^5*d^2 + 15120*q7)/15120
newr7p = (-(b1^6*d) + 15120*r7)/15120
