Найти ошибку в методе Симпсона

using System;
 
namespace Simpson
{
   public class Program
    {
      public  static void Main(string[] args)
        {
            double x, a, b, h, s,s2,e,n;
            m1:
            a = ((2 - Math.Sqrt(2)) / 2);
        b = ((2 + Math.Sqrt(2)) / 2);
            Console.WriteLine("Введiть парне ціле значення N(кількість частин):");
            n = Convert.ToDouble(Console.ReadLine());
        Console.WriteLine("Введiть e");
            e = Convert.ToDouble(Console.ReadLine());
      if(n%2==1 || n<1)
      {Console.WriteLine("Некоректно введені дані"); goto m1;}
          else{
        double [] r = new double [40];
      for(int i=1; r[i]-r[i-1]<=e; i++)
      {
      h = (b- a) / n;
            s = 0;
            x = a+h;
                            
            while (x < a)
            {
                s += 4 *(1-2*x+x*x)*x ; 
                x += h;
                s += 2 * (1-2*x+x*x)*x;
                x += h;
            }
       
            s2 = 0;
        
            while (x < b)
            {
                s2 += 4 *(2*x-x*x)*x;
                x += h;
                s2 += 2 *(2*x-x*x)*x;
                x += h;
            }
      s = h / 3 * (s + ((1-2*x+x*x) - (1-2*x+x*x)));
            s2 = h / 3 * (s2 + ((2*x-x*x) - (2*x-x*x)));
        double ss = s2-s;
      r[0] = 0;
      r[i]= ss;
     
      if(r[i]-r[i-1]<=e)
      { Console.WriteLine("Интеграл = " + r[i]);
        Console.WriteLine("Для цієї точності N= {0}",n);return;
      } n*=2;
      }}}}}