Получить вещественную квадратную матрицу порядка n

template<size_t n> void foo(const double(&arr)[n], double(&mat)[n][n])
{
    for (size_t i{}; i < n; ++i)
    {
        for (size_t j{}; j < n; ++j)
        {
            if (!i)
                mat[i][j] = 0;
            else
                if (i == 1)
                    mat[i][j] = arr[j];
                else
                    mat[i][j] = mat[i - 1][j] * arr[j];
        }
    }
}

#include <iostream>
#include <iomanip>
using namespace std;
 
int main()
{
    int n;
    cout << "n="; cin >> n;
 
    double*x = new double[n];
    
    double **A = new double*[n];
    for (int i = 0; i < n; i++)
         A[i]=new double[n];
    
    for (int i = 0; i < n; i++)
        {
        cout << "x[" << i+1 << "]="; cin >> x[i];
        }
        
    for (int j = 0; j < n; j++)
    {
        A[0][j]=1.;
        for (int i = 0; i < n-1; i++)       
        A[i+1][j]=A[i][j]*x[j];    
    }  
    
    for (int i = 0; i < n; i++)
    {        
        for (int j = 0; j < n; j++)    
        cout << setw(7) << A[i][j] << " ";      
        cout << "\n";
    } 
    
    for (int i = 0; i < n; i++)
    delete[]A[i];
    delete[]A;
    delete[]x;
system("pause");
return 0;
}