Function File: y = Si (x)

Compute the sine integral defined by:

                        x
                       /
               Si(x) = | sin(t)/t dt
                       /
                       0