%%詳細推導求解過程請移步?微信公眾號?CFD入門指南?%%

clc
clear
L=0.2;?????%桿長
t=10;??????%計算時長
T1=20;?????%桿初始溫度
T2=80;?????%環境溫度
rou=8900;???%密度
c=390;????%比熱容
lan=377;????%導熱系數
h=10000;???%與環境的換熱系數
dx=0.001;???%空間步長
dt=0.0001;??%時間步長
m=t/dt;????%時間單元數量
n=L/dx;????%空間單元數量
a=lan/rou/c;????
Fo=a*dt/dx/dx;???%傅里葉數
Bi=h*dx/lan;?????%畢渥數

v1=ones（m+1n+1）*T1;?????%初始時桿上均勻溫度為T1

if?Fo>0.5?&&?（1-2*Fo*Bi-2*Fo<0）?
????fprintf（‘無法收斂，調整步長‘）;????%判斷是否滿足收斂條件
else
????for?i=1:m
????????for?j=2:n
?????????????v1（i+1j）=