Подразумеваемая волатильность модели SABR в VBA

Question

Я пытаюсь написать код, который будет рассчитывать подразумеваемую волатильность в соответствии с моделью SABR с учетом указанных параметров. Тем не менее, я всегда получаю NA для ATM vol. Используется следующая формула подразумеваемой волатильности:

Function SABRIV(alpha As Double, beta As Double, v As Double, rho As Double, f As Double, k As Double, T As Double) As Double
Dim z As Double
Dim x As Double
Dim numerator As Double
Dim denominator As Double

z = (v / alpha) * ((f * k) ^ ((1 - beta) / 2)) * Application.Ln(f / k)
x = Application.Ln((Sqr(1 - 2 * rho * z + z ^ 2) + z - rho) / (1 - rho))
numerator = ((1 - beta) ^ 2 / 24) * ((alpha ^ 2) / ((f * k) ^ (1 - beta))) + 1 / 4 * (rho * beta * v * alpha) / ((f * k) ^ ((1 - beta) / 2)) + (((2 - 3 * rho ^ 2) * v ^ 2) / 24)
denominator = (((f * k) ^ ((1 - beta) / 2)) * (1 + (((1 - beta) ^ 2) / 24) * (Application.Ln(f / k) ^ 2) + (((1 - beta) ^ 4) / 1920) * (Application.Ln(f / k) ^ 4)))

SABRIV = (alpha * (1 + numerator) * T) / denominator * z / x
End Function