[問題] 請問迭代的問題已回收
以下是我寫的程式,我要算出b^(j+1)=A*b^j
請問我怎麼改寫呢?
clear all
h=input('h=');
k=input('k=');
N=(1/h)+1;
lambda=k/h^2;
A=eye(N,N)*(1-2*lambda);
for i=1:N-1
A(i,i+1)=lambda;
A(i+1,i)=lambda;
end
b=[];
for i=2:N
r(i)=(i-1)*(0.1)*pi;
b(i,1)=sin(r(i));
end
b(1,1)=0;
A
b
A*b
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.128.17.211