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