As far as I see, bvpm2 only handles singularities at the start, is there a way to handle also the case where you have singularities at both sides?