求两个矩阵之和

#include<iostream>
using namespace std;

class Matrix
{
public:
  Matrix();
  friend Matrix operator+(Matrix &a, Matrix &b);
  friend ostream& operator<<(ostream &, Matrix &);
  friend istream& operator>>(istream &, Matrix &);
private:
  int mat[2][3];
};

Matrix::Matrix()
{
  for(int i=0;i<2;i++)
    for (int j = 0; j < 3; j++)
    {
      mat[i][j] = 0;
    }
}

Matrix operator+(Matrix &a, Matrix &b)
{
  Matrix c;
  for (int i = 0; i < 2; i++)
    for (int j = 0; j < 3; j++)
    {
      c.mat[i][j] = a.mat[i][j]+b.mat[i][j];
    }
  return c;
}

ostream& operator<<(ostream &out, Matrix &a)
{
  for (int i = 0; i < 2; i++)
    for (int j = 0; j < 3; j++)
    {
      out<<a.mat[i][j]<<" ";
    }
  out << endl;
  return out;
}
istream& operator>>(istream &in, Matrix &a)
{
  for (int i = 0; i < 2; i++)
    for (int j = 0; j < 3; j++)
    {
      in>>a.mat[i][j];
    }
  return in;
}
int main()
{
  Matrix a, b, c;
  cin >> a;
  cin >> b;
  cout << "a:" << endl;
  cout<<a<<endl;
  cout << "b: " << endl;
  cout<<b<<endl;
  c = a + b;
  cout << "c:" << endl;
  cout<<c;
  system("pause");
  return 0;
}