< M A T L A B > Copyright 1984-2005 The MathWorks, Inc. Version 7.1.0.183 (R14) Service Pack 3 August 02, 2005 To get started, select MATLAB Help or Demos from the Help menu. >> % definiamo la funzione >> f=@(x) 1./(1+25*x.*x) f = @(x) 1./(1+25*x.*x) >> a=-1 a = -1 >> b=1 b = 1 >> n=20 n = 20 >> c=interp1(n,a,b,f) c = 1.0e+05 * 0.0000 -0.0000 -0.0002 0.0000 0.0033 -0.0000 -0.0306 0.0000 0.1717 -0.0000 -0.5858 0.0000 1.2102 -0.0000 -1.4679 0.0000 0.9560 -0.0000 -0.2567 0.0000 >> x=-1.1:0.01:1.1; >> gg=f(x); >> size(x) ans = 1 221 >> for k=1:221 pp(k)=phorn(x(k),c); end >> plot(x,pp,x,gg) >> % che brutto !!! >> x=-1.01:0.01:1.01; >> size(x) ans = 1 203 >> gg=f(x); >> pp=0; >> for k=1:203 pp(k)=phorn(x(k),c); end >> plot(x,pp,x,gg) >> % brutto anche questo ma si vede bene che il polinomio di interpolazione >> % fuori dei punti fa forti oscillazioni