RosettaCodeData/Task/Numerical-integration/MATLAB/numerical-integration-3.m
2023-07-01 13:44:08 -04:00

8 lines
231 B
Matlab

function integral = midPointRectIntegration(f,a,b,n)
format long;
width = (b-a)/n; %calculate the width of each devision
x = linspace(a,b,n); %define x-axis
integral = width * sum( f( (x(1:n-1)+x(2:n))/2 ) );
end