Re: [問題] 這種算式metlab能算嗎?已回收

看板MATLAB作者Leon (Achilles)時間16年前 (2008/10/12 03:23)推噓0(0推 0噓 1→)

留言1則, 1人參與討論串6/6 (看更多)

※ 引述《zupo (幫解MATLAB難題囧)》之銘言： : : n是指算到第1000萬項嗎? : ^^^^^^^^^^^^^^^^^^^^^^ : 是的.數字這麼大的理由.當然是希望能逼近到無窮大. : 但是這個題目項數與計算式有關.用 inf 來當矩陣的 : 位置是不可能的.我採取一個電腦可以承受計算的數字呃, 各位, 在使用計算機的暴力之前, 請先想想問題的本質. 這個級數收斂嘛? 你要先證明收歛才有意義. and it's not so hard to find, S = 1 - 1/3 + 1/5 .... = sum_{k} (-1)^(k+1) * 1/(2*k -1 ) = sum_{k} 2/ (4*k+1)(4*k+3) , Since sum 1/(k*k) converge, the series S should converge. 另外提供一個解析解: S = 1 -1/3 + 1/5 .... 這叫 Leibniz formula for pi Consider 1 - x^2 + x^4 ... = 1/(1+x^2) , geometric series, Integrate both parts x - x^3/3 + x^5/5 + ... = inverse tan(x). let x =1. 1- 1/3 + 1/5 .... = inverse tan(1) = pi/ 4. Q.E.D : : zeros(n+1,1) 這個指令的功能? : ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ : 計算結果矩陣內的項數高達一千萬零一項.不先產生個預設為零 : 一千萬零一*1矩陣的結果.就是計算牛步化. : zeros函數是用來產生零矩陣的MATLAB內建函數 : : f(1) = 4.*1; : : 這個迴圈也看不太懂,不好意思我初學者謝謝 : : 剛剛跑了一下我出現這結果 : : 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. : zupo : .............. : 看來你該是買本入門的書來唸的時候了. : 我的程式直接執行就好.沒有輸出入引數的 : 你把整段程式碼貼到執行視窗內當然無法執行 : 是要用 m-file edit 上面點執行來執行即可 -- 趙客縵胡纓，吾鉤霜雪明。銀鞍照白馬，颯沓如流星。十步殺一人，千里不留行。是了拂衣去，深藏身與名。閑過信陵飲，脫劍膝前橫。將炙啖朱亥，持觴勸侯贏。三杯吐然諾，五嶽倒為輕。眼花耳熱後，意氣素霓生。就趙揮金錘，邯鄲先震驚。千秋二壯士，烜赫大梁城。 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 76.170.72.148