SUBROUTINE BVPSOL(FCN,BC,IVPSOL,N,M,T,X,EPS,IOPT,INFO,IRW,RW,
     $                  IIW,IW)
C*    Begin Prologue BVPSOL
      IMPLICIT DOUBLEPRECISION(S)
      EXTERNAL FCN,BC,IVPSOL
      INTEGER N,M
      DOUBLE PRECISION T(M)
      DOUBLE PRECISION X(N,M)
      DOUBLE PRECISION EPS
      INTEGER INFO,IOPT(5)
      INTEGER IRW
      DOUBLE PRECISION RW(IRW)
      INTEGER IIW
      INTEGER IW(IIW)
C
C     ------------------------------------------------------------
C
C*  Title
C
C     (B)oundary (V)alue (P)roblem (Sol)lver for highly nonlinear
C     two point boundary value problems using a local linear
C     solver (condensing algorithm) or a global sparse linear solver
C     for the solution of the arising 
C     linear subproblems.
C
C*  Written by        P. Deuflhard, G.Bader, L. Weimann
C*  Purpose           Solution of nonlinear two-point boundary value
C                     problems.
C*  Method            Local and Global Nonlinear two-point Boundary Value
C                     Problems solver (Multiple shooting approach)
C*  Category          I1b2a - Differential and integral equations
C                             Two point boundary value problems
C*  Keywords          Nonlinear boundary value problems, Multiple
C                     shooting, Newton methods
C*  Version           1.2
C*  Revision          February 2002
C*  Latest Change     January 2004
C*  Library           CodeLib
C*  Code              Fortran 77, Double Precision
C*  Environment       Standard Fortran 77 environment on PC's,
C                     workstations and hosts.
C*  Copyright     (c) Konrad-Zuse-Zentrum fuer
C                     Informationstechnik Berlin (ZIB)
C                     Takustrasse 7, D-14195 Berlin-Dahlem
C                     phone : + 49/30/84185-0
C                     fax   : + 49/30/84185-125
C*  Contact           Lutz Weimann
C                     ZIB, Division Scientific Computing, 
C                          Department Scientific Software
C                     phone : + 49/30/84185-185
C                     fax   : + 49/30/84185-107
C                     e-mail: weimann@zib.de

The distribution contains the following files:

  bvpsol.f         the solver
  bvpsol.nrm       output of the example
  linalg_bvpsol.f  dense linear algebra (deccon,solcon)
  ma28_bvpsol.f    sparse linear algebra (ma28)
  main_bvpsol.f    main example program
  main_bvpsol.data test data for the example program
  makefile         makefile for compilation
  readme           this file
  zibconst.f       machine dependant constants
  zibsec.f         machine dependant timing routine