%雅克比迭代法：建立函數文件Jacobi.m
function?Jacobi（AbXnmw_error）
%用雅克比迭代法求解方程組Ax=b
%輸入：A為方程組的系數矩陣，b為方程組右端的列向量，X為迭代初值構成的列向量，nm為最大迭代次數，w_error為誤差精度
%輸出：x為求得的方程組的解構成的列向量，n為迭代次數
n=1;
m=length（A）;
D=diag（diag（A））;????????%令A=D-L-U計算矩陣D
L=tril（-A）+D;???????????%令A=D-L-U計算矩陣L
U=triu（-A）+D;??????????%令A=D-L-U計算矩陣U
M=inv（D）*（L+U）;???????%計算迭代矩陣
g=inv（D）*b;????????????%計算迭代格式中的常數項
%