%Hermiteinterpolation, x, y och k är givna data
x=[-1 -0.4 0 0.4 1]';
y=[0.06 0.28 1 0.28 0.06]';
k=[0.11 1.01 0 -1.01 -0.11]';
plot(x,y,'o'),hold on
h=diff(x);dy=diff(y);
g=h.*k(1:4)-dy;
c=2*dy-h.*(k(1:4)+k(2:5));
t=0:0.1:1;
for i=1:4
    xt=x(i)+t*h(i);
    yt=y(i)+t*dy(i)+t.*(1-t)*g(i)+t.^2.*(1-t)*c(i);
    plot(xt,yt)
end
title('Hermiteinterpolation (blå) och vanlig interpolation (röd)')
%Vanlig interpolation med 4:e gradspolynom
A=[ones(size(x)) x x.^2 x.^3 x.^4];
c=A\y;
xp=-1:0.01:1;
yp=c(1)+c(2)*xp+c(3)*xp.^2+c(4)*xp.^3+c(5)*xp.^4;
plot(xp,yp,'r')
hold off

