Goal:  Compute four 9 stage order 6 schemes to arbitrary precision.
====

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

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

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

      Get_s9odr6abcd ---  Driver routine to get the schemes.

      s9odr6Eq       ---  Define determining equations.

      GetCritPt      ---  Compute critical points of a polynomial in
                          two variable x and y defined by Fxy. (see below.)

      Fxy            ---  Define a polynomial in two variable x and y.

      f13            ---  Define a polynomial in three variable a, y, and z.
                          It is linear in z.

      get_s9odr6b    ---  Procedure called for obtaining s9odr6a.

      zfixed         ---  Auxiliary procedure called by get_s9odr6b

      zitern         ---  Auxiliary procedure called by get_s9odr6b

      gets9odr6d     ---  Procedure called for obtaining s9odr6d.

      WriteItOut     ---  for dumping out solutions to s9odr6abcd.

      ToCheck         ---  Check residual errors of claimed solutions

   Output:  s9odr6abcd containng three schemes. The scheme 1 is the best.
   -------

   Two Solution Curves:
   --- -------- ------

      Curves.ps      ---  PostScript file displaying Two Solution Curves
                          to the determining equations for 9 stage order
                          6 schemes.

To Run: 
=======

   After Maple prompt, type   read('Get_s9odr6abcd');

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_s9odr6abcd, where 

              Digits    and   DumpDigits

          are initialized.
