lbvp.m

A simple solver for a 2nd order linear BVP

This code solves and then plots the solution of the linear BVP
y'' + p(x)y' + q(x)y = f(x) for xL < x < xR
where
a0*y(xL) + b0*y'(xL) = c0 and a1*y(xR) + b1*y'(xR) = c1

It is not as sophisticated as bvp4c, but it allows for easy manipulation of the solution procedure and solves the problem relatively quickly.