% Il fenomeno di Runge, si veda Quarteroni-Saleri, pagg. 83, 84.

# Octave 3.2.2, Thu Nov 19 12:44:07 2009 ora solare Europa occidentale <unknown@unknown>
x=-5:0.1:5;
length(x)
y=1./(1+x.^2);
plot(x,y)
axis 'equal'
plot(x,y)
nodi=linspace(-5,5,3)
nodi=linspace(-5,5,9)
f=inline('1./(1+x.^2)')
fnodi=f(nodi);
p=polyfit(nodi,fnodi,8)
length(p)
yp=polyval(p,x);
hold on
plot(x,yp,'r')
nodi=linspace(-5,5,20);
fnodi=f(nodi);
p=polyfit(nodi,fnodi,19);
yp=polyval(p,x);
plot(x,yp,'g')
plot(nodi,fnodi,'o')
quit
# Octave 3.2.2, Thu Nov 19 16:30:56 2009 ora solare Europa occidentale <unknown@unknown>

