function [ c,cnr ] = vandinterp(xdata,ydata)

% [c,cnr]=vandinterp(xdata, ydata)
% palauttaa xdata ja ydata-vektoreihin liittyvän interpolaatiopolynomin
% kertoimet korkeimman asteisesta alkaen ja Vandermonden matriisin 
% häiriöluvun ("cond nr.")
% Esim: 
%             xdata=0:5; ydata=xdata.*sin(xdata);
%             [c,cnr]=vandinterp(xdata,ydata)

xdata=xdata(:);ydata=ydata(:); % Sarakkeiksi
V=vander(xdata);
c=V\ydata;
c=c.'; % vaakavektoriksi (ei välttämätöntä, polyval toimii 
       % myös pystyv:lle.) c.' siltä varalta, että kompleksikertoimet.
cnr=cond(V);
end

