Goal:  Compute two 9 stage order 4 schemes to arbitrary precision.
====

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

Files in this directory:
=======================

   Maple Code:
   ----------

      Get_s9odr4a      ---  Driver routine to get the scheme s9odr4a

      Get_s9odr4b      ---  Driver routine to get the scheme s9odr4b

      s9odr4Eq         ---  Define determining equations.

      cx2fixed         ---  Auxiliary procedure called by get_s9odr4b

      dcx1cy1overdcx2  ---  Compute dcx1/dcx2+dcy1/dcx2, thinking
                            cx1 and cy1 as functions of cx2 defined
                            by {f[1]=0, f[2]=0}

      iterndcx1dcy1    ---  An iterative procedure for finding a zero
                            of dcx1/dcx2+dcy1/dcx2-1

      WriteItOut       ---  for dumping out solutions to s9odr4a and s9odr4b.

   Output:   s9odr4a by Get_s9odr4a,   and
   -------   s9odr4b by Get_s9odr4b

To Run: 
=======

   After Maple prompt, type   read('Get_s9odr4a');
                    or type   read('Get_s9odr4b');

Comments: The code as it is operate use 100 decimal digits, while only
          dump 20 decimal digits of final solutions. To change these,
          please go to Get_s9odr4a and Get_s9odr4b, where 

              Digits    and   DumpDigits

          are initialized.
