def plot_eigs(mesh, eigs): from pylab import plot, show import numpy l = Linearizer(mesh) for E, eig in eigs: if E >= 0: break x, y = l.get_xy(eig, 0, 100) print "plotting E=%f" % E plot(x, y) show()
def plot_eigs(mesh, eigs): try: from jsplot import plot, show except ImportError: from pylab import plot, show import numpy l = Linearizer(mesh) for E, eig in eigs: if E >= 0: break x, y = l.get_xy(eig, 0, 100) print "plotting E=%f" % E plot(x, y) show()
def plot_eigs(mesh, eigs, n=10): try: from jsplot import plot, show except ImportError: from pylab import plot, show import numpy l = Linearizer(mesh) for E, eig in eigs: if E < -6: continue if E >= 6: break x, y = l.get_xy(eig, 0, 100) print "plotting E=%f" % E plot(x, y) show()
def plot_eigs(mesh, eigs): #try: # from jsplot import plot, show #except ImportError: from pylab import plot, show, legend import numpy from _forms import potential_python l = Linearizer(mesh) for E, eig in eigs: if E >= 0: break x, y = l.get_xy(eig, 0, 100) print "plotting E=%f" % E plot(x, y, label="E=%f" % E) y = [potential_python(_) for _ in x] plot(x, y, "k-", label="V(r)") legend() show()
def plot_eigs(mesh, eigs, n=10): try: from jsplot import plot, show except ImportError: from pylab import plot, show import numpy l = Linearizer(mesh) count = 0 c = 137.036 n = 4 for E, eig in eigs: if E < -2*c**2: # skip the positron states continue if E > 0: break count += 1 if count > n: break x, y = l.get_xy(eig, 0, 100) print "plotting E=%f" % E plot(x, y) show()