(* pl1.m: S. A. Fulling, Nov. 1995 (first version summer 1994) *) Off[General::spell1] Share[] qmax = 10 max = 2(qmax+1) point = Sum[t^(2j) c[j] x1^(2j), {j,0, qmax}] + O[t]^max identity = point pl1temp = Normal[identity] pl1 = Collect[pl1temp, Prepend[Table[c[qmax-i], {i,0,qmax}], t]] pl1point = pl1 //. c[n_] -> 1 pl1line = pl1 //. x1->1 pl1total = pl1line //. c[n_] -> 1 Save["pl1.out", pl1, pl1point, pl1line, pl1total] Quit[]