A=[1010;0120;1210;0001];
s=10;??%循環次數
d=zeros（s+11）;
d（1）=dot（diag（A）diag（A））;
[mn]=size（A）;
for?k=1:s???%?采用Jacobi法迭代的次數
????a_max=0;??%初設A的各元素中a12最大。
????for?i=1:m
????????for?j=1:n
????????????if?abs（A（ij））>a_max&&i~=j????%?循環比較，找出最大元素
????????????????a_max=A（ij）;
????????????????i_max=i;
????????????????j_max=j;
????????????end
????????end
????end
????t1=（A（i_ma