Purpose: Provide three different versions of codes for solving
=======  Lorenz Strange Attractor by Palindromic Schemes. The 
         reader can easily modified them to solve any system
         of the following kind: y'=f(y), where each component
         of f(y) is quadratic.


Subdirectories:
==============

      DOUBLE/       Codes in DOUBLE PRECISION

      DOUBLE.CS/    Codes in DOUBLE PRECISION, but simply compensated
                    summation technique is used to suppress rounding
                    errors

      REAL*16/      Codes in REAL*16. It should be runnable on
                    machines with FORTRAN that allows REAL*16, e.g.,
                    IBM's FORTRAN, Sun's FORTRAN.
