diff --git a/LinearAlgebra/Matrix.cpp b/LinearAlgebra/Matrix.cpp
index f2daff0..c7e53c7 100644
--- a/LinearAlgebra/Matrix.cpp
+++ b/LinearAlgebra/Matrix.cpp
@@ -139,7 +139,7 @@ Matrix2 Matrix2::Identity(int size) {
 }
 
 Matrix2 Matrix2::Diagonal(float f, int size) {
-  Matrix2 r = Matrix2(size, size);
+  Matrix2 r = Matrix2::Zero(size, size);
   float* data = r.data;
   int valueIx = 0;
   for (int ix = 0; ix < size; ix++) {
@@ -206,16 +206,16 @@ Matrix2 LinearAlgebra::Matrix2::operator*(const Matrix2& B) const {
     int BColOffset = i * BCols;  // BColOffset is constant for each row of B
     for (int j = 0; j < BCols; ++j) {
       float sum = 0;
-      // std::cout << " 0";
+      std::cout << " 0";
       int BIndex = j;
       for (int k = 0; k < ACols; ++k) {
-        // std::cout << " + " << this->data[ARowOffset + k] << " * "
-        //           << B.data[BIndex];
+        std::cout << " + " << this->data[ARowOffset + k] << " * "
+                  << B.data[BIndex];
         sum += this->data[ARowOffset + k] * B.data[BIndex];
         BIndex += BCols;
       }
       r.data[BColOffset + j] = sum;
-      // std::cout << " = " << sum << " ix: " << BColOffset + j << "\n";
+      std::cout << " = " << sum << " ix: " << BColOffset + j << "\n";
     }
   }
   return r;