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

read('s7odr4Eq');
readlib(write):
Digits:=100: DumpDigits:=20:
eps:=10^(-Digits+20):

with(grobner): 
q:=gbasis([a3,b3],[cx[1],cy[1]], plex):
readlib(realroot): intv:=realroot(q[2],1/1000):
cy0[1]:=fsolve(q[2],cy[1]=intv[1][1] .. intv[1][2]):
cx01:=solve(subs(cy[1]=cy0[1],q[1])):

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