ghmatrix.m

来自「jaccobi迭代法和GS迭代法求解线性方程组」· M 代码 · 共 11 行

%Hilbert矩阵生成
%著名的病态正定对称阵
function [h,b]=gHmatrix(n)
h=zeros(n);b(n)=0;
for i=1:n
    for j=1:n
        h(i,j)=1/(i+j-1);
        b(i)=b(i)+j/(i+j-1);
    end
end