如何提高矩阵相乘的效率，马尔科夫矩阵相乘

//M是一个100X100的方阵
double[] Metrix = new double[aM.GetLength(0)];//稳态分布已核对
                for (int c = 0; c <100; c++)
                {
                    for (int d = 0; d < 100; d++)
                    {
                        if (M[0, c] != M[d, c])
                        {
                            M = multip(M, traMetrix);         //multip()是实现矩阵相乘的方法                  
                        }
                    }
                }
               

运行太慢了
要判断100x100次，还要乘好多次
每乘一次就是100x100次运算