Transposicion de una matriz MxN

[mikelromero]

Buenas, estoy desarrollando un método para calcular la traspuesta de una matriz MxN. Cuando lo hago para una matriz cuadrada funciona perfectamente, pero para la matriz que quiero, que tiene que hacer 4 filas y 10 columnas no. Creo que es un problema con la alocación de memoria, pero no sé como solucionarlo.
Este es mi código:

double** transposee(double **chol,int col, int fil){

double **cholin;
cholin = (double**)malloc(fil * sizeof(*cholin));
for(int m = 0; m < fil; m++)
{
cholin[m] = (double*)malloc(col * sizeof(**cholin));
}

/* Initialisation */
for (int i=0; i < fil; i++){
for (int j=0; j < col; j++){
cholin[i][j] = chol[j][i];
}
}

return cholin;

}