试试下面的代码,有问题可追问
A=zeros(9,9);
b=zeros(9,1);
for i=1:9
if i-1>=1
A(i,i-1)=0.85;
end
A(i,i)=-1.98;
if i+1<=9
A(i,i+1)=1.15;
end
b(i)=1+0.1*i;
end
y0=1;y10=6;
b(1)=b(1)-0.85*y0;
b(9)=b(9)-1.15*y10;
y=inv(A)*b
后一个差分方程中,计算 y(i+1) 需要用到其前两步的值,即 y(i) 和 y(i-1),你只给了y(0),好像条件不足吧?
另外,y(10)=6 是什么意思?要求最后计算出来的结果 i=9 时,y(10) 要等于6吗?