void bisection(vector<double>& results,double a, double b)
{
    int iteration_count = 0;
    double middle_point;
    while (fabs(a-b) > EPSILON)
    {
        if (iteration_count == NMAX)
            break;
        middle_point = (a + b) / 2;
        if (function(a) * function(middle_point) < 0)
            b = middle_point;
        else
            a = middle_point;
        iteration_count++;
        if (function(middle_point) == 0)
            break;
    }
    results.push_back(middle_point);

}

main:
    while (b != B_POINT)
    {
        bisection(results, a, b);
        a += 0.5;
        b += 0.5;
    }