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

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

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

appendto(s15odr6abc);
writeln():
close();

for ii from 1 to rowdim(Cs) do
    for i from 1 to 6 do
        x[i]:=Cs[ii,i]:
    od:
    NewtonItern(x):

    ResiErr(x);

    cx0:=vector([0,x[1],x[2],x[3]]):
    cy0:=vector([x[4],x[5],x[6],1/2]):
    appendto(s15odr6abc);
    write(`Solution `); write(ii); write(`:`); writeln():	
    write(`================`); writeln():
    read('WriteItOut');
    writeln():
    close();
od:
