jacobiexample.m

来自「数值线性代数的Matlab应用程序包 共13个程序函数」· M 代码 · 共 12 行

A=rand(3,3);
A=triu(A);
[m,n]=size(A);
for p=1:m
    for q=p+1:n
        A(q,p)=A(p,q);
    end
end
B=Jacobi(A);
A
B