11
$\开始组$

这个帖子,我提供了积分的解决方案$$\int_0^{\infty}\left[\int^1_0J_0(v\rho)\cos(\rho^2/2)\rho\,\mathrm{d}\rho\right]^2v\mathrm{d} v(v)$$

然而,我发现了一个错误总和在我的程序中,以下代码显示:

exp=((-1)^(j+k)2^(-1-2j-2k))/((1+2j+2k)伽马[1+2j]伽马[1x2k]);总和[exp,{j,0,\[Infinity]},{k,0,\[Infinity]}]%//否(*Sin[1]/4*)(*0.210368*)NSum[exp,{j,0,\[Infinity]},{k,0,\[Infinitiy]}](*0.460368-1.99638*10^-45 I*)

那么这是一个错误吗?

$\端组$
1
  • 1
    $\开始组$ 发件人漏洞tag wiki:“请不要用这个标签回答新问题。”但也就是说,我认为每个人都会认为这是一个错误,标签最终会被贴上。所以,除非有人反对,否则别管它。 $\端组$ 评论 5月26日19:13

2个答案2

重置为默认值
7
$\开始组$

在中多次总和[],n=0似乎是一个与bug有关的特例。分开j==0k==0:

经验=((-1)^(j+k)2^(-1-2j-2k))/((1+2 j+2 k)伽马[1+2 j]伽马[1+2 k]);求和[exp,{j,1,\[Infinity]},{k,1,\[Infinity]}]+2和[exp/.k->0,{j,0,无限}]-经验/。{j->0,k->0}//展开(*1/4+正弦[1]/4*)
$\端组$
2
  • $\开始组$ 谢谢你对我错误答案的评论。根据你的回答,我发现Sum[exp,{j,0,[Infinity]},{k,1,[Infinity]{]+Sum[exp/.k->0,{j、0,[Infinity]}]//展开也有效。尽管j和k可以互换,但外面的j=0并不能分开。因此,“bug”似乎来自第一个和(求和k)。 $\端组$ 评论 5月27日13:24
  • $\开始组$ @metroidman是的,有几个变种有效。IIRC,我试过了总和[exp,{j,1,无限},{k,0,无限}]+总和[exp/.k->0,{j首先,这与您的示例类似,但不起作用,并且支持您对第一个和的依赖性。然后我尝试了上面的方法。然后我摆弄各种组合。上面的和类似的求和[exp,{j,1,\[Infinity]},{k,1,\[Infinity]}]+求和[ex/.k->0,{j、1,Infinity}]+exp/。{j->0,k->0}//展开看起来速度最快。在这两个问题中,我不知道我为什么选择这个答案;事后看来,这似乎更为复杂。 $\端组$ 评论 5月27日14:26
5
$\开始组$

我想MMA会先做替换$j至n-k$,然后有一个错误总和对于$n=0$

expr=z^(2j+2 k)/(伽马[1+2 j]伽马[1+2 k])/.j->n-k//简化术语=总和[expr,{k,0,n}]表[{term,Sum[expr,{k,0,n}]},{n,0,4}]总和[expr,{k,0,n},假设->n==0]总和[expr,{k,0,n},假设->n>=0]/.n->0

在此处输入图像描述

总和[1/(伽马[1+2k]伽马[1-2k+n]),{k,0,n}]求和[z^(2k)/(伽玛[1+2 k]伽玛[1-2 k+n]),{k,0,n}]输出[]/.z->1

在此处输入图像描述

$\端组$

您的答案

单击“发布您的答案”,表示您同意我们的服务条款并确认您已阅读我们的隐私政策.

不是你想要的答案吗?浏览标记的其他问题问你自己的问题.