Re: [問題] 這種算式metlab能算嗎?已回收
※ 引述《zupo (幫解MATLAB難題囧)》之銘言:
: ※ 引述《sale58 (便宜大碗)》之銘言:
: : f(x) = 4*(1 - 1/3 + 1/5 - 1/7 + 1/9 - 1/11 + ...)
: : 無限級數有辦法算出收斂值嗎?謝謝
: 我這算是硬幹法吧.
: 若板上有先進有更好的方法.歡迎提出
: -----------------------------------------------------
: function pttex113
: format long
: n = 10000000;
: f = zeros(n+1,1);
: f(1) = 4.*1;
: for h = 1:n
: f(h+1) = 4.*((-1).^(h)).*(1./(2.*h+1));
: end
: f = sum(f)
: ------------------------------------------------
: f 就是收斂值.
: 理論上會收斂到圓周率的值.
: 我是觀察到的啦
: 以下是結果
: f =
: 3.141592753589781
不好意思想問一下 pttex113是內建的m函數嗎?
n是指算到第1000萬項嗎?
zeros(n+1,1) 這個指令的功能?
f(1) = 4.*1;
: for h = 1:n
: f(h+1) = 4.*((-1).^(h)).*(1./(2.*h+1));
: end
這個迴圈也看不太懂,不好意思我初學者謝謝
剛剛跑了一下我出現這結果
>> function pttex113
format long
n = 10000000;
f = zeros(n+1,1);
f(1) = 4.*1;
for h = 1:n
f(h+1) = 4.*((-1).^(h)).*(1./(2.*h+1));
end
f = sum(f)
??? function pttex113
|
Error: Function definitions are not permitted at the prompt or in scripts.
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.152.193
討論串 (同標題文章)
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章