矩阵的加法和乘法

编写一个矩阵的加法函数、一个矩阵的乘法函数，然后调用函数计算矩阵的和与积，并输出。
为了使函数更加通用，本次课程设计使用一维数组来存放矩阵。问题的关键是如何将一维数组与二维矩阵的元素之间建立对应关系。
问题拓展：
1、增加求矩阵的转置、求矩阵的秩、求正方阵的行列式，求矩阵的鞍点元素的位置与值。
2、将矩阵进行初等变换，行初等变换、列初等变换，变换成上三角矩阵或下三角矩阵。

#include<iostream>
using namespace std;
#include<cstdlib>
#include<ctime>
#include<iomanip>

void MatAdd(double a[],double b[],double c[],int n) //矩阵的加法
{
    for(int i=0;i<n;i++)
        c[i]=a[i]+b[i];
}

void MatMull(double a[],double b[],double c[],int m,int n,int p)//矩阵的乘法
{
    int i,j,k;
    for(i=0;i<m;i++)
        for(j=0;j<p;j++)
        {
            c[i*p+j]=0;
            for(k=0;k<n;k++)
                c[i*p+j]+=a[i*n+k]*b[k*p+j];
        }
}

void ShowMat(double a[],int m,int n)  //输出矩阵
{
    for(int i=0;i<m;i++)
    {
        for(int j=0;j<n;j++)
            cout<<setw(8)<<fixed<<setprecision(2)<<a[i*n+j];
        cout<<endl;
    }
    cout<<endl;
}

void MatRand(double a[],int n)  //使用随机数产生一个矩阵
{
    for(int i=0;i<n;i++)
        a[i]=(rand()%100)/10.0;
}

void fun1()  //测试矩阵的加法
{
    double a[4][5],b[4][5],c[4][5];
    MatRand(a[0],20);
    MatRand(b[0],20);
    MatAdd(a[0],b[0],c[0],20);

    ShowMat(a[0],4,5);
    ShowMat(b[0],4,5);
    ShowMat(c[0],4,5);
}

void fun2() //测试矩阵的乘法
{
    double x[4][5],y[5][8],z[4][8];
    MatRand(x[0],4*5);
    MatRand(y[0],5*8);
    MatMull(x[0],y[0],z[0],4,5,8);

    ShowMat(x[0],4,5);
    ShowMat(y[0],5,8);
    ShowMat(z[0],4,8);
}

void main()
{
    srand(time(0));
    fun1();
    fun2();
}