function?[yn]=gauseseidel（Abx0ep）
if?nargin==3
????ep=1.0e-6;
elseif?nargin<3
????error
????return
end
D=diag（diag（A））;????????%求A的對角矩陣
L=-tril（A-1）;??????????%求A的下三角矩陣
U=-triu（A1）;?????