pub fn matrix_mul(a: &Matrix, b: &Matrix) -> Result<Matrix>
Multiplies the a matrix by the b matrix.
a
b