125 SUBROUTINE slahilb( N, NRHS, A, LDA, X, LDX, B, LDB, WORK, INFO)
133 INTEGER n, nrhs, lda, ldx, ldb, info
135 REAL a(lda, n), x(ldx, nrhs), b(ldb, nrhs), work(n)
149 INTEGER nmax_exact, nmax_approx
150 parameter(nmax_exact = 6, nmax_approx = 11)
162 IF (n .LT. 0 .OR. n .GT. nmax_approx)
THEN
164 ELSE IF (nrhs .LT. 0)
THEN
166 ELSE IF (lda .LT. n)
THEN
168 ELSE IF (ldx .LT. n)
THEN
170 ELSE IF (ldb .LT. n)
THEN
173 IF (info .LT. 0)
THEN
174 CALL
xerbla(
'SLAHILB', -info)
177 IF (n .GT. nmax_exact)
THEN
199 a(i, j) =
REAL(M) / (i + j - 1)
205 CALL
slaset(
'Full', n, nrhs, 0.0,
REAL(M), b, ldb)
212 work(j) = ( ( (work(j-1)/(j-1)) * (j-1 - n) ) /(j-1) )
218 x(i, j) = (work(i)*work(j)) / (i + j - 1)