pow1p function

Accurate (1+x)y(1+x)^y, notably for small x|x|