series sinh(log(a/b)*p)/(a-b) at a=b