public double Integral(double[] x, double intPointOne, double intPointTwo)
{
    double integral = 0;
    double i = intPointOne;
    do
    {
        integral += Function(x[i])*.001;
        i = i + .001;
    }
    while (i <= intPointTwo);
    return integral;
}

Здесь функция, которую я должен интегрировать из x1-x2, просто используя суммирование частей. Как я могу сделать этот цикл более эффективным (используя меньше циклов), но более точным?

Where Function changes every iteration, but it should be irrelevant as it's order of magnitude (or boundary) should stay relatively the same...