x=6.35:0.1:7.05;
y=[316.29 330.36 344.85 359.76 375.09 390.85 407.06 423.70];
for i=1:8
    x(i)=x(i)-0.01*alpha;
    y(i)=y(i)+0.1*beeta;
end
%interpoleerimine
%splainidega lähendamiseks on käsk interp1(x,y,xvork,'meetod')
xvork=x(1):0.0001:x(8);
lineaar=interp1(x,y,xvork,'linear');
splain=interp1(x,y,xvork,'spline');
%joonis

x0=6.726+0.01*gamma;
vastus1=interp1(x,y,x0,'linear')
vastus2=interp1(x,y,x0,'spline')
%lähendamine vähimruutudega
kordajad=polyfit(x,y,3); %otsime kuupfunktsiooni
vastus3=polyval(kordajad,x0)
ykuup=polyval(kordajad,xvork);
plot(x,y,'r*',xvork,lineaar,xvork,splain,xvork,ykuup)
legend('algandmed','lineaarsplain','kuupsplain','vähimruutude kuupfunktsioon')