Решение интеграла методом средних прямоугольников с заданной точностью

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Windows.Forms;
 
namespace my_kursach
{
    public partial class Form1 : Form
    {
 
        static double f(double x)
        {
            return 1 / (x * x + 1);
        }
        static double f2(double x)
        {
            return (x*x/(8+x*x*x));
        }
        static double f3(double x)
        {
            return Math.Cos(x);
                
        }
        //static double rect(double a, double b, double n)
        //{            
 
        //double h, i, value = 0;
 
        //h = (b - a) / n;
 
        //for (i = a + h / 2; i<b; i += h)
 
        //value += f(i)*h;
 
        //return value;
        //}
        
  
        public Form1()
        {
            InitializeComponent();
        }
 
        private void button1_Click(object sender, EventArgs e)
        {
            double a, b, n,c;
            a = Convert.ToDouble(textBox1.Text);
            b = Convert.ToDouble(textBox2.Text);         
            n= Convert.ToDouble(textBox3.Text);
            const int eps = 5;
            const int ep = 3;
 
            
            if (a > b)
            {
 
                c = a;
 
                a = b;
 
                b = c;
            }
            // Настраиваем оси графика
            chart1.ChartAreas[0].AxisX.Minimum = a;
            chart1.ChartAreas[0].AxisX.Maximum = b;            
            chart1.ChartAreas[0].AxisX.MajorGrid.Enabled = false;//убираем сетку
            chart1.ChartAreas[0].AxisY.MajorGrid.Enabled = false;
            chart1.ChartAreas[0].AxisX.Crossing = 0;//пересечение осей
            chart1.ChartAreas[0].AxisY.Crossing = 0;
            chart1.Series[0].Points.Clear();
            chart1.Series[1].Points.Clear();
 
 
        if (radioButton1.Checked == true)
        {
 
            double s = 0;
            int i;
            double RndS = Math.Round(s, eps);
            double h = (b - a) / n;
 
 
            for (i = 0; i < n; i++)
            {
 
                s += f2(a + h * (i + 0.5));
                chart1.Series[0].Points.AddXY((Math.Round(a + h * (i + 0.5), ep)), f2(a + h * (i + 0.5)));
                chart1.Series[1].Points.AddXY((Math.Round(a + h * (i), ep)), f2(a + h * (i + 0.5)));
            }
 
            s *= h;
            chart1.Series[0].Points.AddXY((Math.Round(a + h * (i + 0.5), ep)), f2(a + h * (i + 0.5)));
            chart1.Series[1].Points.AddXY((Math.Round(a + h * (i), ep)), f2(a + h * (i + 0.5)));
            RndS = Math.Round(s, eps);
            textBox4.Text = RndS.ToString();
 
        }
        
        }
 
        private void button2_Click(object sender, EventArgs e)
        {
           this.Close();
        }
       
    }
}