A=rand(5,5);[x,kerA]=linsolve(A,b), -inv(A)*b //x is unique
// Benchmark with other linear sparse solver:
-[A,descr,ref,mtype] = ReadHBSparse(SCI+"/modules/umfpack/examples/bcsstk24.rsa");
+[A,descr,ref,mtype] = ReadHBSparse(SCI+"/modules/umfpack/demos/bcsstk24.rsa");
b = zeros(size(A,1),1);