%% Teht. 11 (DokuT) Vandermonde-interpolaatio
%
%%  Funktio vandinterp
type vandinterp
%% Funktion käyttöesim.
%
help vandinterp
xdata=0:5; ydata=xdata.*sin(xdata);

[c,cnr]=vandinterp(xdata,ydata)
format long   
%%
c6=c(6)       % Onko c(6)=0 ? Siltä näyttää,
              %        ja niin on, sillä (0,0) on datapiste => p(0)=0.
format        % Palaa alkup. näyttötarkkuuteen.
x=linspace(0,5);
y=polyval(c,x);

plot(xdata,ydata,'o',x,y,'r')
xlim([-.1 5.1])
shg
%%
%{ 
Huom 1: Tehtävän ohjeesta poiketen kertoimia käsitellään korkeimmasta alk.
Huom 2: Funktion vandinterp käyttö on aivan sama kuin polyfit:n, paitsi
        sikäli helpompaa, että polynomin asteluku tulee automaattisesti. 
Huom 3: Suorita type polyfit.
%}
%%
% publish('H2T11R','pdf')