[C/C++] Matrix addition mit Funktionen

#include <iostream>
using std:: cout;

int i, j;

int main ()
{

float Matrix1 [3][3]=
{
 {6.87f ,4.56f ,28.0f},
 {87.3f ,7.61f ,34.1f},
 {8.7f ,1.92f ,7.83f}
};
float Matrix2 [3][3]=
{
 {73.0f ,23.5f ,12.0f},
 {1.93f ,9.0f ,93.2f},
 {64.74f ,82.1f ,562.0f}

};

float Matrix3[3][3]={0};

void addition (float Matrix3,float Matrix2,float Matrix1);

for (i=0; i<3; i++) //Addition Matrix1 und 2
{
 for(j=0; j<3; j++)
 {
 Matrix3[i][j]=Matrix1[i][j]+Matrix2[i][j];
 }
}
for (i=0; i<3; i++) //Ausgabe Matrix1
{
 for(j=0; j<3; j++)
 {
 if (i==1 && j==0 || i==2 && j==0)
 {
 cout <<"\n" <<Matrix1[i][j]<< " ";
 }
 else
 {
 cout <<Matrix1[i][j]<<" ";
 }
 }
}

cout <<"\n" <<"+"<<"\n";

for (i=0; i<3; i++) //Ausgabe Matrix2
{
 for(j=0; j<3; j++)
 {
 if (i==1 && j==0 || i==2 && j==0)
 {
 cout <<"\n" <<Matrix2[i][j]<< " ";
 }
 else
 {
 cout <<Matrix2[i][j]<<" ";
 }
 }
}

cout <<"\n" <<"="<<"\n";

for (i=0; i<3; i++) //Ausgabe Matrix3
{
 for(j=0; j<3; j++)
 {
 if (i==1 && j==0 || i==2 && j==0)
 {
 cout <<"\n" <<Matrix3[i][j]<< " ";
 }
 else
 {
 cout <<Matrix3[i][j]<<" ";
 }
 }
}
 system ("pause");
 return 0;
}