clear
t=[1 2 5 8 11 12];
y=[5 8 14 26 50 98];

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

%lähendamine
ylineaar=interp1(t,y,tvork,'linear');
kordajad=polyfit(t,y,2) %saame ruutfunktsiooni kordajad
yruut=polyval(kordajad,tvork); %arvutab ruutfunktsiooni väärtused
%joonis
plot(t,y,'k*',tvork,ylineaar,tvork,yruut)
xlabel('t')
ylabel('y')
legend('lähteandmed','lineaarsplain','vähimruutude ruutfunktsioon')

%arvutamine
yvastus1=polyval(kordajad,3) %leiab ruutfunktsiooni väärtuse kui t=3
yvastus2=polyval(kordajad,9)