Код нахождения неопределенного интеграла

double IntervalWidth = 0.0023;
            // интеграл от
            double a = 0;
            // вывод на экран
           
            // интеграл до
            double b = 2.0;
            // вывод на экран
           
            // количество интервалов   
            int numbIntervals = (int)((b - a) / IntervalWidth);
            // высота интервала
            double IntervalHeight = 0.0;
            // координата текущей точки на оси x
            double x = 0.0;
            // текущий интервал
            int Interval = 0;
            // локальная сумма (на одном компьютере)
            double summaLocal = 0.0;
            // цикл от первого до последнего интервала (от точки а на оси х до точки // b) с шагом 1
            for (Interval = 0; Interval <= numbIntervals; Interval++)
            {
                x = a + (IntervalWidth * ((double)Interval - 0.5));
                IntervalHeight += (3.17 * x + 4.2) /Math.Pow(Math.Cos(x),2);
            }
            summaLocal = IntervalWidth * IntervalHeight;
 
            Console.Write("Integral is approximately  "+summaLocal);
            Console.ReadKey();