This script solves a second order differential equation with boundary conditions at the endpoints of an interval, [a, b].