# si(x) is the sin integral, integral of sin(x)/x; si(0):=0; si(inf):= pi/2; si(a_*inf):=pi*sgn(a)/2; si(1,0):=1;