The bandmatrix administration comes from a Direct Solver
which is the worst thing you can implement for this problem !!

Next sheet ==>