In reply to Ian Wakeling.

Yes, that's another approach. It uses the "row=col" trick that I wrote about in the article "Defining Banded Matrices." I discuss the efficiency of that method in the article.