ex2.asv

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

n=input('输入矩阵阶数：');
a=zeros(n,1);b=zeros(n,1);c=zeros(n-1,1);d=zeros(n,1);
for i=1:n-1;
    a(1)=2;a(i+1)=2;
    b(i)=1;c(i)=1;
    d(1)=-7;d(i+1)=-5;
end
for j=2:n;
    q(1)=b(1);
    p(j)=a(j)/q(j-1);
    q(j)=b(j)-p(j)*c(j-1);
    y(1)=d(1);
    y(j)=d(j)-p(j)*y(j-1);
end
x(n)=y(n)/q(n);
for k=n-1;-1:1;
    x(k)=(y(k)-c(k)*x(k+1))/q(k);
end
x