out2]=4x∥Re[]>0 4235e2b 输出的含义是:当ReA>0时,计算结果为x/3,否则, Mathematica无法求解,将返回输 入形式4[e2br2b。完整的结果应是: 下一步,我们将借助引入的“试验”波函数求动能项的期望值。由于我们只讨论基态的能 量本征值,而对基态量子数l=0,此时在径向中心力场势情况下可采用拉普拉斯算子形式 为 △=V2d2,2d 其期望值为 =24xa(2x2-22)2 我们可以看到这里的量纲检验仍然是正确的。我们在下式中省略了Dm[…]符号: 平平→E-E32E2E3=E2→22 动能项的期望值为 p=(dNl, H(,) (82.12) 相应的 Mathematica v4.0计算过程为 MATHEMATICA V4. 0 In[3]:=vx= (*定义“试验”波函数*) In4]:=s,]=D,2]+2pr,l明 (*轨道角动量为零的有效拉普拉斯算符* 4[rwr, a][r, a]drOut[2]= , ] 4 1 4 [Re[ ] 0, 0 2 2 3 If e r dr r ∫ ∞ − > λ λ π λ 输出的含义是：当Reλ > 0 时，计算结果为 3 π λ ,否则，Mathematica 无法求解，将返回输 入形式 dr 。完整的结果应是： 2 e r r 0 2 4 ∫ ∞ − λ π 1 4 4 3 2 = λ π N ⇒ π λ 3 N = . 下一步，我们将借助引入的“试验”波函数求动能项的期望值。由于我们只讨论基态的能 量本征值，而对基态量子数l ，此时在径向中心力场势情况下可采用拉普拉斯算子形式 为 = 0 dr d dr r d 2 2 2 2 ∆ ≡ ∇ = + . 其期望值为 ( ) ( ) r r e dr d dr r d d x r r d xe λ λ π λ λ λ − −         Ψ ∆Ψ = + ∫ 2 , , 2 2 3 3 3 * ∫ ( ) r dr r r e λ π λ λ π λ 2 0 2 2 3 4 2 − ∞ ∫ = − ( ) ( ) ( ) ( )       Γ − Γ = 3 2 3 2 2 2 2 2 3 4 λ λ λ λ λ 3 2 4 1 4 λ λ λ  = −      = − . 我们可以看到这里的量纲检验仍然是正确的。我们在下式中省略了 Dim[…]符号: 3 3 / 2 2 3 / 2 2 2 2 3 1 Ψ Ψ → = → λ − E E E E E x x . 动能项的期望值为 ( ) ( ) µ λ λ µ λ µ 2 , 2 , 2 2 3 * 2 Ψ =         ∆ = Ψ − ∫ d x r r p G . (8.2.12) 相应的 Mathematica V4.0 计算过程为 MATHEMATICA V4.0 In[3]:= [ ] r r E λ π λ ψ λ − = 3 / 2 _, _ : （* 定义“试验”波函数 *） In[4]:= [ ] [ ] [ ] [ , ],{ ,1} 2 _, _ : [ , ],{ ,2} D r r r g r λ = Dψ r λ r + ψ λ （* 轨道角动量为零的有效拉普拉斯算符 *） In[5]:= 4 r [r, ]g[r, ]dr 0 2 π ψ λ λ ∫ ∞