Jacobi法求特征值的matlab程序

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