矩阵相乘

#include usingnamespace std; constint a = 4, b = 4, c = 2; double C[a][c];

void maray_multiply(double A[][b], double B[][c]); int main() { int i, j;

double X[a][b] = {

{ 0.0,1.0,1.0,1.0 }, { 1.0,0.0,1.0,1.0 }, { 1.0,1.0,0.0,1.0 }, { 1.0,1.0,1.0,0.0 }, };

double Y[b][c] = {

{ 1.0, 10.0 }, { 2.0, 20.0 },

{ 3.0, 30.0 }, { 4.0, 40.0 } };

maray_multiply(X,Y); for (i = 0; i < a; i++) {

for (j = 0; j < c; j++) cout << C[i][j] <<\"\\"; cout <<\"\\n\"; } return 0; }

void maray_multiply(doubleA[][b], doubleB[][c]) {

int i, j, k;

for (i = 0; i < a; i++) {

for (j = 0; j < c; j++) {

C[i][j] = 0.0;

for (k = 0; k C[i][j] = C[i][j] + A[i][k] * B[k][j]; } } }