from graphenemodeling.graphene import monolayer as mlg
from graphenemodeling.graphene import _constants as _c
import matplotlib.pyplot as plt
E = np.linspace(-3,3,num=200) * _c.g0
DOS_low = mlg.DensityOfStates(E,model='LowEnergy')
DOS_full = mlg.DensityOfStates(E,model='FullTightBinding')
plt.plot(E/_c.g0,DOS_full/np.max(DOS_full),'k-',label='FullTightBinding')
# [<...
plt.plot(E/_c.g0,DOS_low/np.max(DOS_full),'k-.',label='LowEnergy')
# [<...
plt.xlabel('$E/\gamma_0$')
# Text...
plt.ylabel('DOS (a.u.)')
# Text...
plt.legend()
# <...
plt.show()