Add a system test to verify that analysis and code give a similar result

[valentinaschueller]

The absolute value of analysis.get_S() and process_results.compute_convergence_rate() should be (almost) equal, especially for homogeneous Dirichlet boundary conditions.