Определить столбец матрицы, которому принадлежит максимальный элемент главной диагонали

#include<iostream>
#include<vector>
#include<cmath>
 
 
using namespace std;
 
void Formula(  double**arr, int n, int m)
{
     double PI = 3.14;
 
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < m; j++)
        {
            arr[i][j] = pow(2.5, i + j)*pow(log(pow(8,((i+2)/2))),3.0)+pow((-3.2),i+j+4)+pow(PI,5*i+3*j);
        }
    }
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < m; j++)
        {
            cout << arr[i][j] << " ";
        
        }
        cout << endl;
    }
 
 
}
 
void Vector(double** arr, int n, int m)
{
 
    for (int i = 0; i < n; ++i)
    {
        int max = arr[i][i];
        for (int j = 0; j < n; ++j) {
            if (arr[j][j] > max)
            {
                max = arr[j][j];
 
            }
        }
        cout << max<< " "<<i<<" "<< endl;
        break;
    }
 
 
}  
int main(int argc, char* argv[])
{
    int n, m;
 
    cout << "N: "; cin >> n;
    cout << "M: "; cin >> m;
 
     double** arr = new   double* [n];
    for (int i = 0; i < n; i++)
    {
        arr[i] = new   double[m];
    }
 
    Formula(arr, n, m);
    Vector(arr, n, m);
 
 
    system("pause");
    return EXIT_SUCCESS;
}

#define _USE_MATH_DEFINES
#include <iostream>
#include <cstdlib>
#include <cmath>
using namespace std;
 
double* CreatVector (double**a, int n, int jmax) 
{
    double* b = new double[n];
    for (int i = 0; i < n; i++)    
          b[i]=a[i][jmax];    
    return b;
}
 
double MultNegative (double*b, int n)
{
    double mult=1.;
    for (int i = 0; i < n; i++)    
        if (b[i]<1e-9) mult*=b[i];        
    return mult;
}    
 
int main()
{
    int n, m, jmax=0;
    
    cout << "n="; cin >> n;
    cout << "m="; cin >> m;
 
      double** a = new double* [n];
    for (int i = 0; i < n; i++)    
        a[i] = new double[m];
        
    double lg;
    cout << "Matrix:\n";
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < m; j++) 
        {
            lg=log(pow(8,(double)(i+3)/(j+1)));        
            a[i][j]=pow(2.5,i+j+2)*lg*lg*lg+pow(-3.2,i+j+6)*pow(M_PI,5*(i+1)+3*(j+1)); 
            cout << a[i][j] << "  ";
        }         
    cout << "\n";
    }
    
    for (int j = 1; j < n; j++)        
        if (a[j][j]>a[jmax][jmax]) jmax=j;    
    
    double* b = CreatVector (a,n,jmax);    
    
    cout << "Vector:\n";
    for (int i = 0; i < n; i++)
        cout << b[i] << "\n"; 
        
    cout << "mult negative = " << MultNegative(b,n) << "\n"; 
 
    for (int i = 0; i < n; i++)
    delete[]a[i];
    delete[]a;
    delete[]b;
system("pause");
return 0;
}