ex1_1.m

来自「重点内容是计算中的误差、函数方程求根、插值呆逼近、数值积分和微分、线性代数方程组」· M 代码 · 共 13 行

n=input('输入矩阵阶数：');
A=zeros(n,n);b=zeros(n,1);
for i=1:n;
    for j=1:n;
        A(i,j)=(i+j-1)^2;
    end
    b(i)=sum(A(i,:));
end
x=GaussSolve(A,b)
[L,U]=lu(A);
y=L\b;
z=U\y