Multiplicación

Producto escalar de vectores:
double v1[M], v2[M];
...
prod_escalar = 0.0;
for (i = 0; i < M; i++)
prod_escalar += v1*v2;
/* prod_escalar contiene el producto escalar */
producto de matrices:
double A[M][N], B[N][O], R[M][O];
for (i=0;i<M;i++)
for (j=0;j<O;j++) {
R[j] = 0.0;
for (k=0;j<N;j++)
R[j] += A[k]*B[k][j];
}