用matlab编写积分函数

function y = integral(f,a,b,N)
%输入参数
%f：被积函数
%a：积分下限
%b：积分上限
%N：划分子区间数
%输出参数
%y：积分值

%计算步长
dx = (b-a)/N;

%初始化积分结果
y = 0;

%循环计算每个子区间的积分值，并将结果累加到总积分值中
for i=1:N
x1 = a + (i-1)*dx;
x2 = a + i*dx;
y = y + dx*(f(x1)+f(x2))/2;
end

end