H3teht2

%{
Muodosta interpolaatiopolynomi pisteistölle, joka saadaan laskemalla
funktion f(x)=cos(1+x^2) arvot tasavälisessä x-pisteistössä, jossa on 7
pistettä välillä [0,3]. Piirrä samaan kuvaan funktio, datapisteet
(rinkuloilla) ja interpolaatiopolynomi. Sinun on tiedettävä, mikä on
polynomin asteluku.
Tarkistus: Kulkeeko polynomi kaikkien datapisteiden kautta.
%}
xd=linspace(0,3,7);   % xdata
yd=cos(1+xd.^2);       % ydata

x=linspace(0,3);      % Tiheä pisteistö, jossa lasketaan interp:n arvot.
kertoimet=polyfit(xd,yd,6);  % 7 pisteen kautta kulkee 6.asteen polynomi.
parvot=polyval(kertoimet,x);

plot(xd,yd,'o')     % datapisteet rinkuloilla
hold on
plot(x,parvot)
grid on
title('Polynomi-interpolaatio')