Trapezoidal Method Matlab Code