How to plot energy spectra for SSH model in presence of on-site potential?