Как найти все корни уравнения?

#include <stdio.h>
#include <math.h>
 
const int A = 5;
const int B = 1;
const int C = -3;
const int D = 3;
const double EPS = 0.000001;
 
 
 
double f(double x)
{
   return A * pow(x,3) + B * x * x + C * x + D;
}
 
double diff(double x)
{
    return 15 * x * x + 2 * x - 3;
}
double bisection(double xl, double xr, double EPS)
{
    int n;
    if (f(xl)*f(xr)>0)
                  printf("\nError! No roots in this interval\n");
    else
    {
        int n = 0;
        double xd = xr - xl;
        while ( fabs(f(xl))>EPS || fabs(f(xr))>EPS )
        {
        n = n + 1;
        xd = xd / 2;
            double xm = xl + xd;
        if ( f(xl) *  f(xm) <= 0 )
             xr = xm;
        else
             xl = xm;
        //printf ("\nNumbers of iterations equal: %d x0 = %lf\n", n, xm);
        }
    }
 
    return (xl + xr) / 2;
}
int main ( )
{
   double xl, xr, xm;
   int n;
   printf("Enter the interval");
   printf("\nLeft  xl = "); scanf("%lf", &xl);
   printf("\nRight xr = "); scanf("%lf", &xr);
   if (xl > xr || xl == xr)
   {
       printf("Invalid data");
   }
   double r = bisection(xl, xr, EPS);
   if (f(xl)*f(xr)>0)
                 printf("");
   else
   {
   if (diff(xl) < 0 || diff(xr) > 0)
       printf ("Infimum. Middle of line segment: %.10lf\n", r);
   else
       printf ("Supremum. Middle of line segment: %.10lf\n", r);
   printf ("Value of function: %.10lf\n", f(r));
   }
}

double s = 1; //шаг
double i, j;
for (i = xl; i <= xr; i++)
    {
      for(j = xl + s; j <= xr; j++)
       {
          if (f(i) * f(j) < 0)
            {
                double r = bisection(i, j, EPS);
                printf ("Middle of line segment: %lf\n", r);
                printf ("Value of function: %lf\n", f(r));
            }
       }
    }

double i, j;
int d = 0;
 
for (i = xl; i <= xr; i++)
    {
      for(j = xl + 1; j <= xr; j++)
       {
          if (f(i) * f(j) < 0)
            {
                d = 1;
                double r = bisection(i, j, EPS);
            }
       }
if (d == 1)
printf ("Middle of line segment: %lf\n", r);
else и т. д.
 
 
    }
 
return 0;
}

for (i = xl; i < xr; i++)
{
   if (f(i) * f(i+1) < 0)
      double r = bisection(i, i+1, EPS);
}