1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
|
/* 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;
}
|
