※ 引述《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));