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

read('s9odr4Eq'): readlib(realroot): 
readlib(write):
Digits:=100: DumpDigits:=20:
eps:=10^(-Digits+30):
#
#------------------------------------------------------------------
#  We like to compute the point within the unit cube
#  and minimize
#
#          max(abs(xi[i]'s) and abs(gamma[i]'s)) 
#
#  It turns out that point is not a critical point but
#  satisfies
#
#            1-2*cx[2]=cy[1]
#
   solmax:=fsolve({f[1]=0,f[2]=0, 1-2*cx[2]-cy[1]=0},{cx[1],cx[2],cy[1]},
               cx[1]=0.18 .. 0.181,cx[2]=0.124 .. 0.128,
	       cy[1]=0.7442 .. 0.7572);
   assign(solmax);
   cx0[1]:=cx[1]: cx0[2]:=cx[2]: cy0[1]:=cy[1]: cy0[2]:=1/2:
   err:=abs(f[1]) + abs(f[2]):
    if err<eps then
       appendto(s9odr4a);
       read('WriteItOut');
       close();
    else
       print(`Fatal Error: Attention Please`):
    fi:
