It looks like the method you used has a lot in common with a finite elements analysis. At least the simple version of it that structural engineers use in software like sap2000 and RISA. Do you think it would work if the matrix in the equation were a lot bigger?