function?[ynH]=gauseidel（Abx0epsN）
%Gauss-Serdel迭代法的MATLAB函數文件gauseidel.m如下：
%A=[7?2?1?-2;9?15?3?-2;-2?-2?11?5;1?3?2?13];
%b=[4?7?-1?0]‘;x0=[0000]‘;
%[xn]=gauseidel（Ab[000，0]‘1.0e-6）
format?long?g;
if?nargin==5
????digits（abs（log10（eps）））;
elseif?nargin==4
????digits（abs（log10（eps）））;
????N=100;
elseif?nargin==3
????digits（6）;
????eps=1.0e-6;
????N=100;
elseif?nargin<3
????error
????return
end??????
D=diag（diag（A））;????%求A的對角矩陣
L=-tril（A-1）;??????%求A的下三角陣
U=-triu（A1）;???????%求A的上三角陣
G=（D-L）\U;
f=（D-L）\b;
y=G*x0+f;
n=1;??????????????????%迭代次數
while?norm（y-x0）>=eps&n????x0=y;
????y=G*x0+f;
????n=n+1;
end
y=vpa（y）;
if?n????H=‘迭代收斂‘;
else
????H=‘迭代不收斂‘;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????727??2012-06-02?12:01??楊勇—設計—雅克比\gauseidel.m

?????文件????????685??2012-06-02?11:44??楊勇—設計—雅克比\jacobi.m

?????文件???????3896??2012-06-02?11:46??楊勇—設計—雅克比\untitled.fig

?????文件??????13812??2012-06-02?11:47??楊勇—設計—雅克比\untitled.m

?????目錄??????????0??2012-06-17?20:40??楊勇—設計—雅克比

-----------?---------??----------?-----??----

????????????????19120????????????????????5