# ci(x) is the cos integral, integral of cos(x)/x; ci(0):= -inf; ci(inf):=0;