>
|
with(LinearAlgebra):
M := Matrix(4,4,[[5.0],[0.2,-1.0],[0.0,3.0,0.0],[-2.0,0.0,0.0,1.0]],
shape=symmetric,storage=sparse[lower],datatype=float):
b := Vector(4,1.0,datatype=float):
LinearSolve(M,b,method=SparseIterative,methodoptions=[itermethod=CG,levelfill=1]);
|