x=[1 2.2 2.4 2.7 3.1 3.5 4.5 5];
y=[9.054 15.077 15.754 18.3 17.984 15.852 1.772 -13.042];
%Вычисление вектора коэффициентов полинома y=a1*x^3+a2*x^2+a3*x+a4
a=polyfit(x,y,3)
%Вычисление значений полиномов на интервале от 10 до 20
x1=0.9:0.5:5.1;
y1=polyval(a,x1);
%Построение графика полинома и экспериментальных точек в одной
%графической области
plot(x1,y1,'-k',x,y,'ok')
grid %Сетка