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

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

NumIG:=rowdim(InitialCs): 
x0:=vector(15): 

Digits:=100: DumpDigits:=20:

#  Get s31odr10a

print(`Compute s31odr10a`):
for j from 1 to 15 do
    x0[j]:=InitialCs[7,j]:
od:
err:=NewtonItern(x0): 
print(`Residual error of determining equations`);
print(err):

# Dump out the solution to s15odr8a

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

#  Get s31odr10b

print(`Compute s31odr10b`):
for j from 1 to 15 do
    x0[j]:=InitialCs[5,j]:
od:
err:=NewtonItern(x0): 
print(`Residual error of determining equations`);
print(err):

# Dump out the solution to s15odr8b

appendto(s31odr10b):
writeln():
read('WriteItOut'):
