integral_richardson('I_S','sin(x)', 0, pi, 4, 3);