Beispiel - SAM - ETH Zürich

7 err = [err;erri]; % assemble matrix of error values
8 leg{l} = s p r i n t f(’y0 = %3.2f’,y0);
9 l = l+1;
0 end
1
2 % Plot error curves in log-log scale to discern alegebraic convergence →
Def. 1.4.5
3 f i g u r e(’name’,’erreul’);
4 ts = (tspan(2)-tspan(1))./N;
5 loglog(ts,err,’-+’); hold on;
6 loglog(ts,10*ts,’k-’);
7 x l a b e l(’{\bf timestep h}’,’fontsize’,14);
8 y l a b e l(’{\bf error (Euclidean norm)}’,’fontsize’,14);
9 leg{l} = ’O(h)’; legend(leg,’location’,’southeast’);
Numerische
Mathemtik
R. Hiptmair
rev 35327,
25. April
2011
1.4
p. 77