% Harjutustund 13, ylesanne 1
% tuletiste ligikaudne leidmine
x=[1:0.1:2];
h=0.1;  % yhtlane vork
y=[4,7,2.1,15,3,5.2,6,8.5,14,15.1,17.8];
plot(x,y,'o')
%tuletis=?
tuletis(1)=(y(2)-y(1))/h % valem sammuga ette
% tuletis keskmistes solmedes
% kasutame keskmistatud valemit
for i=2:(length(x)-1)
  tuletis(i)=(y(i+1)-y(i-1))/(2*h);
end
tuletis
% tuletis viimases solmes
tuletis(11)=(y(11)-y(10))/h % valem sammuga taha

% leiame 2. j2rku tuletise keskmistes solmedes
% tuletis2=?
for i=2:(length(x)-1)
  tuletis2(i)=(y(i-1)-2*y(i)+y(i+1))/(h^2);
end
tuletis2

