Python中径向波函数的解析解图？

Question

我正在尝试在Python中绘制以下等式的图。

二维量子环的径向微分方程的解

beta参数是

<< img src =“ https://image.soinside.com/eyJ1cmwiOiAiaHR0cHM6Ly9pLnN0YWNrLmltZ3VyLmNvbS9aaHBubS5wbmcifQ==” alt =“在此处输入图像描述”>“ >>

这是我的尝试

import numpy as np
from scipy.special import gamma, eval_genlaguerre
from scipy import exp
import matplotlib.pyplot as plt


#Dados experimentais
m= 0.067*511
R=5
r=np.linspace(0, 40)
ro = r/R
flux = np.linspace(0,1.0)


def beta(gama): 
   return np.sqrt((m-flux)**2+(gama**4)/4)


def amplitude(gama):
    return (np.sqrt(gamma(1)/(2**beta(gama)*gamma(beta(gama)+1))))


def Rn(gama):
    eval_g = eval_genlaguerre(0,beta(gama),(((gama * ro)**2)/2))
    exp_g = exp(-((gama * ro)**2)/4)
    return amplitude(gama) * (gama*ro)**beta(gama)*exp_g*eval_g

print(beta(1.5)) 
print(amplitude(1.5))
print(Rn(1.5))


sol1 = Rn(1.5)
sol2 = Rn(2.0)
sol3 = Rn(2.5)
sol4 = Rn(3.0)


fig, ax = plt.subplots()
ax.plot(ro, sol1, color = 'red', label = '$\gamma$ = 1.5')
ax.plot(ro, sol2, color = 'green', label = '$\gamma$ = 2.0')
ax.plot(ro, sol3, color = 'blue', label = '$\gamma$ = 2.5')
ax.plot(ro, sol4, color = 'black', label = '$\gamma$ = 3.0')
ax.legend()
ax.set_xlabel('R/r')
ax.set_ylabel('$R_0(r)$')
该代码的输出甚至与我应该重现的情节略有不同，但是没有语法错误。我认为我没有正确使用Python函数。如何给不同的gamma值并将其传递给beta参数和径向波函数？

我正在尝试在Python中绘制以下方程式的图形。 2d量子环径向微分方程的解beta参数是这是我的尝试，将numpy导入为np ...

Answer 1

我不是这个话题，但是狂热中至少存在以下错误（我故意设置为n = 0时是无声的）：