integrate 1+x^2 from -1 to 3 using right endpoint method with 4 intervals