ex1.m

重点内容是计算中的误差、函数方程求根、插值呆逼近、数值积分和微分、线性代数方程组解法、常微分方程初值问题的数值解法、数学软件。

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)
A\b