bvp1d solves systems of boundary value differential equations in a single spatial variable. It is mostly compatible with the MATLAB function bvp4c, one of the main differences being that it currently ...