return (x - loc) * c / expm1(-c * logsf)