Class и operator[] для N мерных массивов (требуется совет от опытных)

 /* В VS2017 работает на Ура! и CodeBlocks компилирует без "варнингов" */
template<typename T> class Matrix {
private:
    T **matrix;
    size_t N, M;
public:
    Matrix(size_t _N, size_t _M) : N(_N), M(_M) {
        this->matrix = new T*[this->N];
        for (size_t i(0); i < this->N; i++)
            this->matrix[i] = new T[this->M];
    };
    T* operator[](size_t index) { return this->matrix[index]; }
    .
    .
    .
    ~Matrix(void) {
        for (size_t i(0); i < this->N; i++)
            delete this->matrix[i];
        delete this->matrix;
    };
};

#include <iostream>
using namespace std;
 
template<size_t N, size_t M, typename T>
struct array_type{
    struct primarry_array{
    T A[M];
    T &operator[](int i){return A[i];}
    };
 
primarry_array B[N];
    
    primarry_array &operator[](size_t n){
        return B[n];
    }   
 
};
 
int main(int argc, char* argv[])
{
const int N=3, M=3;
array_type<N,M, int> a;
for(int i=0; i<N; ++i){
    for(int j=0; j<M; ++j){
a[i][j]=i*j+j;
    }
}
for(int i=0; i<N; ++i){
    for(int j=0; j<M; ++j){
cout<<a[i][j]<<' ';
    }
    cout<<endl;
}
 
cout<<endl;
system("pause");
return 0;
}

T end_of_secondary_array;
T* end_of_primary_array=&end_of_secondary_array;