% Yksinkertainen interpolaatio
x = linspace(0,3,7);
y = cos(1+x.^2);
plot(x,y,'ro')
grid on 
p = polyfit(x,y,6)
X = linspace(0,3,100);
plot(x,y,'ro',X,polyval(p,X),'b')
% Tässähän ei ole koko totuus -- kannattaa piirtää vähän isommalla välillä
% ja katsoa mitä käy. 

