Многочлен Лагранжа

 public double lagr(double x)
        {
            double L = 0;
            double[] l = new double[Convert.ToInt32(textBox1.Text) + 1];
            for (int i =0; i<= Convert.ToInt32(textBox1.Text); i++)
            {
                l[i] = 1;
                for (int j = 0; j < Convert.ToInt32(textBox1.Text) + 1; j++)
                {
                    if(j != i)
                    {
                        l[i] = l[i] * ((x - X[j]) / (X[i] - X[j]));
                    }
                }
                L = L + func(X[i]) * l[i];
            }
            return L;
        }

     public static double LagrangeInterpolation(double[] x, double[] y, double xval)
        {
            double yval = 0.0;
            double Products = y[0];
            for (int i = 0; i < x.Length; i++)
            {
                Products = y[i];
                for (int j = 0; j < x.Length; j++)
                {
                    if (i != j)
                    {
                        Products *= (xval - x[j]) / (x[i] - x[j]);
                    }
                }
                yval += Products;
            }
            return yval;
        }
 
        public static double[] LagrangeInterpolation(double[] x, double[] y, double[] xvals)
        {
            double[] yvals = new double[xvals.Length];
            for (int i = 0; i < xvals.Length; i++)
                yvals[i] = LagrangeInterpolation(x, y, xvals[i]);
            return yvals;
        }