%Kubisk spline, beräkning av k-värdena
x=[-1 -0.5 0 0.5 1]';
y=[0.06 0.28 1 0.28 0.06]';
dely=diff(y);
h=0.5;
n=length(y);
d=4*ones(n,1);d(1)=2;d(n)=2;
p=ones(n-1,1);
q=ones(n-1,1);
b=(3/h)*[dely(1);dely(1:n-2)+dely(2:n-1);dely(n-1)];
k=tridia(d,p,q,b)
%fortsätt med Hermiteinterpolation, openf51.m
%eller använd Matlabs spline funktion
xx=-1:0.01:1;
yy=spline(x,y,xx);
plot(x,y,'o',xx,yy)


