% Tund 12, ülesanne 1
% andmed joonisel - funktsioon tabelina
x=[1:0.1:2];
y=[4,7,2.1,15,3,5.2,6,8.5,14,15.1,17.8];
h=0.1; % sammu pikkus
plot(x,y,'o')
% leiame esimest järku tuletised iga x-i korral
% y'(1), y'(1.1), ...
% leiame tuletise esimese sõlme korral
%tuletis=[]; ?
tuletis(1)=(y(2)-y(1))/h % diferentsvalem sammuga ette
% leiame tuletised vahepealsetes sõlmedes
for i=2:10
  tuletis(i)=(y(i+1)-y(i-1))/(2*h); % keskmistatud valem
end
tuletis
% leiame tuletise viimase sõlme korral
tuletis(11)=(y(11)-y(10))/h
tuletis

% leiame teise tuletise y''
% leiame tuletise sõlmedes 2 kuni 10
% jätame välja esimese ja viimase sõlme
% teinetuletis=[]; ?
for i=2:10
  teinetuletis(i)=(y(i+1)-2*y(i)+y(i-1))/(h^2);
end
teinetuletis % esimeses sõlmes tuletist ei leidnud
% Matlab lisab massiivi vektori esimesse kohta 0-i