PRODUCTO DE DOS MATRICES

#include <stdlib.h>
#include <stdio.h>
# define DIM 3

int main()
{
int i,j,k;
float A[DIM][DIM], B[DIM][DIM], C[DIM][DIM];

for(k=1;k<=2;k++) // Entrada de datos
{
printf("Dame los datos de la matriz %d: n",k);
for(i=0;i<DIM;i++)
{
printf("Dame los datos de la fila %d: n",i+1);
for(j=0;j<DIM;j++)
{
if(k==1) scanf("%f", &A[i][j]);
else scanf("%f", &B[i][j]);
}
}
}
for(i=0;i<DIM;i++) // Producto de las matrices
for(j=0;j<DIM;j++)
{
C[i][j]=0;
for(k=0;k<DIM;k++)
C[i][j] += A[i][k]*B[k][j];
}
printf("La matriz resultado es: n"); // Salida de datos
for(i=0;i<DIM;i++)
{
for(j=0;j<DIM;j++)
printf("%f ",C[i][j]);
printf("n");
}
system("PAUSE");
return 0;
}