Add more system tests
Test more cases that have a known output. For example, the FitzHugh-Nagumo, and the Grey-Scott models.
How to test the implementation?
If an analytical solution is available, construct a grid function to be compared with. If not, use a reference test to know when a modification is changing the results of the solver.