Задать элемент массива функций как параметр функции

#include <iostream>
#include <stdio.h>
using namespace std;
int const n = 1000; // Число узлов.
double const PI = 3.14159265358979323846264;
double Fi1(double x)
{
    return 1;
}
 
double Fi2(double x)
{
    return 2*x;
}
 
double Fi3(double x)
{
    return 4*x*x-1;
}
 
double K(double x, double t)
{
    return t*(t-x);
}
 
//double InFunction(double x) // Подинтегральная функция.
//{
//  return (x - sin(x))*sin(x);
//}
//double InFunction2(double x) // Подинтегральная функция.
//{
//  return sin(x)*sin(x);
//}
//double Proverka(double x) // Подинтегральная функция проверки.
//{
//  return x*sin(x);
//}
double CalcIntegral(double a, double b, double h, double(*Function)(double x))// Подсчет интеграла по методу средних прямоугольников.
{
    double result = 0;
    for (int i = 0; i < n; i++)
        result += Function(a + h * (i + 0.5));
    result *= h;
    return result;
}
 
 
double CalcIntegralDv(double a, double b, double h, double(*K)(double x, double t), double x[], double t[], double(*FiJ)(double x), double(*FiI)(double x))// Подсчет интеграла по методу средних прямоугольников.
{
    double **q = new double*[n]; // две строки в массиве
    for (int i = 0; i < 2; i++)
        q[i] = new double[n]; // и пять столбцов
    double result = 0;
    
            for (int i = 0; i < n; i++)
            for (int j = 0; j < n; j++)
            {
                if (i == 0 && j == 0 && i == n - 1 && j == n - 1)
                    q[i][j] = 1 / 4;
                if (i == 0 && j != 0 && i == n - 1 && j != n - 1 || j == 0 && i != 0 && j == n - 1 && i != n - 1)
                    q[i][j] = 1 / 2;
                if (i != 0 && j != 0 && i != n - 1 && j != n - 1)
                    q[i][j] = 1;
                result += q[i][j] * K(x[i], t[j])*FiJ(x[i])*FiI(x[j]);
            }
            result *= h*h;
}
 
 
double Y(double x, double ba) // Нахождение функции y(x).
{
    double result = ba*sin(x) + x - sin(x);
    return result;
}
double Pogreshnost(double(*Y)(double x, double ba), double(*CalcIntegral)(double a, double b, double h, double(*Function)(double x)), double ba, double a, double b, double h, double x)// Определение погрешности счета.
{
    return result;
}
int main()
{
    double(*FI[3])(double) = { Fi1, Fi2, Fi3 };// объявления массива функций.
    double  a = -1, b = 1, h = (b - a) / n, *x = new double[n], *t = new double[n], *d = new double[3], *c = new double[9];
    //ba = CalcIntegral(a, b, h, &InFunction) / (1 - CalcIntegral(a, b, h, &InFunction2));// Нахождение коэффициента А.
    x[0] = a;
    t[0] = a;
    for (int i = 0; i < n; i++)
        x[i + 1] = x[i] + h;// Разбиение сетки по оси ОХ.
    for (int i = 0; i < n; i++)
        t[i + 1] = t[i] + h;// Разбиение сетки по оси Оt.
    for (int i = 0; i < 3; i++)
        d[i] = CalcIntegral(a, b, h, &Fi1);
    for (int k = 0; k < n; k++)
        for (int m = 0; m < n; m++)
            for (int i = 0; i < 9; i++)
                c[i] = CalcIntegral(a, b, h,(*FI[k])(x[i]), (*FI[m])(x[i])) - CalcIntegralDv(a, b, h, &K, x, t, (*FI[k])(x[i]), (*FI[m])(x[i]));// тут пытаюсь задать элемент массива как параметр функции CalcIntegralDv
    //for (int i = 0; i < n; i++)
    //  cout << Pogreshnost(&Y, &CalcIntegral, ba, a, b, h, x[i]) << endl;
    //cout << -Y(x[i], ba) + sin(x[i])*CalcIntegral(a, b,h, &Proverka ) + x[i] - sin(x[i]) << endl; // Подсчет погрешности внутри функции main().
    return 0;
}

CalcIntegral(a, b, h, (*FI[k])(x[i]), (*FI[m])(x[i]))

(*FI[k])(x[i])

FI[k](x[i])