Composite Trapezoid Rule: Difference between revisions

m
No edit summary
Line 5:
 
function integral = trapezoidal(f, a, b, n)
hx = (b-a)/n;
result = 0.5*f(a) + 0.5*f(b);
for i = 1:(n-1)
result = result + f(a + i*hx);
end
integral = hx*result;
end