function y=lagrange(xdata,xev,j)
% Lagrangen kertojafunktio ("Cardinal function")  HA maalisk. -95
% y=lagrange(xdata,xeval,j); 
% Lasketaan xdata:aan liittyv{ j:s l-funktio pisteiss{ xeval 
% esim: y=lagrange([1 2 3 4],[1 2 3 4],0) antaa 1 0 0 0
% 
j=j+1; % Matlab-vektorissa indeksi alkaa 1:st{
y=1; tuloind=1:length(xdata);ind=tuloind~=j;tuloind=tuloind(ind);
for i=tuloind;y=y.*(xev-xdata(i))./(xdata(j)-xdata(i));end;