Некорректно работает код

#define _CRT_SECURE_NO_WARNINGS
#include <cstdlib>
#include <math.h>
#include <cstdio>
#include <conio.h>
#include <Windows.h>
#include <iostream>
 
using namespace std;
 
double abc(double, double);
 
 
 
int main()
{
    double x_1, x_2, shag, pi, x_id, test;
    double myfunc, func, E;
 
    pi = 3.14;
 
    printf("Vvedite nachal'noe znachenie X: \n");
    cin >> x_1;
    while (x_1 * x_1 > pi * pi)
    {
        printf("Oshibka, povtorite vvod \n");
        cin >> x_1;
    }
 
    printf("Vvedite konechnoe znachenie X: \n");
    cin >> x_2;
    while (x_2 < 0 || x_2 <= x_1)
    {
        printf("Oshibka, povtorite vvod \n");
        cin >> x_2;
    }
 
    printf("Vvedite shag dlya X: \n");
    cin >> shag;
    while (shag < 0)
    {
        printf("Oshibka, povtorite vvod \n");
        cin >> shag;
    }
 
    printf("Vvedite toshnost': \n");
    cin >> E;
    while (E = 0 || E < 0)
    {
        printf("Oshibka, povtorite vvod \n");
        cin >> E;
    }
    printf(" _______________________________________________________________\n");
    printf("|               |               |               |               |\n");
    printf("|       X       |    My Func    |     Func      |     Sigma     |\n");
    printf("|_______________|_______________|_______________|_______________|\n");
 
    for (x_1; x_1 <= x_2; x_1 += shag) {
 
        x_1 = 10000000 * x_1;
        x_1 = round(x_1);
        x_1 = x_1 / 10000000;
        myfunc = x_1 * x_1 * abc(x_1, E) * pi / 180;
 
        if (fmod(x_1, 2.0) == 1.0) {
            cout << "|               |               |               |               |\n";
            printf("|%-15lf|%-15lf|error          |error          |\n", x_1, myfunc);
            printf("|_______________|_______________|_______________|_______________|\n");
 
        }
        else {
            func = tan((pi * x_1) / 2) * tan((pi * x_1) / 2) * pi / 180;
            cout << "|               |               |               |               |\n";
            printf("|%-15lf|%-15lf|%-15lf|%-15lf|\n", x_1, myfunc, func, sqrt(fabs(myfunc * myfunc - func * func)));
            printf("|_______________|_______________|_______________|_______________|\n");
 
        }
 
    }
    printf("\nVvedite ideal'noe znachenie X: ");
    cin >> x_id;
    printf(" _______________________________________________________________\n");
    printf("|               |               |               |               |\n");
    printf("|Eps            |My Func        |Func           |sigma          |\n");
    printf("|_______________|_______________|_______________|_______________|\n");
    for (E = 0.1; E >= 0.0000001; E /= 10)
    {
        cout << E << "\n";
        myfunc = x_id * x_id * abc(x_id, E) * pi / 180;
        func = tan((pi * x_id) / 2) * tan((pi * x_id) / 2) * pi / 180;
        printf("|               |               |               |               |\n");
        printf("|%.7lf      |%-15lf|%-15lf|%-15lf|\n", E, myfunc, func, sqrt(fabs(myfunc * myfunc - func * func)));
        printf("|_______________|_______________|_______________|_______________|\n");
    }
    return 0;
}
 
    double abc(double x, double E)
    {
        double k, i = 1;
        double bigsum = 0;
        double curzn;
        double chislitel;
        double znamenatel = 1;
        double drob = 1000.1;
        double prevdrob = 0;
 
        for (k = 1; fabsl(drob - prevdrob) >= E; k++) {
            chislitel = 2 * (2 * k - 1) * (2 * k - 1) - x * x;
            while (i <= (2 * k - 1)) {
                if (i == 1) {
                    curzn = i * i - x * x;
                }
                else {
                    curzn = (i * i - x * x) * (i * i - x * x);
                }
                znamenatel = znamenatel * curzn;
                i = i + 2;
            }
            drob = chislitel / znamenatel;
            bigsum = bigsum + drob;
            prevdrob = drob;
        }
 
        return bigsum;
    }