0ut[16]={1.96556,{A->1.14471} 将n=a=1代入公式(8.2.15),(8.2.16)和(8.2.17),可对以上结果进行检验。通常将 计算结果绘图显示,有助于进行比较。令2=1,我们分别绘制函数E=2.3812和 E=2.4764a23,然后借助 Mathematica v3.0的Show指令将两图形合在一起进行比较。 MATHEMATICA V3. 0 定义函数 In[17]:= etrue[a]:=2.3381a^(2/3) In[18]:= upper[a]:=2.4764a^(2/3) In[19]: =plot=Plot [etrue[a], fa, 0, 3), AxesLabe1->("a","E", Textstyle-> Fontslant-> Italic”, Fontsize->14}](*绘制E(plot1)*) (* AxesLabe1:定义坐标轴的表述; Textstyle定义字型和字体大小。* 0.511.522.53 图82.2线性势基态能量真值E In [20]: =plot2=Plot [eupper[a],(a, 0, 31, AxesLabel->("a""E"I, TextStyle-> FontSlant ->"Italic", FontSize->14] (*绘制E(plot2)*) 0.511.5 图823线性势基态的变分上限能量EOut[16]= {1.96556,{ λ -> 1.14471}} 将 代入公式(8.2.15),(8.2.16)和(8.2.17)，可对以上结果进行检验。通常将 计算结果绘图显示，有助于进行比较。令 n = a = 1 2µ = 1，我们分别绘制函数 和 ，然后借助 Mathematica V3.0 的 Show 指令将两图形合在一起进行比较。 2 / 3 E 2.3381a true = 2 / 3 2.4764a Evar = MATHEMATICA V3.0 定义函数 In[17]:= etrue[a_] := 2.3381 a^(2/3) In[18]:= eupper[a_] := 2.4764 a^(2/3) In[19]:=plot1=Plot[etrue[a],{a,0,3}, AxesLabel->{”a”,”E”}, TextStyle-> {FontSlant->”Italic”,FontSize->14}] (* 绘制 (plot1) *) Etrue （* AxesLabel: 定义坐标轴的表述； TextStyle 定义字型和字体大小。*） 图8.2.2 线性势基态能量真值 。 In[20]:=plot2=Plot[eupper[a],{a,0,3},AxesLabel->{”a”,”E”},TextStyle-> 2) *) 图8.2.3 线性势基态的变分上限能量 。 true E {FontSlant ->”Italic”,FontSize->14}] (* 绘制 Evar (plot Evar