% Harjutustund 11, ylesanne 4
t=[0:0.4:3.6];
x=[3.1,3.8,4.9,6.5,8.4,9.9,11,11.5,11.7,11.8];
h=0.4;
% funktsioon x(t)
plot(t,x,'o')

% l2hendame andmeid kolmanda astme polynoomiga v2himruutude mottes
% kaalud on vordsed 
k=polyfit(t,x,3)
% f=x(t)=-0.4058t^3+1.6776t^2+1.6078t+2.9754
f=@(t) k(1)*t^3+k(2)*t^2+k(3)*t+k(4);

% leiame x''(t)
syms t;
t2=matlabFunction(diff(f,t,2))

% kiirenduse nullpunkt, x''(t)=0
% fsolve(funktsioon,algl2hend)
hold on
fplot(t2)
yline(0)
grid on
% algl2hend teise tuletise korral 1.5
nullpunkt_t=fsolve(t2,1.5)
nullpunkt_x=f(nullpunkt_t)
[nullpunkt_t nullpunkt_x]