Function s=traprl(f, a, b, M) oInput-f is the integrand input as a string f a and b are upper and lower limits of integration M is the number of subintervals %Output -s is the tra pezoidal rule sum h=(b-a)/M: 0: for k=1: (M-1) x=a+h*k s=s+feva(f,×) en s=h*(feval(f, a)+feval(, b)/2+h*s