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