Add more system tests

Proposal

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.

Related issues

See #