%algandmed
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];

tvork=0:1e-3:3.6;

%interpoleerimine
%1) ruutfunktsioon
kordajadruut=polyfit(t,x,2); %leiame ruutfunktsiooni kordajad
xruut=polyval(kordajadruut,tvork); %leiame ruutfunktsiooni väärtused tvork 
%igas punktis

%2) kuupfunktsioon
kordajadkuup=polyfit(t,x,3);
xkuup=polyval(kordajadkuup,tvork); %leiame kuupfunktsiooni väärtused tvork
%igas punktis

%Joonis
plot(t,x,'r*',tvork,xruut,'b',tvork,xkuup,'g')
legend('algandmed','ruutlähendaja','kuuplähendaja')
xlabel('aeg (s)')
ylabel('vahemaa (dm)')
grid on

%arvutamine
%1) t=1,8
vastus1=polyval(kordajadruut,1.8)
vastus11=polyval(kordajadkuup,1.8)

%2) t=3.5
vastus2=polyval(kordajadruut,3.5)
vastus22=polyval(kordajadkuup,3.5)