< 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. >> % verifica delle formule di interpolazione di newton, e cioe' >> % nwtcoeff che calcola le differenze finite e nwtpol che calcola >> % il valore del polinomio in un punto >> % interpoliamo la funzione usata come esempio precedentemente >> f=@(x) 1./(1+25*x.*x) f = @(x) 1./(1+25*x.*x) >> x=-1:2/19:1; >> size(x) ans = 1 20 >> gg=f(x); >> plot(x,gg); >> % un po' pou' brutto il grafico, ma corretto >> % calcolo differenze successive di gg >> c=nwtcoeff(gg) c = Columns 1 through 12 0.0385 0.0091 0.0036 0.0021 0.0017 0.0018 0.0022 0.0019 -0.0112 -0.0872 0.2613 0.6558 Columns 13 through 19 -6.6904 26.4527 -73.1371 159.9380 -286.8552 413.7725 -413.7725 >> size (c) ans = 1 19 >> % salvo x(1) e il passo 2/19 >> x0=-1 x0 = -1 >> h=2/19 h = 0.1053 >> % calcolo il polinomio di interpolazione nei punti x e controllo che passi per i punti gg >> for k=1:20 pp(k)=nwtpol(x(k),x0,h,c); end >> plot(x,pp) >> % grafico piu' accurato della funzione e del polinomio >> x=-1:0.01:1; >> gg=f(x); >> size(x) ans = 1 201 >> for k=1:201 pp(k)=nwtpol(x(k),x0,h,c); end >> plot(x,gg,x,pp) >> % corrisponde al polinomio gia' tracciato con l'altro metodo. >>