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)
subroutine slaset(UPLO, M, N, ALPHA, BETA, A, LDA)
SLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values...
subroutine xerbla(SRNAME, INFO)
XERBLA
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine slahilb(N, NRHS, A, LDA, X, LDX, B, LDB, WORK, INFO)
SLAHILB