Значение обращается в NaN

            double[,] A = new double[3, 3];
            double[] B = new double[3];
            double[] C;
            for (int i = 2; i < ApprX.Count; i += 2)
            {
                
                for (int j = 0; j < A.GetLength(0); j++)
                {
                    for (int v = 0; v < A.GetLength(1); v++)
                    {
                        A[j, v] = Math.Pow(ApprX[i - 2 + j], 2 - v);
                    }
                    B[j] = ApprY[i - 2 + j];
                }
                C = Matrix.GaussJordan(A, B);
                double d=C[1];
 
                for (int j = Array.IndexOf(X, ApprX[i - 2]); j < Array.IndexOf(X, ApprX[i]); j++)
                {
                    YResult.Add(C[0] * Math.Pow(X[j], 2) + C[1] * X[j] + C[2]);
                }
            }
            YResult.Add(ApprY.Last());
            return YResult;

public static double[] GaussJordan(double[,] system, double[] freeMembers)
        {
            var result = new double[freeMembers.Length];
            double Multi1, Multi2;
            for (int k = 0; k < system.GetLength(0); k++)
            {
                for (int j = k + 1; j < system.GetLength(0); j++)
                {
                    Multi1 = system[j, k] / system[k, k];
                    for (int i = k; i < system.GetLength(0); i++)
                        system[j, i] = system[j, i] - Multi1 * system[k, i];
 
                    freeMembers[j] = freeMembers[j] - Multi1 * freeMembers[k];
                }
            }
            for (int k = system.GetLength(0) - 1; k >= 0; k--)
            {
                Multi1 = 0;
                for (int j = k; j < system.GetLength(0); j++)
                {
                    Multi2 = system[k, j] * result[j];
                    Multi1 += Multi2;
                }
                result[k] = (freeMembers[k] - Multi1) / system[k, k];
            }
            return result;
        }