#define LDK 100
void dfull(int n, double (*dx)[LDK])
{
  int i,j,n1;
   
  n1=n-1; 
  for (i=0; i<=n1-1; i++) 
     for (j=i+1; j<=n-1; j++) 
        dx[j][i]=dx[i][j];
}