k = 0.69502
constant = 0.3751349
zfsD = Op / T / k
topp = 0.75 * (1 + 9 * Exp(-2 * zfsD) + 25 * Exp(-6 * zfsD))
toppZ = 0.75 * (9 + 8 / zfsD - 11 * Exp(-2 * zfsD) / (2 * zfsD) - 5 * Exp(-6 * zfsD) / (2 * zfsD))
bott = 1 + Exp(-2 * zfsD) + Exp(-6 * zfsD)
CHI = (topp / bott + 2 *toppZ / bott) / 3
CHI = CHI * constant * g * g / (3 * (T - th))+tip