#include <iostream> |
#include <iomanip> |
using namespace std; |
void prod( double a[4][4], double b[4][4]) //计算4*4矩阵乘积函数 |
{ |
int i,j; |
double c[4][4]; |
for (i=0;i<4;i++) |
{ |
for (j=0;j<4;j++) |
{ |
c[i][j]=a[i][0]*b[0][j]+a[i][1]*b[1][j]+a[i][2]*b[2][j]+a[i][3]*b[3][j]; |
} |
} |
for (i=0;i<4;i++) |
{ |
for ( j=0;j<4;j++) |
{ |
cout<<setw(5)<<c[i][j]; |
} |
cout<<endl; |
} |
} |
int main() |
{ |
int i ,j; |
double a[4][4]; |
double b[4][4]; |
double c[4][4]; |
//输入第一个4*4矩阵 |
cout<< "请输入第一个4*4矩阵:" <<endl; |
for (i=0;i<4;i++) |
{ |
for ( j=0;j<4;j++) |
{ |
cin>>a[i][j]; |
} |
} |
//输入第二个4*4矩阵 |
cout<< "请输入第二个4*4矩阵:" <<endl; |
for (i=0;i<4;i++) |
{ |
for (j=0;j<4;j++) |
{ |
|
cin>>b[i][j]; |
} |
} |
//计算4*4矩阵乘积 |
cout<< "两个4*4矩阵乘积:" <<endl; |
prod(a,b); |
|
return 0; |
} |