|
MATHEMATICA
|
sol=Solve[{A==va (z^2+z A+z C), B==vb (z^2+z A+z C), C==vc (z^2+z B+z D), D==vd (z^2+z B+z D)}, {A, B, C, D}];
S=1/(1-2 z-A-B-C-D);
vsub={va->ua-1, vb->ub-1, vc->uc-1, vd->ud-1};
Fz[z_, ua_, ub_, uc_, ud_]=Simplify[S/.sol/.vsub];
G[z_]=Simplify[Fz[z, 1, 1, 1, 0]+Fz[z, 0, 1, 1, 1]+Fz[z, 1, 0, 1, 1] +Fz[z, 1, 1, 0, 1] -Fz[z, 1, 1, 0, 0] -Fz[z, 1, 0, 1, 0]-Fz[z, 1, 0, 0, 1]-Fz[z, 0, 1, 1, 0] -Fz[z, 0, 1, 0, 1] -Fz[z, 0, 0, 1, 1]+Fz[z, 1, 0, 0, 0]+Fz[z, 0, 1, 0, 0] +Fz[z, 0, 0, 1, 0] +Fz[z, 0, 0, 0, 1] -Fz[z, 0, 0, 0, 0]];
Drop[Flatten[CoefficientList[Series[1/(1-2z)-G[z], {z, 0, 40}], z]], 5]
CoefficientList[Series[-2x^5(-2+x+2x^2)/((2x-1)(x^2+x-1)(x-1)^2), {x, 0, 50}], x] (* Harvey P. Dale, May 30 2018 *)
|