?%test?計算PE（0）
Fb0=0.4412;
P0=0.897;
T0=650;?%初始溫度
beta=2000;?%升溫速率
%不穩橋在全部鏈接中的概率
PB=Fb0*P0;
%計算反應物中不穩定橋的概率
w=[P0Fb0];?%w（1）---p（0）--完整鏈接在全煤中的分數???w（2）--Fb（0）--不穩定橋在完整鏈接中所占比例
Prb0=PRb（w）;
ww=[Prb0P0];
BR_0=BR（ww）;
%計算末端含有官能團的概率
%x=[275101Fb0*P01.8560.9110.8387]?%x（1）--MWmon-單體單元分子量???x（2）--MWA--核的分子量??x（3）--pb--不穩定橋占完整鏈接比例??x（4）--MWB/MWA???x（5）--p????x（6）--MWC/MWA%
%Pe0=pe（x）
Pe0=1;
%計算反應物中的外圍結構
y=[Pe0P0];?%y（1）--pe（0）??y（2）--p（0）
SR0=SR（y）;?
%計算后形質體的摩爾數

m5_0=m5（P0）;
m4_0=m4（P0）;
m3_0=m3（P0）;
m2_0=m2（P0）;
m1_0=（1-P0）^2;

%計算后形質體的不穩定橋

z=[Fb0P0];
BM5_0=BM5（z）;
BM4_0=BM4（z）;
BM3_0=BM3（z）;
BM2_0=BM2（z）;
BM1_0=0;

%計算SM
SM0=SM（y）;

%計算中間體
n10_0=n10（P0）;
n9_0=n9（P0）;
n8_0=n8（P0）;
n7_0=n7（P0）;
n6_0=n6（P0）;

%計算中間體的不穩定橋
BN10_0=BN10（z）;
BN9_0=BN9（z）;
BN8_0=BN8（z）;
BN7_0=BN7（z）;
BN6_0=BN6（z）;
%計算中間體的外圍結構
SN_0=SN（y）;

%初始非冷凝氣體?表III??fuel?and?energy?55-1991663
YG_0=0;
%計算KB不穩定橋斷裂速率常數的表達式及值?在T溫度的值
t=[0:0.014:0.275];
%T=t*beta+T0;
%kB=kB_splinet（T）;
%外圍結構脫除反應速率
%kG=1d15*exp（（-230000）*power（8.314*T-1））;?%kG--氣體逸出速率常數表達式