#include "f2c.h"
void dlinds(int n, double *a, int lda, double *ainv, int ldainv)
#define A(I,J) a[(I) + (J)* ( lda)]
#define AINV(I,J) ainv[(I) + (J)* ( ldainv)]
{
  int i,j;
  char uplo='u';
  int info;
  
  dlftds(n,a,lda,ainv,ldainv);   
  dpotri_(&uplo, &n, ainv, &ldainv, &info);
  for (i=0; i<=n-1; i++) 
     for (j=0; j<=i-1; j++) 
        AINV(i,j)=AINV(j,i); 
}