高斯-塞德?tīng)柕ǖ牡仃囋趺辞?？在Jacobi迭代中，a=d-l-u，ax=B，DX=（Lu）x，BX=d^{-1}（Lu）x，d^{-1}B所以B=d^{-1}（Lu）在Gauss-Seidel

高斯-塞德?tīng)柕ǖ牡仃囋趺辞螅?/h2>

在Jacobi迭代中，a=d-l-u，ax=B，DX=（Lu）x，BX=d^{-1}（Lu）x，d^{-1}B

所以B=d^{-1}（Lu）

在Gauss-Seidel迭代中，a=d-l-u，但是ax=B（d-l）x=UX，BX=（d-l）^{-1}UX（d-l）^{-1}B

所以B=（d-l）^{-1}u

%Gauss-Seidel函數(shù)%a是未知量的系數(shù)矩陣；B是未知量的系數(shù)矩陣方程的右常量列向量，x0取1（m，1）m為未知數(shù)，EPS為精度，如果不輸入EPS，則默認(rèn)為1.0e-6函數(shù)[x，n]=Gauss-Seidel（a，B，x0，EPS）如果nargin==3 EPS=1.0e-6elseif nargin<3 error returnendD=diag（diag（A））L=-tril（A，-1）U=-triu（A，1）G=（D-L）Uf=（D-L）bx=G*x0 fn=1而norm（x-x0）>=EPS x0=x x=G*x0 f n=n1