#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;
}