//@code objref m, v1, minv m = new Matrix(7,7) v1 = new Vector(7) for i=-1,1 { v1.fill(2 - 3*abs(i)) m.setdiag(i, v1) } minv = m.inverse() m.printf minv.printf m.mulm(minv).printf