a=matrix([[5,0,0,0],[0,5,0,0],[0,0,5,0],[0,0,0,5]])
b=matrix([[1,1,1,1],[1,1,1,1],[1,1,1,1],[1,1,1,1]])
b
a
a-b
c = (a-b)**-1
c #a=2/5,b=1/5;

Strang-ex-1.6-32-U (last edited 2010-12-17 12:15:46 by prasad)