log((1-b^2)^0.5 * ((1+b)/(1-b))^d), for 0<b<0.1, 0<d<0.5