renamed getNN to getNeighbors
plot dispersion relation along 1D path, plot 2D dispersion relation
Use new functions for getting (next) nearest neighbors
Correctly use solve() on masked arrays (vectors outside the Brillouin zone)