Hb
Hb = H - Ha;
if (Hb > 0.0) {
Hb -= 1.0/NXT; /*0.0625L;*/
z = Hb * __polevll(Hb, R, 6); /* z = 2**Hb - 1 */