[xmuhjh]

可以参考如下代码：
给出三种方法，前两种计算的是积分结果；
第三种 cumtrapz 计算积分累积分布。

也可以这样考虑：用trapz或者 quad quadl quad8逐点计算出积分结果，然后把结果画出来。

clear;
x=linspace(0,3,30);
y=x;
inter1=trapz(x,y);

inter2=quadl(@integrand,0,3);

inter3=cumtrapz(x,y);

plot(x,inter3);
hold on;

plot(x,y);


function y=integrand(x)

y=x;