4-1.m

来自「matlab具有强大的数值计算能力」· M 代码 · 共 27 行

n=1; k=1;
while n<=100
    a=ones(n,n);
    for i=1:n
        for j=1:n
            a(i,j)=1/k;
            k=k+1;
            j=j+1;
        end
        i=i+1;
        k=i;
    end
    A(n)=log(cond(a,2));
    n=n+1; 
end
plot(linspace(1,100),A);=1:n
            a(i,j)=1/k;
            k=k+1;
            j=j+1;
        end
        i=i+1;
        k=i;
    end
    A(n)=log(cond(a,2));
    n=n+1; 
end
plot(linspace(1,100),A);