(* pu1.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

pu1series = identity
pu1temp = Normal[pu1series]
pu1 = Collect[pu1temp, Prepend[Table[c[qmax-i], {i,0,qmax}], t]]

pu1point = pu1 //. c[n_] -> 1
pu1line = pu1 //. x1->1
pu1total = pu1line //. c[n_] -> 1

Save["pu1.out", pu1, pu1point, pu1line, pu1total]
Quit[]