Деление многочлена на многочлен

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
 
namespace Алгоритм_Штурмана
{
    class Program
    {
        static void Main(string[] args)
        {
            double k=0,q;
            double[] e = new double[2];
            bool[] j = { true, true, true, true, true, true };
            double[] c = { 0, 0, 0, 0, 0, 0 };
            Console.Write("Aralyk\ne1=");//границы
            e[0] = double.Parse(Console.ReadLine());
            Console.Write("e2=");
            e[1] = double.Parse(Console.ReadLine());
            Console.Write("Write the maximum degree - ");
            int f = int.Parse(Console.ReadLine());
            double[,] d = new double[f + 1,f];
            for (int i = 0; i <= f; i++) {
                Console.Write((i + 1) + "th multiplier - ");
                d[i, 0] = int.Parse(Console.ReadLine());
                d[i, 1] = d[i, 0] * (f - i);
                c[0] = c[0] + d[i, 0] * Math.Pow(e[0], (f - i));
                c[1] = c[1] + d[i, 0] * Math.Pow(e[1], (f - i));
                c[4] = c[2] + d[i, 1] * Math.Pow(e[0], (f - i - 1));
                c[5] = c[3] + d[i, 1] * Math.Pow(e[1], (f - i - 1));
                c[2] = c[2] + d[i, 1] * Math.Pow(e[0], (f - i - 1));
                c[3] = c[3] + d[i, 1] * Math.Pow(e[1], (f - i - 1));
            }
            for (int i = 0; i < 6; i++)
                if (c[i] < 0)
                    j[i] = false;
            Console.Write("f0=" + d[0, 0] + "*x^" + f);//функция 1ая
            for (int i = 1; i < f; i++)
                Console.Write("+" + d[i, 0] + "*x^" + (f - i));
            Console.Write("+" + d[f, 0] + "\nf1=" + d[0, 1] + "*x^" + (f - 1));//ее производная, 2 функция
            for (int i = 1; i < f-1; i++)
                Console.Write("+" + d[i, 1] + "*x^" + (f - i));
            Console.Write("+" + d[f - 1, 1]);
            for (int i = 0; i < f; i++) {//деление многочлена на многочлена и она же 3 функция
                q = d[0, i] / d[0, i + 1];
                for(int l = 0; l <= f-i; l++) {
                    d[l, i + 2] = d[l, i] - d[l, i + 1] * q;
                    k = d[l, i] - d[l, i + 1] * q;
                }
            }
            Console.Write("f2=" + k);
        }
    }
}

for (int i = 0; i < f; i++) {//деление многочлена на многочлена и она же 3 функция