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