/*  File: matrixv.c
    Calculate matrix equation b = Ax */
    
#include <stdio.h>

#define M 3
#define N 3

int main()
{
    double a[M][N] = {{ 3,5,6},
                      { 4,2,1},
                      { 0,7,1}};
    double v[N] = {2,1,-2};
    double b[M];
    
    int i,j;
    
    for(i = 0;i<N;i++)
    {
        b[i] = 0;
        for(j=0;j<N;j++)
        {
            b[i] += a[i][j] * v[j];
        }
    }
    printf("b=%g,%g,%g\n",b[0],b[1],b[2]);
    return 0;
    
}