mag.m

来自「信号与系统分析以matlab实现」· M 代码 · 共 29 行

function f=mag(n)
if n>2
    x=magic(n)
    for j=1:n
        s1=0;
        for i=1:n
            s1=s1+x(j,i);
        end
        s1
    end
    for i=1:n
        s2=0;
        for j=1:n
            s2=s2+x(j,i);
        end
        s2
    end
    s3=0;
    for i=1:n
        for j=1:n
            if (i==j)
                s3=s3+x(j,i);
            end
        end
    end
    s3
else
    break
end