#include <iostream.h>


class Matrix3x3
{
public:
   double A[3][3];


   Matrix3x3 addMatrix(Matrix3x3 *m2)
   {
      Matrix3x3 out;
      int i, j;
   
      for (i = 0; i < 3; i = i+1)
         for (j = 0; j < 3; j = j+1)
            out.A[i][j] = this->A[i][j] + m2->A[i][j];
   
      return(out);
   }
};

void print(Matrix3x3 *m)
{
   int i, j;

   for (i = 0; i < 3; i = i+1)
   {
      for (j = 0; j < 3; j = j+1)
         printf("%5.2lf ", m->A[i][j]);
      cout << endl;
   }
}

int main(int argc, char *argv[])
{
   Matrix3x3 A, B, C;

   A.A[0][0] = 1.0; A.A[0][1] = 0.0; A.A[0][2] = 0.0;
   A.A[1][0] = 1.0; A.A[1][1] = 1.0; A.A[1][2] = 0.0;
   A.A[2][0] = 1.0; A.A[2][1] = 1.0; A.A[2][2] = 1.0;

   cout << "Initial matrix A:" << endl << endl;
   print(&A);

   B.A[0][0] = 1.0; B.A[0][1] = 2.0; B.A[0][2] = 3.0;
   B.A[1][0] = 4.0; B.A[1][1] = 5.0; B.A[1][2] = 6.0;
   B.A[2][0] = 7.0; B.A[2][1] = 8.0; B.A[2][2] = 9.0;

   cout << "Initial matrix B:" << endl << endl;
   print(&B);

   C = B.addMatrix(&A);     // <--- Member function

   cout << endl << "After C = B.addMatrix(&A), matrix C:" << endl << endl;
   print(&C);
}