用Python作图法求解薛定谔方程

我想画出有限势阱的薛定谔方程的解。这是一个超越方程 f（E）=2 cot（2mE）1/2L/h）+（U-E）/E）1/2-（E/（U-E））1/2

其中E是变量

我的代码给了我一个参差不齐的图形

import pylab as pi
import pylab as p
import numpy as np


L = 10e-9 # meters
U = 1.602e-20 # joules
m = 10e-9 # kilograms
h = 1.055e-34 # meters squared * kilograms / seconds


x = p.linspace(0, 1.601999e-20, 1000) # before negative in square root
y = 2/p.tan(2*p.pi*((2*m*x)**(1/2))*(L/h)) + ((U-x)/x)**1/2 - (x/(U-x))**1/2


tol = 100
y[y > tol] = np.nan
y[y < -tol] = np.nan


p.plot(x, y, 'g-', lw=1)
p.show()

将pylab导入为pi
将pylab作为p导入
将numpy作为np导入
L=10e-9米
U=1.602e-20焦耳
m=10e-9千克
h=1.055e-34平方米*千克/秒
x=p.linspace（0，1.601999e-20，1000）#在平方根为负数之前
y=2/p.tan（2*p.pi*（（2*m*x）**（1/2））*（L/h））+（（U-x）/x）**1/2-（x/（U-x））**1/2
tol=100
y[y>tol]=np.nan
y[y<-tol]=np.nan
p、 绘图（x，y，'g-'，lw=1）
p、 show（）

---

有人能帮忙找出问题吗

您是否使用

Python

语言？你应该给它贴上

Python

的标签来吸引更多的注意力。对此，我很抱歉。是的，我知道了，我给你贴上了

python

标签。我可能会看到你有时有

1/2

没有括号，有时有括号。你应该注意这一点……为什么要导入两次pylab呢？它应该看起来像一个余切图。即使是平方根。