Urgent help with C++ program that uses the trapezoidal method for finding integrals.?

integral(sin(t)/t) is the "Sine Integral." The expression used in the trapezoidal method is relatively straight forward, I'd just hard-code the double derivative of the sine integral.