function [kertoimet, kuntoisuus] = vandinterp(xdata,ydata)
% Funktio laskee vektoreihin xdata ja ydata liittyvän
% interpolaatiopolynomin kertoimet Vandermonden systeemin ratkaisulla ja
% palauttaa myös cond-luvun.
% Huomaa kertoimien MATLAB-tyylinen suuntaus 
% a(1)x^n + a(2)x^(n-1)+...+a(n-1)x +a(n)
% Käyttöesimerkki 
% xdata = 0:5; ydata = xdata.*sin(xdata);
% [c kunto] = vandinterp(xdata,ydata);
x = xdata(:); y = ydata(:);
V = vander(x);
% Opetellaan itse tekemään Vandermonde eri tehtävässä.
% Tässä nyt eleganteimmasta päästä oleva, harjoitellaan
% hiukan "maanläheisempiäkin" tapoja.
V  = bsxfun(@power,x*ones(1,length(x)),5:-1:0);
kuntoisuus = cond(V);
kertoimet = V\y;
