R=matrix( [[1,-0.25,-0.25,0],[-0.25,1,0,-0.25],[-0.25,0,1,-0.25],[0,-0.25,-0.25,1] ])
print R
p=vector( [50,50,25,25])
x=R.solve_right(p)
print x

}}}

Kreyszig-18.3 1 (last edited 2010-12-18 05:15:15 by 10)