Вычислить среднее арифметическое элементов строк матрицы

#include <math.h>
#include <iostream>
#include <iomanip>
#include <time.h>
#include <Windows.h>
 
using namespace std;
 
void crtmtx(int n1, int m1, double mt1[8][4])   //Функція створення матриці за формулою
{
    srand(time(NULL));
    int i1, j1; //Лічильники
    cout << "Матриця: " << endl;
    for (i1 = 0; i1 < n1; i1++)
    {
        for (j1 = 0; j1 < m1; j1++)
        {
            mt1[i1][j1] = (-1)*pow((2+i1/3+j1),i1)-exp(cos(j1))-3*3.14;
            cout << setw(10) << setprecision(5) << mt1[i1][j1];
        }
        cout << endl;
    }
}
void vctrmx(int n2, int m2,double mt2[8][4])   //Среднее арифметическое элементов строк матрицы
{
    int i2,j2;
    double sum_str,sr_arif;
 
    for (i2=0;i2<=n2;i2++)
{
        sum_str=0;
        sr_arif=0;
 
    for (j2=0;j2<=m2;j2++)
{
    sum_str=sum_str+mt2[i2][j2];
 
}
    cout<<"\n"<<sum_str;
    sr_arif=sum_str/n2;
    cout<<" ravno "<<sr_arif<<" \n ";
}
    cout << endl;
}
double fndG(int n3, double vc3[])   //Функція знаходження G за елементами масиву
{
    double g3;          //Вихідні данні
    double sum = 0;     //Сума
    double p;           //Добуток
    int k = 1, n = 1;   //Лічильники
    while (k <= 8)
    {
 
        p = 1;
        while (n <= k)
        {
            p *= vc3[n];
            n++;
        }
        sum += vc3[k] * vc3[k] - p;
        k++;
    }
    g3 = sqrt(fabs(sum));
    return g3;
}
int main()
{
    SetConsoleCP(1251);
    SetConsoleOutputCP(1251);
    //srand(time(NULL));
    int n = 8, m = 4;   //Кількість елементів
    double mt[8][4];    //Матриця
    double vc[8];       //Вектор
    double g;           //Вихідна змінна
    crtmtx(n, m, mt);
    vctrmx(n, m, mt);
    g = fndG(n, vc);
    cout << "G= " << g << endl;
}

#define _USE_MATH_DEFINES
#include <iostream>
#include <vector>
#include <cmath>
using namespace std;
 
double G (const vector<double>& x)
{
    double sum=0.0, mult;
    for (int k=1; k<=8; k++)
    {
       mult=1.0;
       for (int n=1; n<=k; n++)
       mult*=x[n-1];       
       sum+=x[k-1]*x[k-1]-mult;
    }   
    return sqrt(fabs(sum));
}
 
int main()
{
    const int N=8, M=4;
    double a[N][M];
    vector<double> x(N);
    
    cout << "Matrix:\n";
    for (int i = 0; i < N; i++)
    {
       x[i]=0.0;
       for (int j = 0; j < M; j++)
       {
       a[i][j]=pow(-(double)(2+i)/(3+j),i)-exp(cos(j))-3.*M_PI;
       x[i]+=a[i][j];
       cout << a[i][j] << "  ";
       }
       cout << "\n";
       x[i]/=M;
    }
    
    cout << "Vector:\n";
    for (int i = 0; i < N; i++) 
    cout << x[i] << "  ";
    cout << "\n";   
    
    cout << "G=" << G(x) << "\n";
 
system("pause");
return 0;
}