clear

%algandmed
aeg=[30 33 34 35 39 44 45];
viga=[1.1 1.21 1.23 1.25 1.3 1.4 1.42];
kaalud=[1 1 1 1.5 1.5 2 2];

%võrk ja kuupsplain
vork=30:10^-2:45;
ysplain=interp1(aeg,viga,vork,'spline');

%vähimruutude mõttes lineaarfunktsioon P=c_1x+c_2
A(1,1)=sum(kaalud.*aeg.^2);
A(1,2)=sum(kaalud.*aeg);
A(2,1)=sum(kaalud.*aeg);
A(2,2)=sum(kaalud);
b(1)=sum(kaalud.*aeg.*viga);
b(2)=sum(kaalud.*viga);
c=inv(A)*b'
ylineaar=polyval(c',vork);

%joonis
plot(aeg,viga,'r*',vork,ysplain,vork,ylineaar)
legend('algandmed','kuupsplain','lineaarne lähend vähimruutude mõttes')

%arvutame vea, kui töötunde 50
viga1=interp1(aeg,viga,50,'spline')
viga2=polyval(c',50)