mus a FORTRAN77 code which implements the multiple shooting method for two point boundary value problems, for linear or nonlinear cases, by Robert Mattheij and G Staarink.
The two point boundary value problem is posed on an interval [a,b]. Over this interval, the solution y(x) satisfies some first-order differential equation. At the endpoints, some conditions may be imposed on some of the values of y.
Cases in which the differential equation is of higher order are reduced to a vector of first order equations in the standard way.
The routine MUSL() is suitable for problems in which the underlying differential equation is linear, having a form such as
y' = A(x) y(x) + q(x)
with linear boundary constraints:
Ba y(a) + Bb y(b) = beta
The routine MUSN() handles the more general case where the differential equation is nonlinear:
y' = f(x,y)
and the boundary constraints may be nonlinear as well:
g ( y(a), y(b) ) = 0
The computer code and data files described and made available on this web page are distributed under the MIT license
mus is available in a FORTRAN77 version.
bvpsol, a FORTRAN77 library which solves highly nonlinear two point boundary value problems using a local linear solver (condensing algorithm) or a global sparse linear solver for the solution of the arising linear subproblems, by Peter Deuflhard, Georg Bader, Weimann.
COLNEW, a FORTRAN77 library which solves a mixed-order system of ordinary differential equations (ODE's) subject to separated, multipoint boundary conditions, using collocation at Gaussian points, by Uri Ascher and Georg Bader.
FD1D_BVP, a FORTRAN77 program which applies the finite difference method to a two point boundary value problem in one spatial dimension.
FEM1D_BVP_LINEAR, a FORTRAN77 program which applies the finite element method, with piecewise linear elements, to a two point boundary value problem in one spatial dimension, and compares the computed and exact solutions with the L2 and seminorm errors.
Robert Mattheij, G Staarink