apro

clear all
X = [ 1 2 3 4 ];
Y = [ 1 1 3 6 ];

for i = 1:4
    for j = 1:3
        A(i,j) = X(i) ^ (j-1);
    end
end

disp(A);


%M = A' * A
%b = A' * Y'

%L = M / b'