[討論] trapzoidal
小弟我上課作業要用trapzoidal來算積分
計算∫sin(x)/x 上下限是1和0
打出來的M.file
function I=trap(f,a,b,n)
x=a;
h=(b-a)/n;
s=feval(f,a);
for j=1:n-1
x=x+h;
s=s+2*feval(f,x);
end
s=s+2*feval(f,b);
I=(b-a)*s/(2*n);
在指令欄輸入
f=@(x)sin(x)/x
trap(f,0,1,4)
然後會出現
ans =
NaN
我不確定是我哪邊打錯了,我把function改成x或x^2/2等等比較簡單的函數,
都有辦法求解,
請問有神手能幫我看一下嗎?
感謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.124.45.196
推
04/15 12:53, , 1F
04/15 12:53, 1F
→
04/15 12:55, , 2F
04/15 12:55, 2F
→
04/15 12:56, , 3F
04/15 12:56, 3F
推
04/15 13:54, , 4F
04/15 13:54, 4F
推
04/15 14:40, , 5F
04/15 14:40, 5F
→
04/15 14:40, , 6F
04/15 14:40, 6F
→
04/15 15:16, , 7F
04/15 15:16, 7F
→
04/15 16:39, , 8F
04/15 16:39, 8F
→
04/17 05:34, , 9F
04/17 05:34, 9F
→
04/17 05:35, , 10F
04/17 05:35, 10F
推
04/17 12:26, , 11F
04/17 12:26, 11F
推
04/18 20:42, , 12F
04/18 20:42, 12F
→
04/18 20:43, , 13F
04/18 20:43, 13F
→
04/21 10:11, , 14F
04/21 10:11, 14F
→
04/21 10:12, , 15F
04/21 10:12, 15F
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
