作者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