minor changes
added plot_3D_spectrum
Hamiltonian H now takes 3 arguments; minor changes