SUBROUTINE DLINDS(N,A,LDA,AINV,LDAINV) IMPLICIT DOUBLE PRECISION(A-H,O-Z) DIMENSION A(LDA,*),AINV(LDAINV,*) CALL DLFTDS(N,A,LDA,AINV,LDAINV) CALL DPOTRI('U',N,AINV,LDAINV,INFO) DO 30 I=1,N DO 40 J=1,I AINV(I,J)=AINV(J,I) 40 CONTINUE 30 CONTINUE RETURN END