Solving Boundary Value Problems in R