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

c0:=[w1,z1,y1,x1]:

delta0:=vector(5):
delta0[1]:=c0[1]:
for i from 2 to 4 do delta0[i]:=c0[i]-c0[i-1]: od:
delta0[5]:=1-2*c0[4]:

writeln():
write(`Delta's:`): writeln():
write(`=======`): writeln():
for i from 1 to 5 do
    write(evalf(delta0[i],DumpDigits)); writeln();
od:
writeln():
write(`C's:`): writeln():
write(`===`): writeln():
for i from 1 to 4 do
    write(evalf(c0[i],DumpDigits)); writeln();
od:
