% Veden viskositeetti- johdantoesimerkki
% Tiedosto viskositeetti.m    HA 21.8.04
% Johdantoesimerkki

clear all; close all
T=[2 5 7 15];   
V=[1.67 1.519 1.43 1.14];
plot(T,V,'xk') % Piirretään taulukkopisteet erillisinä
% x-merkkeinä värillä k, niinkuin black
xlim([0 15.5]); % Asetellaan koordinaatistoa paremmaksi.
c=polyfit(T,V,3); % Interpolaatiopolynomin kertoimet
                  % suoraan Matlabin valmiilla funktiolla.
                  % Jos halutaan piirtää polynomin kuvaaja,
x=linspace(0,15); % lasketaan sen arvot vektorin x 
                  % pisteissä komennolla:
v=polyval(c,x)    
hold on
%plot(x,v)        % Solmun kuvasta jätetään pois.

y=polyval(c,10)   % Lasketaan polynomin arvo pisteesä 10.
%  1.31           % Tämä saadaan
plot(10,y,'ro')   % Piirretään laskettu arvo punaisella
                  % (red) rinkulalla (o)

