linearregression.m

clear *; close all; clc;

x = [2,3,2,3]';
y = [4, 6, 6, 8]';
plot(x,y, '*'), hold on;

axis([0,4,0,10]);

sx = sum(x);
sy = sum(y);
sxx = sum(x.^2);
sxy = sum(x.*y);
A = [size(x, 1) sx
    sx, sxx];

B = [sy sxy]';

C = A\B

xn = 0:0.1:10;
yn = C(1) + C(2)*xn;

plot(xn, yn);