Разложение функции в ряд Фурье

void calculation()
        {
            //Разложение в ряд Фурье
            int n = inArray.Length; //всего точек
            int z = tb_K.Value; //расчет до k-й гармоники
 
            //нахождение A0
            double a0 = 0;
            for (int i = 0; i < n; i++)
                a0 += inArray[i].F;
            a0 = a0 / n;
 
            //находим коэффициенты A и уголы psi
            double[] a = new double[z]; // A(k)
            double[] p = new double[z]; // Psi (k)
            for (int k = 0; k < z; k++)
            {
                double b = 0;
                double c = 0;
                for (int i = 0; i < n; i++)
                {
                    double arg = (k + 1) * inArray[i].x;
                    arg *= Math.PI / 180;
                    b += inArray[i].F * Math.Sin(arg);
                    c += inArray[i].F * Math.Cos(arg);
                }
                b = 2 * b / n;
                c = 2 * c / n;
                a[k] = Math.Sqrt(b * b + c * c);
                p[k] = Math.Atan(c / b);
            }
 
            //ВЫВОД
            // ...
        }