Нахождение корня уравнения методом половинного деления

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
 
namespace ConsoleApplication1
{
    class Program
    {
        delegate double fx(double r, double[] Q, double[] P, double[] V, double fs, int n, double s, double Io);
        static double f(double r, double[] Q, double[] P, double[] V, double fs, int n, double s, double Io)
        {
            double G;
            for (int t = 1; t <= n; t++)
            {
                G = (Q[t - 1] * (P[t - 1] - V[t - 1]) - 800 - 300) * (1 - 0.6) + 300;
                s = s + G * Math.Pow(1 + r, n - t);
            }
            Console.WriteLine(s);
            return fs = s + 200 / Math.Pow(1 + r, n) - Io * Math.Pow(1 + r, n);
        }
        static double Poldel(fx f, double a, double b, double[] Q, double[] P, double[] V, double fs, int n, double s, double Io)
        {
            double xt;
            do
            {
                xt = (a + b) / 2.0;
                if (f(a, Q, P, V, fs, n, s, Io) * f(xt, Q, P, V, fs, n, s, Io) < 0)
                    b = xt;
                else
                    a = xt;
            }
            while (b - a > 1e-5);
            return xt;
        }
        static void Main(string[] args)
        {
            double[] Q1 = new double[6] { 150, 170, 185, 210, 240, 280 };
            double[] Q2 = new double[6] { 153, 177, 209, 228, 265, 288 };
            double[] Q3 = new double[6] { 158, 181, 189, 215, 265, 285 };
            double[] P1 = new double[6] { 35, 37, 41, 44, 48, 52 };
            double[] P2 = new double[6] { 35, 39, 46, 49, 52, 54 };
            double[] P3 = new double[6] { 38, 40, 47, 49, 50, 52 };
            double[] V1 = new double[6] { 25, 29, 31, 33, 35, 40 };
            double[] V2 = new double[6] { 25, 29, 31, 34, 38, 40 };
            double[] V3 = new double[6] { 25, 32, 34, 36, 38, 40 };
            double Io1 = 5500, Io2 = 8000, Io3 = 7800;
            double s = 0, fs = 0;
            int n = 6;
            double r1, r2, r3;
            r1 = Poldel(f, -1000, 1000, Q1, P1, V1, fs, n, s, Io1);
            r2 = Poldel(f, 0, 1, Q2, P2, V2, fs, n, s, Io2);
            r3 = Poldel(f, 0, 1, Q3, P3, V3, fs, n, s, Io3);
            Console.WriteLine("Коэффициенты дисконтирования r1 = {0:f2}%    r2 = {1:f2}%    r3 = {2:f2}% ", r1 * 100, r2 * 100, r3 * 100);
            Console.ReadKey();
        }
    }
}