Note: This post is part of my scientific computing course (UW AMATH481) in autumn 2016. The course web page is here.
The code that you wrote for the project solves ODEs with Dirichlet boundary conditions, but restricted to the case of homogeneous (zero) boundary values. You may find that this is a pain when it comes to generating an interesting looking example. Luckily, there is a simple trick — essentially, adding a linear function to your unknown — which allows you to use the same solver to solve the same ODE with general Dirichlet boundary conditions. You can find a write up and example code for this here. Note that the example code will require your matrix and right hand side builders from the project to run.
Email: [my last name] [at] uw [dot] edu
Copyright © 2017 Travis Askham.