#	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;