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

read('s9odr4Eq'):
read('cx2fixed'); 
readlib(write):
read('dcx1cy1overdcx2');
read('iterndcx1dcy1');
#
Digits:=100: DumpDigits:=20:
eps:=10^(-Digits+30):

x0L:=vector([0.17390938084483, 0.62187148834834]);
cx20L:=0.052:
x0U:=vector([0.17412936515402, 0.62567214188176]);
cx20U:=0.056:

xM:=vector(2): 
cx2M:=NewtonItern(x0L,cx20L,x0U,cx20U,xM):
cx0[1]:=xM[1]: cx0[2]:=cx2M: cy0[1]:=xM[2]: cy0[2]:=1/2:

err:=abs(subs(cx[1]=cx0[1],cx[2]=cx0[2],cy[1]=cy0[1],f[1]))
     +abs(subs(cx[1]=cx0[1],cx[2]=cx0[2],cy[1]=cy0[1],f[2])):
if err < eps then
       appendto(s9odr4b);
       read('WriteItOut');
       close();
else
   print(`Fatal Error: Attention Please`):
fi:
