5
I think that the main problem might be with the solver you are using. The Hamiltonian (matrix) in this case is Hermitian, it is even symmetric since it is purely real. You could use eigh instead of eig to take advantage of this.
Furthermore, you are not removing only the first and last points but intervals of size 1 at each end. Following, I show you a ...
Only top voted, non community-wiki answers of a minimum length are eligible
