why my loop doesn't run or produce an error at all?

Question

Xi Luan 2022-9-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1806905-why-my-loop-doesn-t-run-or-produce-an-error-at-all

评论： Xi Luan 2022-9-19

采纳的回答： Abolfazl Chaman Motlagh

cranknicholsen.m

for NT = 2:60 %number of time steps

for n = 533:475 %number of grid points

maindiag = 2*(1+r)*ones(n+1,1);%n+1 because including one ghost point below i = 0

lowdiag = -r*ones(n+1,1);

updiag = -r*ones(n+1,1);

S = spdiags([lowdiag maindiag updiag],-1:1,n+1,n+1);

S = full(S);

S(1,1) = 1;

S(1,2) = -0.007875;

S(1,3) = -1;

S(n+1,n-1) = -1;

S(n+1,n) = 0;

S(n+1,n+1) = 1;

b(1,1) = 0;%first row of RHS vector b

%interior rows of RHS vector b

for i = 1:n-1

b(i+1,1) = r*C(i,1)+(2-2*r)*C(i+1,1)+r*C(i+2,1);

end

b(n+1,1) = 0;%final row of RHS vector b

C(1:n+1)= linsolve(S,b);

plot((0:dz:totz)',C)

end

I'm trying to use crank nicholson to solve for fick's 2nd law diffusion equation. The above code is a section of it. The fourth line from bottom is how I try to save the calculated result/concentration in each time step. r is a constant. Initial C matrix has been calculated and is a column vector at the first time step (NT = 1).

After run, n = [] and NT = 60. but C matrix didn't update at all and remained as initial at NT = 1.