static public RegressionData polinomialAutoRegression(double[,] X, double[] Y)
       {
           RegressionData data = new RegressionData(X, Y);
           double[,] XData = new double[(X.Length / Y.Length)*3-2, Y.Length];
           double[] YData = new double[Y.Length];
           int k = 1;
           for (int j = 0; j < Y.Length; j++)
           {
               XData[0, j] = X[0, j];
               YData[j] = Y[j];

           }
           for (int i = 1; i < X.Length / Y.Length; i++)
           {
               for (int j = 0; j < Y.Length; j++)
               {
                   XData[k, j] = X[i, j];
                   XData[k + 1, j] = Math.Pow(X[i, j], 2);
                   XData[k + 2, j] = Math.Pow(X[i, j], 3);
               }
               k += 3;
           }

               data = baseRegression(XData, YData);
               int s = 1;
           for (int i = 0; i < data.Y.Length; i++)
           {
               data.Y[i] = data.Koef[0];

           }

           for (int i = 1; i < X.Length / Y.Length; i++)
           {
               for (int j = 0; j < Y.Length ; j++)
               {
                   data.Y[j] += data.Koef[s] * (X[i, j]) + data.Koef[s + 1] * Math.Pow((X[i, j]), 2) + data.Koef[s + 2] * Math.Pow((X[i, j]), 3);
               }
               s += 3;
           }
               return data;

       }