{
#include "TMatrixD.h"    
#include "TDecompSVD.h"
#include "TRandom.h"

    Int_t np = 5;
    Int_t ns = 50000;
    
    TMatrixD A(ns, np);
    
    TRandom exp;
    Int_t ivar;
    for (Int_t ipat = 0; ipat < ns; ipat++)
    {
        A(ipat,0) = 1.0;
        for(ivar=0; ivar < np; ivar++)
            A(ipat,ivar) = 1.0 + exp.Uniform();
    }
   
    TDecompSVD svd(A);
   
}