% Comparison between Euler and true solution

[t, yeu] = eulode(@example,[0 4],2,1);
ytrue = solution([0:1:4]);

disp('     t       y_Euler   y_true   relative error (%)')
disp([t', yeu', ytrue', (abs(ytrue-yeu)./ytrue*100)'])