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];
kaalud=[1 1 2 2 3 3 2 2 1 1];
%xvork
xvork=-4:10^-3:5;
%vähimruutude mõttes lähendamine kaalusid arvestades
%maatriks A
A(1,1)=sum(kaalud.*x.^4);
A(1,2)=sum(kaalud.*x.^3);
A(1,3)=sum(kaalud.*x.^2);
A(2,1)=sum(kaalud.*x.^3);
A(2,2)=sum(kaalud.*x.^2);
A(2,3)=sum(kaalud.*x);
A(3,1)=sum(kaalud.*x.^2);
A(3,2)=sum(kaalud.*x);
A(3,3)=sum(kaalud);
%vektor b
b(1)=sum(kaalud.*y.*x.^2);
b(2)=sum(kaalud.*y.*x);
b(3)=sum(kaalud.*y);
%lineaarse süsteemi Ac=b' lahendamine
c=inv(A)*b'
%ruutfunktsiooniga lähendamine
yvork=polyval(c',xvork);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%võrdsete kaaludega lähendamine
kordajad=polyfit(x,y,2);
yruut=polyval(kordajad,xvork);

%joonis
plot(x,y,'r*',xvork,yruut,xvork,yvork)
legend('algandmed','ruutfunktsioon vähimruutude mõttes võrdsete kaaludega','ruutfunktsioon vähimruutude mõttes erinevate kaaludega')
title('2. ülesanne')
xlabel('x')
ylabel('y')

%arvutamine, kui x=2,25
vastus1=polyval(kordajad,2.25);
vastus2=polyval(c',2.25);
vastus=[vastus1,vastus2]