Parallel MPI program to solve ODEs of the form u" + r(x)u = f(x), for a parallel computing course