function matrix = ls_matrix(x,d)

% Inputs:
% 
% x: vettore 1 x N di nodi
% d: grado massimo nella matrice
%
% Output:
% 
% matrix: matrice N x (d+1) per l'approssimazione ai minimi quadrati

if d+1 > length(x)
    error('Pochi punti, o grado troppo alto!')
end

matrix = zeros(length(x),d+1);

for i=1:d+1
    
    matrix(:,i) = x.^(d+1-i)';

end