ex2_2.m

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

n=input('输入矩阵阶数：');
t=cputime;
a=zeros(n-1,1);b=zeros(n,1);c=zeros(n-1,1);d=zeros(n,1);
a=linspace(1,1,n-1);b=linspace(2,2,n);c=linspace(1,1,n-1);
for i=1:n;
d(i)=-5;d(1)=-7;
end
A=diag(a,-1)+diag(b)+diag(c,1);
[L,U]=lu(A);
y=L\d;
x=U\y;
e=cputime-t