% Vand der Pol Epsilon
epsilon = 0.1;

% Refinements
Ns = [14, 18, 22, 26, 30, 34, 38, 42];

% Parameter for c*dt^(p/3)
c = 0.1;

diary('ode_convergence_log.txt');

methods_to_test = { '6s3k4p', '7s4k5p', '8s3k6p', '8s3k7p', '9s4k8p', '9s4k9p' };

for method=methods_to_test

	method = method{:};
	dudt = SSP_Tools.Integrators.MSRK('coefficients', method, 'initial_integrator', 'use-exact' );
	
	
	problem = SSP_Tools.TestProblems.Vanderpol('epsilon', epsilon, 'integrator', dudt);
	Test = SSP_Tools.Tests.ConvergenceODE('problem', problem, 'refinements', Ns, 't', 4.0);

	Test.run_test();
	
	% Print the test buffer to the console 
	cellfun(@fprintf, Test.output_buffer)
	results = Test.results;
	
	errors = [];
	
	for i=1:numel(results)
		errors(:,end+1) = results(i).error(:,end);
	end
	
	format long
	errors = reshape(errors, 1, 2*numel(Ns))'
	
end

diary('off');