Зацикливание программы, найти ошибку в коде

#define EE 1E-80
#define MAX_ROOTS 15
 
 
typedef struct
{
    Complex root[MAX_ROOTS];
    int nor; // number of root
    int deg; // degree
} Roots;
 
// возвращает номер корня полинома для заполнения области
int process_root(Complex z, Roots *p, float EPS)
{
    int i;
    for (i = 0; i < p->nor; i++)
        if (mod(csub(z, p->root[i])) < EPS)
            return i + 1;
    p->root[p->nor] = z;
    p->nor += 1;
    return p->nor;
}
 
// инициализация корней
void Rootinit(Roots *roots)
{
    register int i;
    for (i = 0; i<MAX_ROOTS; i++)
        roots->root[i] = comp(0, 0);
    roots->nor = 0;
    roots->deg = 0;
}
 
// возвращает значение полинома в точке z
Complex f(Complex z, int deg, float poly[])
{
    register int i;
    register Complex f;
    f = comp(poly[deg], 0);
    for (i = deg - 1; i >= 0; i--)
        f = cadd(cmult(f, z), comp(poly[i], 0));
 
    return f;
}
 
// возвращает значение первой производной полинома в точке z
Complex df(Complex z, int deg, float poly[])
{
    register int i;
    register Complex df;
    df = comp(deg*poly[deg], 0);
    for (i = deg - 1; i>0; i--)
        df = cadd(cmult(df, z), comp(i*poly[i], 0));
 
    return df;
}
 
int drawFrac(float dx, float per, float poly[], int deg, int NIT, float EPS)
{
    float dy, xmin, ymin, x, y;
    int X, Y;
    int ITEST;
    int i, j, k, kit, col; // счетчики, кол-во итераций, цвет
    Complex z, w, fz, dfz;
    Roots roots;
    FILE *fp;
 
 
    Rootinit(&roots);
    roots.deg = deg;
 
 
    per = sqrt(per / 18);
    //per *= 0.01;
    // масштаб
    X = getmaxx()*per;
    Y = getmaxy()*per;
    dy = (dx*Y) / X;
    xmin = -dx;
    ymin = dy;
    for (i = 0; i <= Y; i++)
    {
        y = ymin - 2 * i*dy / Y;
        for (j = 0; j <= X; j++)
        {
            x = xmin + 2 * j*dx / X;
            z = comp(x, y);
            ITEST = 1;
            // цикл до тех пор, пока не достигнем нужного числа итераций,
            // или не найдем максимально близкий корень
            for (k = 1; k <= NIT; k++)
            {
                kit = k;
                fz = f(z, roots.deg, poly);
                dfz = df(z, roots.deg, poly);
                if (mod(dfz) <= EE) { ITEST = 0; break; }
                w = z;
                z = csub(z, cdiv(fz, dfz)); // zn+1 = zn - f(zn)/f'(z)
                if (mod(csub(z, w)) <= EPS) break;
            }
            if (kit == NIT) ITEST = 2;
            if (ITEST != 1) putpixel(j, i, BLACK);
            else
            {
                col = process_root(z, &roots, EPS);
                putpixel(j, i, col);
            }
        }
        if (kbhit() != 0) break;
    }
    while (!kbhit()){}
 
    return 0;
}
 
 
 
 
int main()
{
 
 
 
    setlocale(LC_ALL, "Russian");
    int number = 0;//Номер фрактала
    while (!false)
    {
        cout << "\n1)Фрактал \"HenonIFS\"\n"
             << "2)Дерево Пифагора\n"
             << "3)Фрактал zUnknown\n"
             << "4)Фрактал Ньютона\n";
        number = 0;
        cout << "Выберите фрактал: ";
        while (!(cin >> number) || (cin.peek() != '\n') || (number <= 0) || (number >= 5))
        {
            cin.clear();
            while (cin.get() != '\n');
            cout << "\a";
            cout << "Выберите фрактал: ";
        }
        switch (number)
        {
            case 1:
            {
                int answ = 0;
                float poly[MAX_ROOTS + 1] = { -1, 0, 0, 0, 0, 1, 0 };
                float dx = 2, per = 20;
                float EPS = 1E-10;
                int i, deg = 5, NIT = 40;
                do
                {
                    printf("INPUT TO DRAW FRACTAL\n");
                    printf("1) Polynominal degree (max is %d) : \t%d\n", MAX_ROOTS, deg);
                    printf("   Coeficients : \n");
                    for (i = deg; i >= 0; i--) printf("\t\tA%d= %f\n", i, poly[i]);
                    printf("2) Width of the area on the cartesian surface (dx) : \t%f\n", dx);
                    printf("3) Percentage of the area of the screen to be covered: \t%f\n", per);
                    printf("4) Maximum Number of Iterations : \t\t\t%d\n", NIT);
                    printf("5) Precision : \t\t\t\t\t\t%G \n", EPS);
                    printf("6) Draw Fractal\n");
                    printf("7) Exit  \n");
                    printf("\nWhat is your choice ? ");
                    cin >> answ;
                    switch (answ)
                    {
                    case(1) :
                        do
                        {
                            printf("\nEnter the degree of the complex polynominal [<=%d]: ", MAX_ROOTS);
                            scanf("%d", &deg);//от 18 до 18
                        } while (deg > 15 || deg < 0);
                        printf("\nEnter the coeficients starting with the greatest power\n");
                        for (i = deg; i >= 0; i--)
                        {
                            printf("A%d = ", i);
                            scanf("%f", &poly[i]);
                        }
                        break;
                    case(2) :
                        do
                        {
                            printf("\nEnter the width (dx>0) = ");
                            scanf("%f", &dx);
                        } while (dx <= 0);
                        break;
                    case(3) :
                        do
                        {
                            printf("\nEnter the percentage of the screen to be covered by the fractal(%)\n");
                            printf("(the smallest the fastest ->useful for a preview) : ");
                            scanf("%f", &per);
                        } while (per <= 0 || per > 100);
                        break;
                    case(4) :
                        do
                        {
                            printf("\nEnter the maximum Number of ITerations (NIT>0) : ");
                            scanf("%d", &NIT);
                        } while (NIT <= 0);
                        break;
                    case(5) :
                        do
                        {
                            printf("\nEnter precision ( >0 ) : ");
                            scanf("%e", &EPS);
                        } while (EPS < 0);
                        break;
                    case(6) :
                        initwindow(700, 700);
                        drawFrac(dx, per, poly, deg, NIT, EPS);
                        WaitEscEnter();
                        closegraph();
                        break;
                    case(7) :
                        exit(0);
                    }
                } while (1);/*end switch */
                break;
            }
            case 2:....
    
        }
    }
    return 0;
}

                    printf("INPUT TO DRAW FRACTAL\n");
                    printf("1) Polynominal degree (max is %d) : \t%d\n", MAX_ROOTS, deg);
                    printf("   Coeficients : \n");
                    for (i = deg; i >= 0; i--) printf("\t\tA%d= %f\n", i, poly[i]);
                    printf("2) Width of the area on the cartesian surface (dx) : \t%f\n", dx);
                    printf("3) Percentage of the area of the screen to be covered: \t%f\n", per);
                    printf("4) Maximum Number of Iterations : \t\t\t%d\n", NIT);
                    printf("5) Precision : \t\t\t\t\t\t%G \n", EPS);
                    printf("6) Draw Fractal\n");
                    printf("7) Exit  \n");
                    printf("\nWhat is your choice ? ");

int answ = 0;
                float poly[MAX_ROOTS + 1] = { -1, 0, 0, 0, 0, 1, 0 };
                float dx = 2, per = 20;
                float EPS = 1E-10;
                int i, deg = 5, NIT = 40;
                do
                {
                    printf("INPUT TO DRAW FRACTAL\n");
                    printf("1) Polynominal degree (max is %d) : \t%d\n", MAX_ROOTS, deg);
                    printf("   Coeficients : \n");
                    for (i = deg; i >= 0; i--) printf("\t\tA%d= %f\n", i, poly[i]);
                    printf("2) Width of the area on the cartesian surface (dx) : \t%f\n", dx);
                    printf("3) Percentage of the area of the screen to be covered: \t%f\n", per);
                    printf("4) Maximum Number of Iterations : \t\t\t%d\n", NIT);
                    printf("5) Precision : \t\t\t\t\t\t%G \n", EPS);
                    printf("6) Draw Fractal\n");
                    printf("7) Exit  \n");
                    printf("\nWhat is your choice ? ");
                    cin >> answ;
                    switch (answ)
                    {
                    case(1) :
                        do
                        {
                            printf("\nEnter the degree of the complex polynominal [<=%d]: ", MAX_ROOTS);
                            scanf("%d", &deg);//от 18 до 18
                        } while (deg > 15 || deg < 0);
                        printf("\nEnter the coeficients starting with the greatest power\n");
                        for (i = deg; i >= 0; i--)
                        {
                            printf("A%d = ", i);
                            scanf("%f", &poly[i]);
                        }
                        break;
                    case(2) :
                        do
                        {
                            printf("\nEnter the width (dx>0) = ");
                            scanf("%f", &dx);
                        } while (dx <= 0);
                        break;
                    case(3) :
                        do
                        {
                            printf("\nEnter the percentage of the screen to be covered by the fractal(%)\n");
                            printf("(the smallest the fastest ->useful for a preview) : ");
                            scanf("%f", &per);
                        } while (per <= 0 || per > 100);
                        break;
                    case(4) :
                        do
                        {
                            printf("\nEnter the maximum Number of ITerations (NIT>0) : ");
                            scanf("%d", &NIT);
                        } while (NIT <= 0);
                        break;
                    case(5) :
                        do
                        {
                            printf("\nEnter precision ( >0 ) : ");
                            scanf("%e", &EPS);
                        } while (EPS < 0);
                        break;
                    case(6) :
                        initwindow(700, 700);
                        drawFrac(dx, per, poly, deg, NIT, EPS);
                        WaitEscEnter();
                        closegraph();
                        break;
                    case(7) :
                        exit(0);
                    }
                } while (1);/*end switch */

setlocale(LC_ALL, "Russian");

case 2:
            {
                int colour3one, colour3two, colour3three, colour3four, colour3five;
                PrintColour();
                EnterColour(colour3one);
                EnterColour(colour3two);
                EnterColour(colour3three);
                EnterColour(colour3four);
                EnterColour(colour3five);
                int cycle = 1, cycle2 = 1;//Изменяем цвет в цикле
                double inclination = -32323.0;//Изменение наклона при каждой ветке
                cout << "Введите изменение наклона (от -1000.0 до 1000.0): ";
                while (!(cin >> inclination) || (cin.peek() != '\n') || (inclination <= -1001.0) || (inclination >= 1001.0))
                {
                    cin.clear();
                    while (cin.get() != '\n');
                    cout << "\a";
                    cout << "Введите изменение наклона (от -1000.0 до 1000.0): ";
                }
                initwindow(700, 700);
                for (int length = 90; length >= 30; length -= 20)
                {
                    if (cycle2 == 1) cycle = colour3one;
                    if (cycle2 == 2) cycle = colour3two;
                    if (cycle2 == 3) cycle = colour3three;
                    if (cycle2 == 4) cycle = colour3four;
                    zUnknown(350, 350, 270, length, 2, 20, cycle, inclination);
                    zUnknown(350, 350, 180, length, 2, 20, cycle, inclination);
                    zUnknown(350, 350, 90, length, 2, 20, cycle, inclination);
                    zUnknown(350, 350, 0, length, 2, 20, cycle, inclination);
                    cycle2++;
                }
                for (int corner = 315; corner >= 45; corner -= 90)
                {
                    zUnknown(350, 350, corner, 15, 2, 20, colour3five, inclination);
                }
                for (int corner = 270; corner >= 0; corner -= 90)
                {
                    zUnknown(350, 350, corner, 15, 2, 20, colour3five, inclination);
                }
                WaitEscEnter();
                closegraph();
                break;
            }

#include <windows.h>
 
int main()
{
  SetConsoleCP(1251);
  SetConsoleOutputCP(1251);
  ...
}