Esempio n. 1
0
        static double[] Da_thuc_Newton(double[,] Tysaiphan, double[] X)
        {
            int n = X.Length;

            double[] Pn = new double[n];
            int      k  = 0;

            for (int i = n - 1; i >= 0; i--)
            {
                if (i == 0)
                {
                    Pn[n - 1] += Tysaiphan[0, 0];
                }

                if (i == 1)
                {
                    Pn[n - 1] -= Tysaiphan[0, 1] * X[0];
                    Pn[n - 2] += Tysaiphan[0, 1];
                }

                if (i > 1)
                {
                    double[] Xi = new double[i];

                    for (int j = 0; j < i; j++)
                    {
                        Xi[j] = X[j];
                    }

                    double[] Horner = HornerPn.Nhandathuc(Xi);

                    for (int j = 0; j < Horner.Length; j++)
                    {
                        Horner[j] *= Tysaiphan[0, i];
                        Pn[j + k] += Horner[j];
                    }
                    k++;
                }
            }

            return(Pn);
        }