function?pde_TwoTemperature_solve

%??????Ce（Te）*DTe/Dt=D（Ke（Te）*DTe/Dx）/Dx-g*（Te-Ti）+S（xt）;?????
%??????Ci（Ti）*DTi/Dt=g*（Te-Ti）;
%
%???In?the?form?expected?by?PDEPE?the?equations?are
%
%????|Ce（Te）|?????????|Te|??????|?Ke（Te）*?DTe/Dx??|????|?-g*（Te-Ti）+S（xt）?|
%????|??????|?.*??D_??|??|?=?D_?|?????????????????|?+??|???????????????????|
%????|??Ci??|?????Dt??|Ti|???Dx?|????????0????????|????|?????g*（Te-Ti）?????|

%????????c???????????????????????????????f????????????????????????s


clear?all
format?long
m?=0;
?x=0:10e-10:400e-9;%步長為1nm
?t=0:10e-15:10e-12;%步長為10fs
%%---------------------------------

sol?=?pdepe（m@pdefun@pdeic@pdebcxt）;
T1=?sol（::1）;
T2=?sol（::2）;?

subplot（221）
mesh（xtT1）;
title（‘電子溫度Te（xt）‘）;
xlabel（‘深度x/m‘）;
ylabel（‘時間t/s‘）;
?%hidden?on
?%shading?interp

subplot（222）
mesh（xtT2）;
title（‘晶格溫度Ti（xt）‘）;
xlabel（‘深度x/m‘）;
ylabel（‘時間t/s‘）;

subplot（223）
plot（tT1（:1）‘r‘tT2（:1）‘b‘）;%T1（:10）表示取T1矩陣的第10列全部元素?即取特定x=10nm處的?溫度隨時間的變化
text（0.35e-122500‘?\leftarrow?Te‘）
text（1e-12380‘?\uparrow?Ti‘）
title（‘電子晶格溫度分布Te（0t）Ti（0t）.‘）;
xlabel（‘時間?t/s‘）;
ylabel（‘溫度?T/K‘）;

subplot（224）
plot（xT1（50:）‘r‘xT2（50:）‘b‘）;%T1（100:）表示取T2矩陣的第100行全部元素?即取特定t=1ps時?溫度隨時間的變化
title（‘電子晶格溫度分布Te（x1ps）Ti（x1ps）.‘）;
xlabel（‘深度?x/m‘）;
ylabel（