w0=1; w1=-2;w2=-1;d1=0.6;d2=-0.5;

ptA=Series[(w0/(1-d1 B)),{B,0,10}]
valA=Map[Coefficient[ptA,B,#]&,Range[0,10]]
outA=ListPlot[valA,Filling->Axis]
Export["prob_8_1_part_a.eps",outA];

ptB=Series[(((w0-w1 B))/(1-d1 B)),{B,0,10}]
valB=Map[Coefficient[ptB,B,#]&,Range[0,10]]
outB=ListPlot[valB,Filling->Axis]
Export["prob_8_1_part_b.eps",outB];

ptC=Series[(((w0-w1 B-w2 (B^2)))/(1-d1 B)),{B,0,10}]
valC=Map[Coefficient[ptC,B,#]&,Range[0,10]]
outC=ListPlot[valC,Filling->Axis]
Export["prob_8_1_part_c.eps",outC];

ptD=Series[(((w0-w1 B))/((1-d1 B-d2 (B^2)))),{B,0,10}]
valD=Map[Coefficient[ptD,B,#]&,Range[0,10]]
outD=ListPlot[valD,Filling->Axis]
Export["prob_8_1_part_d.eps",outD];