function A=laplacian(n)
%ritorna la matrice Laplaciana di dimensione n
%cioe' la matrice con -2 sulla diagonale
%e 1 sulla sopra- e sottodiagonale
%
%e.g.
% -2  1  0  0
%  1 -2  1  0
%  0  1 -2  1
%  0  0  1 -2

A=zeros(n);
for i=1:n
    for j=1:n
        if i==j
            A(i,j)=-2;
        end
        if i==j+1 || i==j-1
            A(i,j)=1;
        end
    end
end