以下函数用fjacobi.m文件名存盘function X=fjacobi(A,b,X0)D=diag(diag(A));% tril(A)求矩阵的左下三角,tril(A,-1)就不包括主对角线。L=-tril(A,-1);U=-triu(A,1);B=【1】;F=D\b;X=B*X0+F;n=1;m=30;while norm(X-X0)>;=1e-5 && n<;=mX0=【2】;X=【3】;n=n+1;end在命令窗口输入a,b,x0,调用fjacobi函数得到结果>;>; a=[10 -2 -1;-2 10 -1;-1 -2 5];>;>; b=[3;15;10];>;>; x0=[0 0 0]';>;>; x=【4】x =1.00002.00003.0000
以下函数用fjacobi.m文件名存盘function X=fjacobi(A,b,X0)D=diag(diag(A));% tril(A)求矩阵的左下三角,tril(A,-1)就不包括主对角线。L=-tril(A,-1);U=-triu(A,1);B=【1】;F=D\b;X=B*X0+F;n=1;m=30;while norm(X-X0)>;=1e-5 && n<;=mX0=【2】;X=【3】;n=n+1;end在命令窗口输入a,b,x0,调用fjacobi函数得到结果>;>; a=[10 -2 -1;-2 10 -1;-1 -2 5];>;>; b=[3;15;10];>;>; x0=[0 0 0]';>;>; x=【4】x =1.00002.00003.0000
1