/* fichero bimatriz.c */
#include <stdio.h>
#define SIZE 5
void main(void)
{
    double matriz1[SIZE][SIZE], matriz2[SIZE][SIZE];
    double solucion[SIZE][SIZE], sum;
    int dimension, i, j, k;
    printf("Introduzca la dimensión de las matrices:\n");
    printf("Dimensión: ");
    scanf("%d", &dimension);
    printf("Introduzca los elementos de la primera matriz:\n");
    for (i=0; i<dimension; i++) {
        for (j=0; j<dimension; j++) {
            printf("a(%d,%d): ", i+1, j+1);
            scanf("%lf", &matriz1[i][j] );
        }
    }
    printf("\n");
    printf("Introduzca los elementos de la segunda matriz:\n");
    for (i=0; i<dimension; i++) {
        for (j=0; j<dimension; j++) {
            printf("b(%d,%d): ", i+1, j+1);
            scanf("%lf", &matriz2[i][j] );
        }
    }
    printf("\n");
    for (i=0; i<dimension; i++) {
        for (j=0; j<dimension ; j++) {
            sum = 0.0;
            for (k=0; k<dimension; k++)
                sum += matriz1[i][k] * matriz2[k][j];
            solucion[i][j] = sum;
        }
    }
    printf("Solución:\n\n\n");
    for (i=0; i<dimension; i++) {
        for (j=0; j<dimension; j++)
            printf("%5.2lf      ", solucion[i][j]);
        printf("\n");
    }
}