clear
x=[-4 -3 -2 -1 0 1 2 3 4 5];
y=[2.03 3 5.2 6.8 7 9.6 11 13.6 14 16.6];
kappa=[1 1 2 2 3 3 2 2 1 1];
%xvork
xvork=-4:0.001:5;
%maatriks A
A(1,1)=sum(kappa.*x.^4);
A(1,2)=sum(kappa.*x.^3);
A(1,3)=sum(kappa.*x.^2);
A(2,1)=sum(kappa.*x.^3);
A(2,2)=sum(kappa.*x.^2);
A(2,3)=sum(kappa.*x);
A(3,1)=sum(kappa.*x.^2);
A(3,2)=sum(kappa.*x);
A(3,3)=sum(kappa);
%vektor b
b(1)=sum(kappa.*y.*x.^2);
b(2)=sum(kappa.*y.*x);
b(3)=sum(kappa.*y)
%lineaarse süsteemi Ac=b' lahendamine
c=inv(A)*b'
%ruutfunktsiooniga lähendamine
yvork=polyval(c',xvork);
%joonis
plot(x,y,'r*',xvork,yvork)
legend('algandmed','ruutfunktsioon vähimruutude mõttes')
title('3. ülesanne')
xlabel('x')
ylabel('y')