作者evovss48 (疯狂假面)
看板MATLAB
标题Re: [运算] 梯形法
时间Sun Dec 4 01:10:18 2011
※ 引述《handsomepow (handsomepow)》之铭言：
: 小弟我刚学MATLAB
: 我想用梯形法求exp(-x^2)积分 [0,2] 切n=500
: 不知道哪里出错了 麻烦高手指导一下 谢谢
: a=0;
: b=2;
: n=500;
: A=0;
: for x=a:(b-a)/n:b
: if x==a
: f1= exp(-x^2);
: else
: f2 = exp(-x^2);
: A = A+0.5*(f1+f2)*(b-a)/n;
: f1 = f2;
: end
: end
: disp('梯形积分法面积=');
: A
可以简化成
x=linspace(0,2,500);
y=exp(-x.^2);
Ans=trapz(x,y);
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 115.43.218.238
1F：推 handsomepow:题目要我把把梯形法写成function所以我想说看能不能跑 12/04 09:49