clear();

a=[1 2 3;4 5 6;7 8 9]


k=10; 
n=10;
D2=[];
D=[];
D(1)=1;
kappa = 10;
for k=2 : 10,
D(k)= kappa^(-1/(n-k+1));
k = k+1;
end
D=diag(D);


U=[];
V=[];
I=eye(n,n);
J=-diag(ones(n-1,1),1);
K=-diag(ones(n-1,1),-1);
H= 2*eye(10,10)+J+K;
[U,S,V]=svd(H);
A=U*D*V';
e=ones(10,1)
b=A*e;
x=linsolve(A,-b);

//Pour vérifier notre résultat on fait cond(A) et on vérifie que cela est bien égale à 10
