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

read('s15odr6Eq'): 
read('InitialCxsCys'): 
read('Newton'): 
readlib(write):

Digits:=100: DumpDigits:=20:
x0:=vector(6): 

eps:=10^(-Digits+30):
#
for j from 1 to 6 do
    x0[j]:=Cs[9,j]:
od:
err:=NewtonItern(x0):

print(`Residual error of determining equations`);
print(err);

# Write it to a file

cx0:=vector([0,x0[1],x0[2],x0[3]]):
cy0:=vector([x0[4],x0[5],x0[6],1/2]):
appendto(s15odr6); writeln():
read('WriteItOut');
writeln():
close();
