2d trapezoidal rule python