Create a UDF
Public Function Qr(Pr, Pbh, Ql, Pbp) As Double
Dim J As Double, Qm As Double, Qbp As Double
If Pr <= Pbp Then
Qm = Ql / (1 - 0.2 * (Pbh / Pr) - 0.8 * (Pbh / Pr) ^ 2)
Qr = Qm * (1 - 0.2 * (Pbh / Pr) - 0.8 * (Pbh / Pr) ^ 2)
ElseIf Pbh >= Pbp Then
J = Ql / (Pr - Pbh)
Else
J = Ql / (Pr - Pbp) + (Pbp / 1.8) * (1 - 0.2 * (Pbh / Pbp) - 0.8 * (Pbh / Pbp) ^ 2)
Qbp = J * (Pr - Pbp)
If Pbh >= Pbp Then
Qr = J * (Pr - Pbh)
Else
Qr = Qbp + (J * (Pbp / 1.8) * (1 - 0.2 * (Pbh / Pbp) - 0.8 * (Pbh / Pbp) ^ 2))
End If
End If
End Function