x=[1 2 5 8 11 12]; %x vektor
y=[5 8 14 26 50 98]; %y vektor

%teeme võrgu
xvork=1:10^-3:12;

%lähendamine splainiga
ylineaar=interp1(x,y,xvork,'linear');
ysplain=interp1(x,y,xvork,'spline');

%lähendamine vähimruutude mõttes
kordajadruut=polyfit(x,y,2) %saame ruutfunktsiooni kordajad
yruut=polyval(kordajadruut,xvork);
kordajadkuup=polyfit(x,y,3) %saame kuupfunktsiooni kordajad
ykuup=polyval(kordajadkuup,xvork);

%joonis
plot(x,y,'r*',xvork, ylineaar,xvork,ysplain,xvork,yruut,xvork,ykuup)
legend('algandmed','lineaarsplain','kuupsplain', 'ruutfunktsioon','kuupfunktsioon')

%arvutamine
vastus1=interp1(x,y,10,'linear'); %lähendamine sirgega
vastus2=interp1(x,y,10,'spline'); %kuupsplaini väärtus
vastus3=polyval(kordajadruut,10); %ruutfunktsiooni väärtus
vastus4=polyval(kordajadkuup,10); %kuupfunktsiooni väärtus
vastus=[vastus1,vastus2, vastus3,vastus4]