％参数
mliulv=20/3600;
Cp=4200;
Di=0.007;
hfi=240;
houdub=0.0005;
daoreb=3.12;
houdut=0.0005;
daoret=398;
 
m=6;
length=2;
n=1000;
 
Tp=[297.7931135 297.773487  297.7575334 297.7488286 297.7457609 297.7467345 297.7506569 297.7568042 297.7646826 297.7739414 297.7843217 297.7956264 297.8077009 297.820422  297.8336888 297.8474185 297.8615417 297.876 297.8907439 297.9057313 297.9209263 297.9362978 297.9518194 297.9674684 297.9832254 297.9990735 298.0149984 298.030988  298.0470316 298.0631205 298.0792471 298.0954049 298.1115885 298.1277935 298.1440159 298.1602527 298.176501  298.1927588 298.2090242 298.2252956 298.241572  298.2578522 298.2741355 298.2904211 298.3067086 298.3229976 298.3392876 298.3555785 298.3718699 298.3881618 298.404454  298.4207464 298.4370389 298.4533315 298.4696241 298.4859166 298.5022091 298.5185015 298.5347939 298.5510861 298.5673782 298.5836702 298.599962  298.6162538 298.6325454 298.6488369 298.6651282 298.6814194 298.6977105 298.7140014 298.7302922 298.7465829 298.7628734 298.7791638 298.795454  298.8117442 298.8280341 298.844324  298.8606137 298.8769032 298.8931927 298.9094819 298.9257711 298.9420601 298.958349  298.9746377 298.9909263 299.0072148 299.0235031 299.0397913 299.0560793 299.0723672 299.088655  299.1049426 299.1212301 299.1375175 299.1538047 299.1700918 299.1863787 299.2026655 299.2189522 299.2352387 299.2515251 299.2678114 299.2840975 299.3003835 299.3166693 299.332955  299.3492406 299.365526  299.3818113 299.3980964 299.4143814 299.4306663 299.446951  299.4632356 299.4795201 299.4958044 299.5120886 299.5283727 299.5446566 299.5609403 299.577224  299.5935075 299.6097908 299.626074  299.6423571 299.65864   299.6749229 299.6912055 299.707488  299.7237704 299.7400527 299.7563348 299.7726167 299.7888986 299.8051803 299.8214618 299.8377432 299.8540244 299.8703055 299.8865864 299.9028672 299.9191478 299.9354281 299.9517083 299.9679882 299.9842679 300.0005472 300.0168261 300.0331046 300.0493825 300.0656597 300.0819361 300.0982114 300.1144854 300.1307577 300.147028  300.1632957 300.1795601 300.1958206 300.212076  300.2283252 300.2445667 300.2607988 300.2770193 300.2932255 300.3094145 300.3255825 300.341725  300.3578368 300.3739117 300.3899424 300.4059204 300.4218355 300.4376761 300.4534284 300.4690766 300.4846022 300.4999838 300.5151965 300.5302117 300.5449959 300.5595105 300.5737101 300.5875421 300.6009443 300.6138432 300.6261509 300.6377609 300.6485425 300.6583317 300.6669189 300.6740274 300.6792803 300.6821437 300.6818295 300.6771415 300.6663272 300.6476346 300.6252002];
a=0.54945/(955.18*Cp);
dy=0.01;
dt=0.1;
Fo=a*dt/dy^2;
 
  for j=2:200
         Tf(j,k+1)=0;
         T2=5;
    while abs(T2-Tf(j,k+1))>=0.001
 Tf(j,k+1)=T2;
qv=(Tp(j)-Tf(j,k))/(1/hfi/Di/pi+houdub/daoreb/dy+houdut/daoret/dy);
               
T2=(1-2*Fo)*Tf(j,k)+Fo*(Tf(j-1,k)+Tf(j+1,k))+dt*qv/(955.18*Cp); 
 end
  end
end
for j=1:200
plot(j,Tf(j,k)),xlabel('direction y'),ylabel('temperature')
hold on
end


大家帮忙看一下程序，这是计算一维非稳态流体温度梯度的程序，计算出来的温度是一样的，与初值相同
帮忙分析一下是怎么回事？
谢谢大家！