第五章.txt

来自「包括牛顿法」· 文本 代码 · 共 22 行

b=2*ones(100,1);
c=-1*ones(100,1);
a=-1*ones(100,1);
d=zeros(100,1);
d(1)=1;
d(100)=1;
u=zeros(100,1);
l=zeros(100,1);
y=zeros(100,1);
l(1)=b(1);
y(1)=d(1)/l(1);
u(1)=c(1)/l(1);
for i=2:100
l(i)=b(i)-a(i)*u(i-1);
y(i)=(d(i)-y(i-1)*a(i))/l(i);
u(i)=c(i)/l(i);
end
x=zeros(100,1);
x(100)=y(100);
for i=99:-1:1
x(i)=y(i)-u(i)*x(i+1);
end