Benchmark_new_step_restart.m

% Solve problem

% this is just an example that we can pass one parameter from the code 
% in the COMSOL script. The parameter which should be passed is the time step size

% this is the time step specified 
interval_length = return_timestep();
step = interval_length/1;
a = 0;
b = interval_length;

% -------- CALL THE SOLVER --------
fem.sol=femtime(fem, ...
    'init',fem0.sol, ...
    'solcomp',{'u','v'}, ...
    'outcomp',{'u','v'}, ...
    'tlist',[ a : step : b ], ...
    'atol',{'0.0000010'}, ...
    'rtol',0.00001, ...
    'tout','tlist', ...
    'maxorder',1, ...    
    'initialstep',0.00005, ...
    'maxstep',0.0005);