Вычисление определенного интеграла методом Монте-Карло

    public partial class Lab3 : Form
    {
      static Random rand = new Random();
       static int numb = rand.Next(5, 30);
        public Lab3()
        {
            InitializeComponent();
            textBox_b.Text = Convert.ToString(Math.PI / 2);
            textBox_n.Text = Convert.ToString(numb);
        }
 
        public double f(double x)
        {
            return Math.Cos(x) / (1 + x * x);    
        }
 
        private void button1_Click(object sender, EventArgs e)
        {
            #region Input
            int n, i;
            double a, b, max, min;
 
            a = double.Parse(textBox_a.Text);
            b = double.Parse(textBox_b.Text);
            n = int.Parse(textBox_n.Text);
 
            #endregion
            #region condition
            if (f(a) >= f(b))
            {
                max = f(a); // max значение функции
                min = f(b); // min значение функции
            }
            else
            {
                max = f(b);
                min = f(a);
            }
            #endregion
            #region ...
            double[] x = new double[n];
            double[] y = new double[n];
            int count = 0;
            Random rand = new Random();
            for (i = 0; i < n; i++)
            {
                x[i] = a + (b - a) *rand.Next(5, 30); ; // random N
             chart1.Series["Series1"].Points.DataBindXY(x, y);  // Точки на граффике
 
               y[i] = min + (max - min) *rand.Next(5, 30);
                if (y[i] > f(x[i]))
                {
                    count++;
                }
            }
            double integral = (count * (b - a) * (max - min)) / n;
            richTextBox1.Text += string.Format("\nmin = {0}", min);
            richTextBox1.Text += string.Format("\nmax = {0}", max);
            richTextBox1.Text += string.Format("\nN = {0}", count);
            richTextBox1.Text += string.Format("\nIntegral is = {0}", integral);
            int A = Convert.ToInt32(a);
            double func;
            for(i = A; i<b; i++)
            {
                func = Math.Cos(i) / (1 + i * i);
            chart1.Series["Series2"].Points.AddXY(i, func);
        }
            #endregion
        }
    }

            Console.WriteLine("Integrating function f(x) = x * x from 0 to 1\n");
            Random rand = new Random();
            int nPts = 1000000;
            double s0 = 0;
            double sigma = 0;
            double x = 0.0;
            double fcn;
            for (int i = 0; i < nPts; ++i)
            {
                x = rand.NextDouble();
                fcn = x * x;
                s0 += fcn;
                sigma += fcn * fcn;
            }
            s0 /= nPts;
            sigma /= nPts;
            sigma = Math.Sqrt(Math.Abs(sigma - s0 * s0) / nPts);
            Console.WriteLine("Analytical result           = 0.3333333....333");
            Console.WriteLine("Monte Carlo Integral result = " + s0 + " +- " + sigma);

sigma = Math.Sqrt(Math.Abs(sigma - s0 * s0) / nPts);