0t=42mR0-z 2e-r25/2 er27 Mathematica表达式Dpsi[r,1 ambda],r,n]]功能为,以r为变量对r,求n次偏导。幂 指数势函数()=ar的期望值为 Jd'rYGr, a)()Y(, )=" dr=24m r(n+3) (》=Jdr(,xy(Pv,4) (n+3) (8.2.13) (2A) 量纲分析要求 即耦合常数a的量纲为 则方程(82.13)的量纲也是正确的,即 E 对应的 Mathematica v4.0的指令为: MATHEMATICA V4. 0 积分*) 0ut6=∥Ra1>0&&Re>-32))am3+n12+"1 结合方程(8.2.12)和(8.2.13),可得能量表示为 E(2) 2 3) (8.2.14 对于任何λ>0的值,这一能量解始终是能量值“真解”的上界:E。≤E()。通过求可使E(a) 取最小值的变分参数λ值,即解出λ-就可以很容易地改进这一能量上限,使其逼近真解 显然,L由下式给出 对公式(8.2.14)求偏导,可得 r(n+3)Out[5]= ] 2 , 4 [Re[ ] 0, 0 5 / 2 7 / 2 3 / 2 2 dr e r e If e r r r r ∫ ∞ − − −         > − − + π λ π λ π λ π λ λ λ λ λ 4 Mathematica 表达式 D[psi[r,lambda],{r,n}]功能为，以 r 为变量对Ψ( )求 n 次偏导。幂 指数势函数V 的期望值为 r,λ ( ) n r = ar ( ) ( ) ( ) ( ) ( ) 3 3 0 2 2 3 3 * 2 3 , , 4 4 + ∞ + − Γ + Ψ Ψ = = ∫ ∫ n n r n d x r V r r a r e dr a λ π π λ π π λ λ λ λ . 即 ( ) ( ) ( ) ( ) ( ) ( ) 3 3 * 3 2 3 , , 4 + Γ + ≡ Ψ Ψ = ∫ n n V r d x r V r r a λ λ λ λ . (8.2.13) 量纲分析要求 V → n = ar n aE− → E 即耦合常数 a 的量纲为 a → n+1 E . 则方程(8.2.13)的量纲也是正确的, 即 n a − λ → E E E n n = +1 − . 对应的 Mathematica V4.0 的指令为： ------------------------------------------------------------------------------ MATHEMATICA V4.0 ------------------------------------------------------------------------------ In[6]:= r E dr (* 积分 *) n 2λr 0 2 − ∞ + ∫ Out[6]= [Re[ ] 0& & Re[ ] 3,2 [3 ], ] 2 0 3 3 2 If n Gamma n e r dr n n r +n ∞ − − − − − ∫ > > − + λ λ λ ------------------------------------------------------------------------------ 结合方程(8.2.12)和(8.2.13)，可得能量表示为 ( ) ( ( ) ) n a n E µ λ λ λ 2 3 2 2 2 Γ + = + . （8.2.14） 对于任何λ > 0 的值，这一能量解始终是能量值“真解”的上界：Etrue ≤ E(λ)。通过求可使E(λ) 取最小值的变分参数λ 值，即解出λ min，就可以很容易地改进这一能量上限，使其逼近真解。 显然，λ min由下式给出 ( ) = 0 ∂ ∂ λ E λ . 对公式（8.2.14）求偏导，可得 ( ) ( ) ( ) 0 2 3 1 = Γ + = − ∂ ∂ n+ n an E µ λ λ λ λ