fplot( @(x) x^2 * cos(x), [-4 4],'r*' )
hold on

p= [1,2,3,4,5]
x=1:0.1:10
y=polyval(p,x)
plot(x,y)


polyfit(x,y,n) // n è il grado