%Splineinterpolation, x, y är givna data, NAM-exemplet
x=[0.5 1.5 2.5 3 4 5]';
y=[3 1.5 1.5 1 1 0]';
plot(x,y,'o')
hold on
xx=[0.5:0.01:5]';
yy=spline(x,y,xx);
plot(xx,yy)
title('splineinterpolation (blå) och vanlig interpolation (röd)')
%Vanlig interpolation med 5:e gradspolynom
%A=[ones(size(x)) x x.^2 x.^3 x.^4 x.^5];
%c=A\y;
p=polyfit(x,y,5);%ger koeff i polynom p(1)*x^n+p(2)*x^(n-1)+...
xp=0.5:0.01:5;
%yp=c(1)+c(2)*xp+c(3)*xp.^2+c(4)*xp.^3+c(5)*xp.^4+c(6)*xp.^5;
yp=polyval(p,xp);%räknar ut p(1)*xp.^n+p(2)*xp.^(n-1)+...
plot(xp,yp,'r')
hold off

