求解一个定积分的题

lxx364282219

我的程序已有，但有错误希望哪位大侠，帮忙修改一下，我是新手，不是数学专业，但要交实验！！！！
：设 y=x.^3-3*x.^2+x+5  ，根据定积分定义编写一段程序，从几何上演示用小矩形面积和逼近曲边梯形面积的过程。
function juxs(fname,a,b,n) % 定积分概念演示，随机分割、随机取近似，并求近似值
xi(1)=a; xi(n+1)=b;
for i=1:n-1
xi(i+1)=a+(i+rand(1))*(b-a)/n;
end
I=0;
hold on;
for i=1:length(xi)-1
sxi=xi(i)+rand(1)*(xi(i+1)-xi(i));
syi=feval(fname,sxi);
I=I+syi*(xi(i+1)-xi(i));
xii=[xi(i) xi(i) xi(i+1) xi(i+1) xi(i)];
yii=[0 syi syi 0 0];
fill(xii,yii,'c');
end
x=a:(b-a)/100:b;
y=feval(fname,x);
plot(x,y,'r','markersize',20);
hold off;
fprintf('n=%6d, I=%12.5f\n',n,I);
function y=fex(x)
y=x.^3-3*x.^2+x+5；
而且，不会调用，帮忙怎么调用，谢谢。

赤血冰霜

我也不会，想学习，望会的人快来光顾啊！