Из матриц А, В и С получить матрицу (А+В)*С

#include <iostream>
#include <cmath>
using namespace std;
 
int main()
{
const int n=2;
double A[n][n], B[n][n], C[n][n];
 
for (int i=1;i<=n;i++)
    {   
        for (int j=1;j<=n;j++)
        {
            cout << "A[" << i <<"," << j << "] = ";
            cin >> A[i][j];
        }
    }
cout<<"\n";
 
for (int i=1;i<=n;i++)
    {   
        for (int j=1;j<=n;j++)
        {
            cout << "B[" << i <<"," << j << "] = ";
            cin >> B[i][j];
        }
    }
cout<<"\n";
 
for (int i=1;i<=n;i++)
    {   
        for (int j=1;j<=n;j++)
        {
            cout << "C[" << i <<"," << j << "] = ";
            cin >> C[i][j];
        }
    }
 
cout<<"\n";
 
for (int i=1; i<=n; i++)
{
        for (int j=1; j<=n; j++)
        {
            A[i][j]+= B[i][j];
           
        }
}
 
for(int i = 1; i<=n ; i++)
    for(int j = 1; j<=n; j++)
    {
      B[i][j] = 0;
      for(int k = 1; k<=n; k++)
        B[i][j] += A[i][k] * C[k][j];
    }
 
    cout<<"(A+B)*C:\n\n";
    for (int i=1;i<=n;i++)
    {   
        for (int j=1;j<=n;j++)
        {
            cout << "(A+B)*C[" << i <<"," << j << "] = "<<B[i][j]<<endl;
        }
    }
 
system ("pause");
 
 
}