clear();
points = [2,3;
4,1;
4,4;
5,3;
6,4];
%disp(points);
[ht,wd] = size(points);
for i=1 : ht
for j=1 : wd
if(j==1)
X(i,j)=1;
else
X(i,j)=points(i,1);
endif
endfor
endfor
disp(X);
for i=1 : ht
Y(i) = points(i,2);
endfor
Y = Y';
disp(Y);
w0=0;
w1=0;
W = [w0;w1];
disp(W);
#W = (X^t * X)^(-1) * X^t * Y
#y = kx + b, где k = w1, b = w0
M_1 = (X' * X);
M = inv(M_1);
disp(M);
#disp(M*M_1);
W = M * X' * Y;
disp(W);
figure;
x = -10:0.1:10;
#y = kx + b, где k = w1, b = w0
plot(x, W(1) + x*W(2),"markersize",15, ".");
hold on
plot(points(:,1), points(:,2),"markersize",15, ".");
hold off