Перестановка элементов матрицы относительно главной диагонали

void diag(double** a,int n)
{int i,j,k;
 double p;
  for(i=0;i< n;i++)
	 for(j=i;j< n;j++)
	 {p=a[i][j];
	  a[i][j]=a[j][i];
	  a[j][i]=p;
     }
}            

Результат работы функции: