function[results] = test(solvers, func, timespan, init)


results = cell(length(solvers),4);

for I = 1:length(solvers)
        odefun = solvers{I};
        [t,y] = odefun(func, timespan, init);
        results{I,1} = func2str(odefun);
        results{I,2} = t;
        results{I,3} = y;
        results{I,4} = numel(t);
   
end
end