Program to solve non-linear boundary value problems