To run the code, you need to have MatLab installed. (https://www.mathworks.com/)

The code was tested on Windows 10, Intel(R) Core(TM) i7-5500 CPU @ 2.40GHz MatLab Version R2019a-Update 5.

You can run the code by opening MatLab, running 'init.m' once and then one of the following programs:
	- Temporal_Error_Analysis.m (error analysis of the four different methods shown in the paper.)
   Error analysis of a single method, using a finer mesh of the same method as a reference solution is performed by
	- Error_ExpEuler.m
	- Error_2ndOrder.m
	- Error_3rdOrder.m
	- Error_4thOrder.m
	- Error_3rdOrderHeun.m (an additional 3rd-order scheme using the Butcher tableau for order three methods
				shown in the paper with c_2 = 1/3)
The code creates figures of the initial data, the function on the domain, the evolution on the boundary as well as the
mentioned error plots.
The file "saveMatrices.mat" contains the mass matrices and stiffness matrices for the one- and two-dimensional Laplace operator.

Additional png-pictures (created by the above code) are in the folder /Figures/.


(Johannes Wiedemann - 2019-09-30)