% vanderinterp.m % Solmukirjoitus 2/2019 Heikki Apiola xd=(1:6)'; % xdata pystyvektorina yd=[16;18;21;17;15;12]; % ydata pystyvektorina N=length(xd); ykkoset=ones(N,1); % Sama kuin [1 1 1 1 1 1]' V=[ykkoset xd xd.^ 2 xd.^ 3 xd.^ 4 xd.^ 5]; a=V\yd; % Ratkaistaan yhtälöryhmä V*a = y a' % Polynomin kertoimet kasvavien potenssien mukaan. % Tilan säästämiseksi näytetään vaakavektorina.