Рекурсивный поиск элементов матрицы

#include <iostream>
 
using namespace std;
 
 
void enterLineMatrix(double** matrix, int i, int M)
{
    static int j = 0;
    if (M != 0)
    {
        std::cout << "[" << i << "]" << "[" << j << "] -> ";
        std::cin >> matrix[i][j];
        if (j != M - 1)
        {
            j++;
            enterLineMatrix(matrix, i, M);
        }
        else
        {
            j = 0;
            return;
        }
    }
}
 
void enterMatrix(double** matrix, int N, int M)
{
    static int i = 0;
    if (N != 0)
    {
        enterLineMatrix(matrix, i, M);
        if (i != N - 1)
        {
            i++;
            enterMatrix(matrix, N, M);
        }
        else
        {
            i = 0;
            return;
        }
    }
}
 
void OutLineMatrix(double** matrix, int i, int M)
{
    static int j = 0;
    if (M != 0)
    {
        std::cout << matrix[i][j] << "\t";
        if (j != M - 1)
        {
            j++;
            OutLineMatrix(matrix, i, M);
        }
        else
        {
            j = 0;
            return;
        }
    }
}
 
void OutMatrix(double** matrix, int N, int M)
{
    static int i = 0;
    if (N != 0)
    {
        OutLineMatrix(matrix, i, M);
        std::cout << std::endl;
        std::cout << std::endl;
        if (i != N - 1)
        {
            i++;
            OutMatrix(matrix, N, M);
        }
        else
        {
            i = 0;
            return;
        }
    }
}
 
double foo_1(double** matrix, int N, int M)// 4 string
{
    double sum = 0;
    if(N >= 4)
    {
        for (size_t i = 0; i < M; i++)
        {
            sum += matrix[3][i];
        }
    }
    return sum;
}
 
double foo_2(double** matrix, int N, int M) // average
{
    double average = 0;
    if (N < 2 && N != 0)
    {
        for (size_t i = 0; i < M; i++)
        {
            average += matrix[0][i];
        }
        average = average / M;
    }
    else
    {
        for (size_t i = 0; i < M; i++)
        {
            average += matrix[N - 2][i];
        }
        average = average / M;
    }
    return average;
}
 
int main() 
{
    int N;
    int M;
    std::cout << "Enter N: ";
    std::cin >> N;
    std::cout << "Enter M: ";
    std::cin >> M;
 
 
    double** matrix = new double* [N];
    for (size_t i = 0; i < N; i++)
    {
        matrix[i] = new double[M];
    }
    std::cout << "Enter matrix: " << std::endl;
    enterMatrix(matrix, N, M);
    std::cout << "Your matrix: " << std::endl;
    OutMatrix(matrix, N, M);
    std::cout << std::endl;
    std::cout << "4-th string sum: " << foo_1(matrix, N, M) << std::endl;
    std::cout << "Pre last string average: " << foo_2(matrix, N, M) << std::endl;
    std::cout << "Result: " << foo_1(matrix, N, M) * foo_2(matrix, N, M) << std::endl;
    return 0;
}