Сложение, вычитание и умножение над матрицами целых чисел

using System;
using System.Linq.Expressions;
 
public static class MatrixExtension
{
    public static T[,] Sum<T>(this T[,] matrix, T[,] otherMatrix)
    {
        if (matrix.GetLength(0) != otherMatrix.GetLength(0))
            throw new Exception();
 
        if (matrix.GetLength(1) != otherMatrix.GetLength(1))
            throw new Exception();
 
        T[,] result = new T[matrix.GetLength(0), matrix.GetLength(1)];
 
        for (int i = 0; i < matrix.GetLength(0); i++)
            for (int j = 0; j < matrix.GetLength(1); j++)
                result[i, j] = Sum<T>(matrix[i, j], otherMatrix[i, j]);
        
        return result;
    }
 
    public static T[,] Subtraction<T>(this T[,] matrix, T[,] otherMatrix)
    {
        if (matrix.GetLength(0) != otherMatrix.GetLength(0))
            throw new Exception();
 
        if (matrix.GetLength(1) != otherMatrix.GetLength(1))
            throw new Exception();
 
        T[,] result = new T[matrix.GetLength(0), matrix.GetLength(1)];
 
        for (int i = 0; i < matrix.GetLength(0); i++)
            for (int j = 0; j < matrix.GetLength(1); j++)
                result[i, j] = Subtraction<T>(matrix[i, j], otherMatrix[i, j]);
 
        return result;
    }
 
    public static T[,] Multiply<T>(this T[,] matrix, T[,] otherMatrix)
    {
        if (matrix.GetLength(0) != otherMatrix.GetLength(1))
            throw new Exception();
 
        T[,] result = new T[matrix.GetLength(0), otherMatrix.GetLength(1)];
 
        for (int i = 0; i < matrix.GetLength(0); i++)
            for (int j = 0; j < otherMatrix.GetLength(1); j++)
                for (int k = 0; k < matrix.GetLength(0); k++ )
                    result[i, j] = Sum<T>(result[i, j], Multiply<T>(matrix[i, k], otherMatrix[k, j]));
        return result;
    }
 
    private static T Sum<T>(T left, T right)
    {
        var p1 = Expression.Parameter(typeof(T));
        var p2 = Expression.Parameter(typeof(T));
        Func<T, T, T> func = Expression.Lambda(Expression.Add(p1, p2), p1, p2).Compile() as Func<T, T, T>;
        return func(left, right);
    }
 
    private static T Subtraction<T>(T left, T right)
    {
        var p1 = Expression.Parameter(typeof(T));
        var p2 = Expression.Parameter(typeof(T));
        Func<T, T, T> func = Expression.Lambda(Expression.Subtract(p1, p2), p1, p2).Compile() as Func<T, T, T>;
        return func(left, right);
    }
 
    private static T Multiply<T>(T left, T right)
    {
        var p1 = Expression.Parameter(typeof(T));
        var p2 = Expression.Parameter(typeof(T));
        Func<T, T, T> func = Expression.Lambda(Expression.Multiply(p1, p2), p1, p2).Compile() as Func<T, T, T>;
        return func(left, right);
    }
}
 
 
class Program
{
    static void Main(string[] args)
    {
        int[,] matrix1 = new int[,] {
            { 1, 2, 3 },
            { 4, 5, 6 },
            { 7, 8, 9 }
        };
        int[,] matrix2 = new int[,] {
            { 1, 2, 3 },
            { 4, 5, 6 },
            { 7, 8, 9 }
        };
        int[,] result = matrix1.Multiply<int>(matrix2);
        for (int i = 0; i < result.GetLength(0); i++)
        {
            for (int j = 0; j < result.GetLength(1); j++)
                Console.Write("{0, 4}", result[i, j]);
            Console.WriteLine();
        }
        Console.ReadKey(true);
    }
}