# shi(x) is the sinh integral, integral of sinh(x)/x; shi(0):=0; shi(inf):= inf;