Транспонировать квадратную матрицу

#include <stdio.h>
#include <stdlib.h>
double** new_matrix(unsigned rows, unsigned cols)
{
  double** dest = new double*[rows];
  for (int row = 0; row != rows; row++)
  {
    dest[row] = new double[cols];
    for (int col = 0; col != cols; col++)
      dest[row][col] = (50 - rand() % 101) / 100.0;
  }
  return dest;
}
 
double** transp(double** src, unsigned rows, unsigned cols)
{
  double** dest = new double*[cols];
  for (int row = 0; row != cols; row++)
  {
    dest[row] = new double[rows];
    for (int col = 0; col != rows; col++)
      dest[row][col] = src[col][row];
  }
  return dest;
}
 
void print_matrix(double** src, unsigned rows, unsigned cols, const char* header)
{
  printf("\n%s [%d x %d]\n", header, rows, cols);
  for (int row = 0; row != rows; row++)
  {
    for (int col = 0; col != cols; col++) printf(" %5.2lg", src[row][col]);
    printf("\n");  
  }
}
 
void free_matrix(double** src, unsigned rows)
{
  for (int row = 0; row != rows; row++) delete []src[row];
  delete []src;
}
 
#define ROWS  4
#define COLS  7
 
int main()
{
  double** x = new_matrix(ROWS, COLS);
  print_matrix(x, ROWS, COLS, " >> source: ");
  
  double** t = transp(x, ROWS, COLS);
  print_matrix(t, COLS, ROWS, " >> dest: ");
  
  free_matrix(x, ROWS);
  free_matrix(t, COLS);
  return 0;
}