implementation rules Midpoint rule M(f)=(b−a)f(2a+b) Trapezoidal rule T(f)=2b−af(f(a)+f(b)) Simpson’s rule S(f)=6b−a(f(a)+5f(m)+f(b)) exact for n node up to n degree