Re: [請益] 請問有辦法這樣算嗎？

看板Mathematica作者chungyuandye (養花種魚數月亮賞星星)時間11年前 (2013/01/28 22:07)推噓1(1推 0噓 0→)

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

※ 引述《van6499 (KIDD)》之銘言： : 標題: [請益] 請問有辦法這樣算嗎？ : 時間: Mon Jan 28 19:19:35 2013 : : : : : 小弟手頭有兩個式子需要作圖 : : 其中一個帶有0到無限大的積分 : : 另外一個帶有Sigma 0到無限大的加總 : : : 這兩條式子已經確定會收斂，但是用PLOT指令下去跑圖真的跑得非常非常久 : : 想了一下，想用其他方式代替完整的函數圖 : : 請問有辦法做出下面兩種嗎？ : : : 1. 給定指定的X值群，讓它跑出相對應的函數解嗎？ : : EX：給定X=0.1 0.2 0.3 0.4 ....... 100.1 ...149.5 150 : : 跳出以上數列相對應的解 : (* 可以觀察出Mathematica Plot的取點方式與上面的方法有何差異 *) pts = Plot[Sin[x], {x, 0, 2 Pi}][[1, 1, -1, 2, 1]]; Plot[Sin[x], {x, 0, 2 Pi}, Epilog -> {Red, Point@pts}] : : 2. 輸入一群X.Y值，做出XY平面的點圖 : : : : 或者請問各位大神有更好的方法嗎？ : : Clear["Global`*"]; part1[x_?NumberQ] := 2 Pi*0.0512*1*195*4*(0.0512 x)/(0.0001*36)* Quiet@NIntegrate[ u*Exp[-0.0512*x/0.0001/36*u^2]*(Pi/2 + ArcTan[BesselY[0, u]/BesselJ[0, u]]), {u, 0, 10000}] Plot[part1[x], {x, 1/24, 250}, PlotRange -> {0, 100}] ff[z_] = 2/Sqrt[Pi]*Integrate[Exp[-a^2], {a, z, Infinity}]; (* 字太小，括號後半段沒打 *) part21[x_] = 195/((1/(2 Pi*0.0512)*(Sum[ ff[2 Pi/Sqrt[4*512*x]]/(2*n*1), {n, 1, 2000}] + Sum[ff[2 Pi/Sqrt[4*512*x]]/(2*n*1), {n, -1, -2000, -1}]) + ff[6/Sqrt[4*512*x]]/6)); -- 養花種魚數月亮賞星星 http://chungyuandye.twbbs.org -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 218.173.134.241