% Plotter kardinalfunksjonene for et sett noder, gitt ved
% vektoren x. 

x = [-1,1,3,4];
t = linspace(min(x),max(x),101);


n = 4;
l = ones(n,101);
for i=1:n
    for j=1:n
        if i ~= j
            l(i,:)=l(i,:).*(t-x(j))/(x(i)-x(j));
        end
    end
end
plot(t,l);
grid
legend('l_1(x)','l_2(x)','l_3(x)','l_4(x)');
title('Kardinalfunksjoner')
