Emphasis is put on detecting possible C rank deficiency.
User input bounds on the uncertainty C in the elements of A are used to detect numerical rank deficiency. The algorithm employs a row and column pivot strategy to C minimize the growth of uncertainty and round-off errors. All input arrays are changed on exit. NB is the C number of M by 1 right hand sides. For ease of use, RE and AE may be input as either C vectors or scalars.