%interpolation med rät linje och 3:e gradspolynom
x=[-2 0 2 4]';
y=[9 2 19 15]';
%rät linje
plot(x,y,'-o')
%3:e gradspolynom
A=[ones(size(x)) x x.^2 x.^3];
c=A\y
xp=[-3:0.1:5]';
yp=c(1)+c(2)*xp+c(3)*xp.^2+c(4)*xp.^3;
hold on
plot(xp,yp)
title('Interpolation av styckvis räta linjer samt 3:e gradspolynom')
xlabel('x'),ylabel('y')
hold off