Inductance Matrix for square array

Calculates the full inductance matrix for an array of interconnected square cells

I wrote this to calculate the full inductance matrix (L) for a rectangular array of square cells which are electrically connected with non-zero width leads. I realize I can do this in fasthenry, but I wanted to do this as part of a larger simulation in which the array dimensions are user-specified, so I need to calculate L for different size arrays. Because of the fairly simple geometry, I felt I could use an analytic solution instead of using an FEA approach to solve for L. Any comments, suggestions or questions are welcom.