Итерации - реализация интерполяции Лагранжа

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;
using System.IO;
 
namespace interpol2
{
    public partial class Form1 : Form
    {
        public Form1()
        {
            InitializeComponent();
            openFileDialog1.DefaultExt = "txt";
            openFileDialog1.Filter = ".txt|*.txt";
            openFileDialog1.Title = "Open file";
            openFileDialog1.Multiselect = false;
            openFileDialog1.FileName = string.Empty;
        }
 
        string FNwithExten = string.Empty;
        string newFNwithExten = string.Empty;
        string FN = string.Empty;
        string Dir = string.Empty;
        int QuantBeat = 0;
        int nom = 0;
        int newQB = 0;
 
        /////////////////////////////////////////
        //обработчик для дроби в правой части
        double drob(double xp, int step)
        {
            int[] time = new int[6];
            double chesl = 1;
            double znam = 1;
            for (int i = 0; i == 5; i++) { if (i == step) { continue; } chesl *= xp - i; }
            for (int i = 0; i == 5; i++) { if (i == step) { continue; } znam *= step - i; }
            return chesl / znam;
        }
        // здесь фактический обработчик массива
        void Execute(double[] arrvx, double[] arrvix, int sizevx, int kolyzlov)  
        {
            int count = 0;
            double itog = 0;
            for (int i = 0; i != sizevx-6; i++)
            {
                if (i < 2) 
                {
                   for (int Nyzla = 0; Nyzla < kolyzlov; Nyzla++)
                    {
                        arrvix[count] = itog;
                        count++;
                   }               
                }
                else{
                  for (int Nyzla = 1; Nyzla < kolyzlov; Nyzla++)
                   {
                        for (int j = 0; j == 5; j++)
                       {
                            itog += arrvx[i + j] * drob(2+(Nyzla / kolyzlov), j);
                        }
 
                    }                           
                    arrvix[count] = itog;
                    count++; 
                }
            }
        }
        ////////////////////////////////////////
        private void bOpen_Click(object sender, EventArgs e)
        {
            if (openFileDialog1.ShowDialog() == DialogResult.OK)
            {
                try
                {
                    FNwithExten = openFileDialog1.FileName;
                    FN = Path.GetFileNameWithoutExtension(FNwithExten);
                    Dir = Path.GetDirectoryName(FNwithExten);
                    StreamReader SR = new StreamReader(FNwithExten);
                    lb1.Text = FNwithExten;
                    string Textfile = SR.ReadToEnd();
                    string[] LinesText=Textfile.Split(new Char[]{'\n',' '},StringSplitOptions.RemoveEmptyEntries);
                    QuantBeat = LinesText.Length;
                    double[] arr1 = new double[QuantBeat];
                    for (int i = 0; i < QuantBeat; i++)
                    {
                        arr1[i] = Convert.ToDouble(LinesText[i]);                    
                    }
                    try
                    {
                        nom=Convert.ToInt16(tb1.Text);
                        newQB = QuantBeat * nom;
                        double[] arr2=new double[newQB];
                        try
                        {
                            Execute(arr1,arr2,QuantBeat,nom);
                            newFNwithExten = Dir +"\\"+ FN + "_inter.txt";
                            using (StreamWriter SW = new StreamWriter(newFNwithExten))
                               {
                                foreach (int i in arr2){ SW.WriteLine(i); }
                               }
                        }
                        catch{lb2.Text ="bad solution";}
                    }
                    catch { lb2.Text = "Vi vvodite ne chislo ili ne celoe chislo"; }
 
                }
                catch { lb2.Text = "not reading"; }
            }
        }
 
        private void bClose_Click(object sender, EventArgs e)
        {
            Close();
        }
 
 
 
    }
}

                    for (int i = 0; i < QuantBeat; i++)
                    {
                        arr1[i] = Convert.ToDouble(LinesText[i]);                    
                    }

                        double[] arr2=new double[newQB];

                            Execute(arr1,arr2,QuantBeat,nom);

        void Execute(ref double[] arrvx, ref double[] arrvix, int sizevx, int kolyzlov)  
        {
            int count = 0;
            double itog = 0;
            //....

                            Execute(ref arr1,ref arr2,QuantBeat,nom);

namespace PolynomialInterpolation
{
 
    public enum InterpolationType
    {
        Equidistant,
        Optimal
    }
 
    public class PolynomialInterpolation
    {
        private void NodesInitialisation(int countofnodes, InterpolationType _type)
        {
            _nodes = new double[countofnodes];
            switch (_type)
            {
                case InterpolationType.Equidistant:
                    for (int i = 0; i < countofnodes; i++)
                    {
                        _nodes[i] = a + i * (b-a) / (countofnodes - 1);
                    }
                    break;
                case InterpolationType.Optimal:
                    for (int k = 0; k < countofnodes; k++)
                    {
                        _nodes[k] = Math.Cos(Math.PI * (2 * k + 1) / (2 * countofnodes)) * (b - a) / 2 + (b + a) / 2;
                    }
                    break;
            }
        }
 
        private PolynomialInterpolation(double a, double b, int n, InterpolationType type)
        {
            this.a = a;
            this.b = b;
            NodesInitialisation(n, type);
        }
 
        public PolynomialInterpolation(double a, double b, int n, InterpolationType type, Func<double,double> f)
            : this(a, b, n, type)
        {
            _sourcefunction = f;
        }
 
        public PolynomialInterpolation(double a, double b, int n, InterpolationType type, double[] values)
            : this(a, b, n, type)
        {
            _values = values;
        }
 
        public double[] Nodes
        {
            get
            {
                return _nodes;
            }
        }
 
        public double LagrangePolynomial(double x)
        {
            double sum = 0;
            if (_sourcefunction != null)
            {
                for (int i = 0; i < _nodes.Length; i++)
                {
                    sum += _sourcefunction(_nodes[i]) * BasisPolynomial(i, x);
                }
            }
            else
            {
                for (int i = 0; i < _values.Length; i++)
                {
                    sum += _values[i] * BasisPolynomial(i, x);
                }
            }
            return sum;
        }
 
        public double BasisPolynomial(int k, double x)
        {
            double result = 1;
            for (int i = 0; i < _nodes.Length; i++)
            {
                if (i != k)
                {
                    result *= (x - _nodes[i]) / (_nodes[k] - _nodes[i]);
                }
            }
            return result;
        }
 
        private double[] _values;
        private Func<double, double> _sourcefunction;
        private double[] _nodes;
        private double a, b;
    }
}

        public PolynomialInterpolation(double a, double b, int n, InterpolationType type, Func<double, double> f)
            : this(a, b, n, type)
        {
            _sourcefunction = f;
        }
 
        public PolynomialInterpolation(double a, double b, int n, InterpolationType type, double[] values)
            : this(a, b, n, type)
        {
            _values = values;
        }
 
        public double[] Nodes
        {
            get
            {
                return _nodes;
            }
        }

for (int i = 0; i < _values.Length; i++)                 {                     sum += _values[i] * BasisPolynomial(i, x);                 }

        double drob(double xp, int step)
        {
            double result = 1;
            for (int i = 0; i == 5; i++) { if (i != step) result *= (xp - i)/(step-i); }
            return result;
        }

for (int i = 0; i == 5; i++)

        //обработчик для дроби в правой части
        double drob(double xp, int step)
        {
            for (int i = 0; i < 6; i++) { 
                if (i != step) 
                    result *= (xp - i)/(step-i);
            }
            return result;
        }
        // здесь фактический обработчик массива
        void Execute(double[] arrvx,double[] arrvix, int sizevx, int kolyzlov)  
        {
            int count = 0;
            double itog = 1;
            for (int i = 0; i != sizevx; i++)
            {
                  for (int Nyzla = 0; Nyzla < kolyzlov; Nyzla++)
                   {
                       for (int j = 0; j < 6; j++)
                       {
                           itog += arrvx[i + j] * drob(2+(Nyzla/kolyzlov), j);
                       }
                    arrvix[count] = itog;
                    count++;
                    }                           
                   
            }
        }

itog += arrvx[i + j] * drob(2.0+Convert.ToDouble(Nyzla)/Convert.ToDouble(kolyzlov), j);