function?onedicon

%?—————————和一維非穩(wěn)態(tài)擴(kuò)散算例一樣的初始和邊界條件?p122頁碼
%?—————————采用乘方格式?時(shí)間步長(zhǎng)為0.001s
%?—————————初始溫度場(chǎng)為200，速度為2m/s，長(zhǎng)度為2cm，t=0s時(shí)刻?hào)|側(cè)溫度突然降至0C。
%?—————————時(shí)間差分采用全隱式格式
clc
disp（‘一維無內(nèi)熱源非穩(wěn)態(tài)全隱?乘方格式離散，初始場(chǎng)為200C，右側(cè)邊界溫度變?yōu)?，左側(cè)邊界條件可選擇‘）;
disp（‘----------------------------------------------------------------------------‘）
a=input（‘選擇左邊界條件：1—第一類邊界條件，2--第二類邊界條件?:???‘）;
q=0;?t1=0;
if（a==0）
????q=input（‘熱流為?：?‘）;
else
????t1=input（‘選擇左邊界第一類條件溫度，取值在0-200內(nèi)：?‘）;
end
t2=input（‘選擇右邊界第一類條件溫度，取值在0-200內(nèi)：??‘）;
u=input（‘please?input?the?velocity?:???‘）;
n=input（‘長(zhǎng)度2cm上劃分網(wǎng)格數(shù)?:???‘）;
dt=input（‘時(shí)間步長(zhǎng)?:???‘）;
%————————————幾何條件，物性參數(shù)
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）;
%————————————沒有內(nèi)熱源
sp=0;
%————————————初始溫度場(chǎng)
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