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

%võrk
tvork=0:1e-3:3.6;

%lähendamine vähimruutude mõttes
%1) ruutfunktsioon
kordajad1=polyfit(t,x,2); %leiame ruutfunktsiooni kordajad
ruut=polyval(kordajad1,tvork); %leiame ruutfunktsiooni väärtused

%2) kuupfunktsioon 
kordajad2=polyfit(t,x,3); %leiame kuupfunktsiooni kordajad
kuup=polyval(kordajad2,tvork); %leiame kuupfunktsiooni väärtused

%lähendamine kuupsplainiga
splain=interp1(t,x,tvork,'spline'); %leiame kuupsplaini väärtused

%joonis
plot(t,x,'r*',tvork,ruut,tvork,kuup,tvork,splain)
legend('algandmed','ruutfunktsioon','kuupfunktsioon','kuupsplain')
xlabel('aeg')
ylabel('vahemaa')

%arvutamine
vastusruut1=polyval(kordajad1,1.8)
vastusruut2=polyval(kordajad1,3.5)
vastuskuup1=polyval(kordajad2,1.8)
vastuskuup2=polyval(kordajad2,3.5)
vastusplain1=interp1(t,x,1.8,'spline')
vastusplain2=interp1(t,x,3.5,'spline')