#                           Ren-Cang Li, June 1, 1996
#                           na.rcli@na-net.ornl.gov

xi0:=vector(m+1): 
xi0[1]:=cx0[1]:
for i from 2 to m do xi0[i]:=cx0[i]-cx0[i-1]: od:
xi0[m+1]:=1-2*cx0[m]:

gamma0:=vector(m):
gamma0[1]:=cy0[1]:
for i from 2 to m-1 do gamma0[i]:=cy0[i]-cy0[i-1]: od:
gamma0[m]:=1/2-cy0[m-1]:

writeln():
write(`xi's:`): writeln():
write(`====`): writeln():
for i from 1 to m+1 do 
    write(evalf(xi0[i],DumpDigits)); writeln();
od:
writeln():
write(`cx's:`): writeln():
write(`====`): writeln():
for i from 1 to m do 
    write(evalf(cx0[i],DumpDigits)); writeln();
od:
writeln():
writeln():
write(`gamma's:`): writeln():
write(`=======`): writeln():
for i from 1 to m do 
    write(evalf(gamma0[i],DumpDigits)); writeln();
od:
writeln():
write(`cy's:`): writeln():
write(`====`): writeln():
for i from 1 to m do 
    write(evalf(cy0[i],DumpDigits)); writeln();
od:
