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