function[A] = matrice(k, N)
Delta=zeros(N,N);
Delta(1,1)=1;

for i=2:N
Delta(i,i)=k^(-(1/(N-i+1)))
end

H = diag(ones(1,N)*2)+diag(ones(1,N-1)*-1,-1)+diag(ones(1,N-1)*-1,1);
[U,S,V]=svd(H);
A=U*Delta*V';
endfunction