перемножение не квадратных матриц

[F1remaster]

Запутался в последней функции перемножения матрицы А 3х4 на матрицу транспонированную АТ 4х3. Должна получится матрица 3х3, первые 3 числа считает норм потом не совпадает.

#include <stdio.h>
#include <stdlib.h>
#include <conio.h>
#define ma 3
#define na 4
#define mc 4
#define nc 3
void at(double *A, int m, int n,double *AT);
void mult (double *A,double *AT, int m, int n,int p,double *M );
main()
{static double A[ma][na]=
{{ 2., 4., 5., 7. },
{ 9., 7., 7., 3. },
{-8., -7., -8., -4.}};
double C[mc][nc],AT[na][ma],M[ma][ma];
int i,j;
/*printf("\n Input C[%d][%d] by rows: \n",mc,nc);

for (i=0; i<mc; i++) {printf("\n");
for (j=0; j<nc; j++)
scanf("%lf", &C[i][j]); }
//Ввод матрицы  C
printf("\n C[%d][%d]: \n",mc,nc);
for (i=0; i<mc; i++){ printf("\n");
for (j=0; j<nc; j++)
printf("%11.4g",C[i][j]); } //end i, j */
printf("\n");
at(*A,na,ma,*AT);
//вывод матрицы АТ:
printf("\n AT[%d][%d]: \n",na,ma);
for (i=0; i<na; i++) {printf("\n");
for (j=0; j<ma; j++)
printf("%11.4g",AT[i][j]); } 
printf("\n");
mult (*A,*AT,ma,ma,na,*M); // ма=3 на=4
//вывод перемножения
printf("\n M[%d][%d]: \n",ma,ma);
for (i=0; i<ma; i++) {printf("\n");
for (j=0; j<ma; j++)
printf("%11.4g",M[i][j]); } 
printf("\n");
getch();
return 0; } 
// функция транспонирования матрицы А
void at(double *A,int m,int n,double *AT) //m=na=4 ; n=na=3;
{ int i,j;
for (i=0; i<n; i++) for (j=0; j<m; j++)
*(AT+j*n+i) =*(A+i*m+j);
}
//Вот тут косяк начинается
void mult(double *A, double *AT, int m, int n,int p,double *M)
{ int i,j,k,;


for (i=0; i<m; i++) //m=3
for (j=0; j<n; j++) //n=3
  {
*(M+i*n+j)=0;
for (k=0; k<p; k++) //p=4
*(M+i*n+j)+=*(A+i*n+k)*(*(AT+k*m+j));
  
  }
}

[F1remaster]

все разобрался =)

void mult(double *A, double *AT, int m, int n,int p,double *M)
{ int i,j,k,;


for (i=0; i<m; i++) //m=3
for (j=0; j<n; j++) //n=3
  {
*(M+i*n+j)=0;
for (k=0; k<p; k++) //p=4
*(M+i*n+j)+=*(A+i*p+k)*(*(AT+k*m+j));