/* fichero vmatriz.c*/
/* Producto de matriz por vector */
#include <stdio.h>
#define  SIZE   3
void main(void)
{
    double matriz[SIZE][SIZE];
    double vector[SIZE];
    double solucion[SIZE];
    double sum;
    int    i,j;
    
    printf("Introduzca los datos de la matriz:\n");
    for(i=0; i<SIZE ; i++) {
        for(j=0; j<SIZE; j++) {
            printf("\nElemento (%d,%d): ", (i+1), (j+1));
            scanf(" %lf", &matriz[i][j]);
        }
    }
    printf("\n\nIntroduzca los datos del vector:\n");
    for(i=0; i<SIZE ; i++) {
        printf("\nElemento %d: ", (i+1));
        scanf("%lf", &vector[i]);
    }
    for(i=0; i<SIZE; i++) {
        sum=0;
        for(j=0; j<SIZE; j++) {
            sum += matriz[i][j]*vector[j];
        }
        solucion[i] = sum;
    }
    printf("\nEl vector solucion es:\n");
    for(i=0; i<SIZE; i++) {
        printf("Elemento %d = %lf\n", i+1, solucion[i]);
    }
}