rank — rang
[i]=rank(X) [i]=rank(X,tol)
rank(X) calcule le rang "numérique" de
X c'est à dire le nombre de ses valeurs
singulières supérieures à norm(size(X),'inf') *
norm(X) * %eps.
rank(X,tol) est le nombre de valeurs singulières de
X supérieures à tol.
Notez que la valeur par défaut de tol est
proportionnelle à norm(X). Par exemple
rank([1.d-80,0;0,1.d-80]) vaut 2 !.