Using finite differences to solve an equations is quick and easy, especially for ODEs. On the boundaries it is not possible to use the same scheme as in the bulk and hence the boundaries need to be treated separately.  I have read several times that, at the boundary, one can use  a  lower order finite difference approximation.…