Beispiel - SAM - ETH Zürich

1 f u n c t i o n errinf = gaussintperr(fun,n)
2 % Error of polynomial interpolation in Gauss points in maximum norm
3 % fun : handle to function to be interpolated
4 % n : number of interpolation nodes
5
6 path(path,’../SupportScripts’);
7 [nodes,weights] = GaussQuad(n);
8 errinf = polyintperr(fun,nodes’);
Listing 2.3: Erzeugen der Plots für Bsp. 2.2.59
1 f u n c t i o n plotgaussintperr
2 % Error of Gaussian interpolation for solution of logistic differentialequation
3
4 rec = []; % Array for recording errors
5 k = 1;
6 f o r lambda=[2,4,7,11,16]
7 sol = @(t) 1./(1+exp(-0.5*lambda*t));
8 errs = [];
9 f o r n=1:10, errs = [errs,gaussintperr(sol,n)]; end
0 rec = [rec;errs];
1 leg{k} = s p r i n t f(’\\lambda = %d’,lambda);
2 k = k+1;
3 end
4
Numerische
Mathemtik
R. Hiptmair
rev 35327,
25. April
2011
2.2
p. 190