function?onedicon

%?—————————和一維非穩態擴散算例一樣的初始和邊界條件?p122頁碼
%?—————————采用乘方格式?時間步長為0.001s
%?—————————初始溫度場為200，速度為2m/s，長度為2cm，t=0s時刻東側溫度突然降至0C。
%?—————————時間差分采用全隱式格式
clc
disp（‘一維無內熱源非穩態全隱?乘方格式離散，初始場為200C，右側邊界溫度變為0，左側邊界條件可選擇‘）;
disp（‘----------------------------------------------------------------------------‘）
a=input（‘選擇左邊界條件：1—第一類邊界條件，2--第二類邊界條件?:???‘）;
q=0;?t1=0;
if（a==0）
????q=input（‘熱流為?：?‘）;
else
????t1=input（‘選擇左邊界第一類條件溫度，取值在0-200內：?‘）;
end
t2=input（‘選擇右邊界第一類條件溫度，取值在0-200內：??‘）;
u=input（‘please?input?the?velocity?:???‘）;
n=input（‘長度2cm上劃分網格數?:???‘）;
dt=input（‘時間步長?:???‘）;
%————————————幾何條件，物性參數
L=0.02;r=1.0;k=0.03;dx=L/n;
F=r*u;D=k/dx;pe=F/D;
TT=ones（n+11）;
%————————————沒有內熱源
sp=0;
%————————————初始溫度場
tp0=200*ones（n1）;
x=linspace（00.02n+1）;
for?i=1:n
?????xx（i）=（x（i）+x（i+1））/2;
end
xxx（n+1）=0.02;
x