wangsheng37 发表于 2012-8-17 10:08:45

请教有关Bessel积分

在论坛上看到有一个求Bessel积分的例子:exp(-x^2)*sin(x)*x*besselj(1,x),当时版主说可以用Mathematica计算比较合适,并给出了若x积分区间为 (0,inf),则结果为0.144742。但是由于本人目前对该软件不熟,试了下出现NIntegrate::nlim: x = inf is not a valid limit of integration.请教各位牛人这个积分应该怎么做?
ps:我用matlab的fplot看了下该函数的图像,应该是可以积出来的

TBE_Legend 发表于 2012-8-17 11:16:14

Infinity,不是inf。

wangsheng37 发表于 2012-8-17 16:49:30

谢谢版主,没想到和matlab还是有很大区别。最后在mathematica里面得到了结果0.144742。
f := Exp[-x^2]*Sin*x*Besselj;
NIntegrate, {x, 0, Infinity}];
1/4 HypergeometricPFQ[{5/4, 7/4}, {3/2, 5/2}, -1]
N]
0.144742
页: [1]
查看完整版本: 请教有关Bessel积分