Как найти все корни нелинейного уравнение f(x)=0 на отрезке x∈[a, b] прямой параболической интерполяцией?

using System;
using System.IO;
using System.Numerics;
 
namespace PriamaiaParabolicheskaiaInterpoliacia
{
    internal class Program
    {
        public delegate double Function(double x);
 
        static double F(double x)
        {
            return Math.Pow(x, 2) + 4*x + 1; // Пример функции
        }
 
        static double DivadedDifference_1(double x0, double x1)
        {
            return F(x1) / (x1 - x0) + F(x0) / (x0 - x1);
        }
 
        static double DivadedDifference_2(double x0, double x1, double x2)
        {
            return F(x2) / ((x2 - x1) * (x2 - x0)) + F(x1) / ((x1 - x2) * (x1 - x0)) + F(x0) / ((x0 - x2) * (x0 - x1));
        }
 
        static Complex Interpolation(Function f, double x0, double x1, double x2, double epsilon)
        {
            Complex root = x2;
 
            while (true)
            {
                
                double a = DivadedDifference_2(x0, x1, x2);
                double b = DivadedDifference_1(x0, x1) + (x0-x1) * DivadedDifference_2(x0, x1, x2);
                double c = f(x0);
 
                
                double discriminant = b * b - 4 * a * c;
 
                Complex z1, z2;
 
                if (discriminant < 0)
                {
                    
                    z1 = (-b + Complex.Sqrt(discriminant)) / (2 * a);
                    z2 = (-b - Complex.Sqrt(discriminant)) / (2 * a);
                }
                else
                {
                    
                    z1 = (-b + Math.Sqrt(discriminant)) / (2 * a);
                    z2 = (-b - Math.Sqrt(discriminant)) / (2 * a);
                }
 
                Complex z = Complex.Abs(z1) < Complex.Abs(z2) ? z1 : z2;
                Complex nextRoot = x2 + z;
 
                
                if (Complex.Abs(nextRoot - root) < epsilon)
                    return nextRoot;
 
                
                x0 = x1;
                x1 = x2;
                x2 = nextRoot.Real;
                root = nextRoot;
            }
        }
 
        static void Main(string[] args)
       {
            double a=-10, b=10;
            Complex root = Interpolation(F, a, (a + b) / 2, b, 0.0001);
 
            Console.WriteLine("Найденный корень: " + root);
            Console.WriteLine("Проверка решения: f(x) = " + F(root.Real).ToString());
       }
}