如何使用matlab编写求积分的程序

假设以矩形或梯形面积来拟合吧带慎键；

function yanshi(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

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))

xx=[xi(i) xi(i) xi(i+1) xi(i+1) xi(i)]

yy=[0 syi syi 0 0]

fill(xx,yy,'c')

hold on

end

x=a:(b-a)/100:b

y=feval(fname,x)

plot(x,y,'r','markersize',20)

fprintf('n=%6d I=%12.5f\n',n,I)

hold off

function y=lirui(x)

y=1+x.^2

t = (1/16:1/8:1)'*2*pi

x = sin(t)

y = cos(t)

fill(x,y,'r')

axis square

感觉不需亮大要使用m文件，根据你洞和给的代码，函数f(x)的参数x没有作用啊。具纳键盯体需求是什么？

下面是求积分

f=@(x) exp(0.5*x).*sin(x+pi/6)

s=quad(f,0,pi)

运用下列代码，就可以得到cos⁴x的积分肢历唯值。

在matlab的命令窗口中，输入

>>syms x

>>int((cos(x))^4)

执历培行后，得到如下结果烂如

注意：不定积分计算后还得加积分常数 C

---