read('s15odr8Eq'): 
read('InitialCs'): 
read('Newton'): 
readlib(write):

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

Digits:=100: DumpDigits:=20:
#
x0:=vector(7): 
for j from 1 to 7 do
    x0[j]:=InitialCs[14,j]:
od:
print(`Start to compute s15odr8`):
err:=NewtonItern(x0):
print(`Residual error of determining equations`);
print(err):

# Dump out the solution to s15odr8

appendto(s15odr8): 
writeln():
read('WriteItOut'):
close():
