return (1.0 - float_power(h, -k))/k