static double[][] MatrixMultiplication(double[][] matrixOne, double[][] matrixTwo)
{
 int aRows = matrixOne.Length; int aCols = matrixOne[0].Length;
 int bRows = matrixTwo.Length; int bCols = matrixTwo[0].Length;
 if (aCols != bRows)
 throw new Exception("Out of shape matrices");
 double[][] result = CreateMatrix(aRows, bCols);
 for (int i = 0; i < aRows; ++i) // each row of MatrixOne
 for (int j = 0; j < bCols; ++j) // each col of MatrixTwo
  for (int k = 0; k < aCols; ++k)
 result[i][j] += matrixOne[i][k] * matrixTwo[k][j];
 return result;
}
static double[][] CreateMatrix(int rows, int cols)
{
 double[][] result = new double[rows][];
 for (int i = 0; i < rows; ++i)
 result[i] = new double[cols]; 
 return result;
}