При выводе результата работы программы пишет #QNAN0

#include <stdio.h>
#include <math.h>
#include <iostream>
using namespace std;
 
 
 
double f(double x)
{
    int A = 2;
    int B = -3;
    int C = -4;
    int D = 2;
    return A * pow(x,3) + B * x * x + C * x + D;
}
 
/*
double df(double x)
{
    int A = 2;
    int B = -3;
    int C = -4;
    int D = 2;
    return A * 3 * x * x + B* 2 * x - C;
}
 
 
double df2(double x)
{
    int A = 2;
    int B = -3;
    int C = -4;
    int D = 2;
    return A * 3 * 2 * x + B * 2;
}
*/
 
double b(double xl, double xr)
{
    if (f(xl) * f(xr) > 0)
        printf("\nError! No roots in this interval\n");
    else
    {
        double xd = xr - xl;
        while ( fabs(f(xl)) > 0.001 || fabs(f(xr)) > 0.001 )
        {
            xd = xd / 2;
            double xm = xl + xd;
            if ( f(xl) *  f(xm) <= 0 )
                xr = xm;
            else
                xl = xm;
          //printf ("\nxm = %lf", xm); //Value of the root after each iteration
        }
    }
    return (xl + xr) / 2;
}
 
 
double koren(double (*f)(double), double *a, double xl, double xr)
{
    for(int j = 0; j < 3; j++)
        {
            for (double i = xl; i < xr; i++)
                {
                    if (f(i) * f(i+1) <= 0)
                        {
                                a[j] = b(i, i + 1);
                        }
                }
        }
}
 
 
int main ( )
{
   double xl, xr;
   double *a;
   a = new double [3];
   printf("Enter the interval");
   printf("\nLeft  xl = "); scanf("%lf", &xl);
   printf("\nRight xr = "); scanf("%lf", &xr);
   printf("\n\nROOTS:\n\n");
   if (xl > xr || xl == xr)
   {
       printf("Invalid data");
   }
    else
        {
            for(int j = 0; j < 3; j++)
            {
                double d = koren(f, a, xl, xr);
                printf ("\nx[%d] = %lf", j, d);
            }
        }
   delete[] a;
   return 0;
}