Вычисление произведения квадратных матриц

/* Вычисление произведения квадратных матриц */

#include <stdio. h>

#include <malloc. h>

main()

{ int i, j,k, n;

float **a, **b, **c;

printf("\n Размер матриц? - ");

scanf("%d",&n);

a = (float**)calloc(n, sizeof(float*));

b = (float**)calloc(n, sizeof(float*));

c = (float**)calloc(n, sizeof(float*));

if(a==NULL || b==NULL || c==NULL)

{ m1: printf("\n Недостаточно памяти");

exit(1);

}

for(i=0; i<n; i++)

{ a[i] = (float*)calloc(n, sizeof(float));

b[i] = (float*)calloc(n, sizeof(float));

c[i] = (float*)calloc(n, sizeof(float));

if(a[i]==NULL || b[i]==NULL || c[i]==NULL) goto m1;

}

printf("\n Введите построчно матрицу а \n");

for(i=0; i<n; i++)

for(j=0; j<n; j++)

scanf("%f", &a[i][j]);

printf("\n Введите построчно матрицу b \n");

for(i=0; i<n; i++)

for(j=0; j<n; j++)

scanf("%f", &b[i][j]);

/* Вычисление произведения матриц */

for(i=0; i<n; i++)

for(j=0; j<n; j++)

{ c[i][j] = 0;

for(k=0; k<n; k++)

c[i][j] += a[i][k] * b[k][j];

}

/* Вывод результата */

for(i=0; i<n; i++)

{ for(j=0; j<n; j++)

printf("%f ", c[i][j] );

printf("\n");

}

}