trapezoidal rule code snippets in matlab