'   n As Integer      число элементов векторов дат и сумм
'   dt() As Date      вектор дат
'   db() As Double    вектор сумм
'   sMethod As String метод расчета процентов
Public Function XIRR(n As Integer, dt() As Date, db() As Double, sMethod As String) As Variant

    Dim t() As Integer ' Число дней от первоначальной даты
    Dim dd() As Double ' Дисконтированные к начальной дате суммы
    
    ' Выделение места в памяти для нужного количества элементов
    ReDim t(0 To n) As Integer
    ReDim dd(0 To n) As Double
    
    Dim iLoop As Integer, iMax As Integer: iMax = 100
    Dim x1 As Double, x2 As Double, xMid As Double, _
        y1 As Double, y2 As Double, yMid
        
    ' Заполнение массива число дней от начальной даты до очередного CashFlow
    Dim i As Integer
    For i =  0  To n
        ' Метод расчета процентов
        '  360 / 360  Месяц/квартал/полугодие/год =  30 / 90 / 180 / 360  дней
        ' 365/360 Фактический интервал / 360 дней в году
        '  365 / 365  Фактический интервал /  365  дней в году
        t(i) = DateDiffMethod(dt( 0 ), dt(i), sMethod)
        '        Debug.Print "i=", i, "t(i)=", t(i)
    Next i
    
    ' Инициализация переменных для поиска решения методом половинного деления
    ' x1 и x2 - текущие значения аргумента, определяющие границы отрезка
    ' y1 и y2 - соотвествующие им значения функции
    ' xMid и yMid - середина отрезка и значение функции в этой точке
    x1 = 0.1: x2 = x1 + (x1 / 2)
    y1 = 0:   y2 = 0
    
    ' Первоначальной расчет значения функции на границе отрезка x1
    For i =  0  To n
        dd(i) = db(i) / (( 1  + x1) ^ (t(i) / DateDiffBasis(sMethod)))
        y1 = y1 + dd(i)
        '        Debug.Print "i=", i, "db(i)=", db(i), "t(i)=", t(i), "dd(i)=", dd(i)
    Next i
    Debug.Print y1
    iLoop = 0
    
NEXT_LOOP:
    ' Расчет значения функции на границе отрезка x2
    y2 =  0 
    For i =  0  To n
        dd(i) = db(i) / (( 1  + x2) ^ (t(i) / DateDiffBasis(sMethod)))
        y2 = y2 + dd(i)
    Next i
    iLoop = iLoop +  1 
    If iLoop >= iMax Then GoTo NOT_FOUND
    ' Расчет значения функции на границе отрезка x2
    ' Если между двумя точками x1 и x2 нет решения,
    ' то двигаться увеличивающимися шагами до тех пор,
    ' пока решение функции не окажется между границами
    If y1 * y2 >  0  And Abs(y2) < Abs(y1) Then
        ' Определение направления движения в строну от х1 к x2
        '         Debug.Print "y1 * y2 > 0 And Abs(y2) < Abs(y1)", _
        '         "R1=" & x1, "R2=" & x2, "S1=" & y1, "S2=" & y2
        Dim x As Double
        x = x2: x2 = x1 + (x2 - x1) * 2: x1 = x
        y1 = y2
        GoTo NEXT_LOOP
    ElseIf y1 * y2 > 0 And Abs(y2) > Abs(y1) Then
        ' Определение направления движения в обратную строну от х2 к x1
        '         Debug.Print "y1 * y2 > 0 And Abs(y2) > Abs(y1)", _
        '         "R1=" & x1, "R2=" & x2, "S1=" & y1, "S2=" & y2
        x2 = x1 - (x2 - x1)
        GoTo NEXT_LOOP
    Else
        ' Поиск решение между границами x1 и x2
        '         Debug.Print "y1*y2 < 0            ", _
        '         "R1=" & x1, "R2=" & x2, "S1=" & y1, "S2=" & y2
NEXT_LOOP2:
        If Abs(x2 - x1) < 0.000000000001 Then
            GoTo FOUND
        End If
        ' Расчет средней точки и значения функции в ней
        xMid = (x1 + x2) /  2 
        yMid =  0 
        For i =  0  To n
            dd(i) = db(i) / (( 1  + xMid) ^ (t(i) / DateDiffBasis(sMethod)))
            yMid = yMid + dd(i)
        Next i
        '         Debug.Print "                     ", _
        '         "R1=" & x1, "R2=" & x2, "S1=" & y1, "S2=" & y2
        ' Выбор половины отрезка, в котором будет производиться следующая итерация
        If y1 * yMid >  0  Then
            x1 = xMid
            y1 = yMid
            GoTo NEXT_LOOP2
        ElseIf y2 * yMid >  0  Then
            x2 = xMid
            y2 = yMid
            GoTo NEXT_LOOP2
        Else
            GoTo FOUND
        End If
    End If

FOUND:
    XIRR = x1
    Exit Function
    
NOT_FOUND:
    XIRR = Null
    Exit Function

End Function