function Xi=Xi_cal(t,E,Xo)
T = t+273.15; % K
R = 8.314; % unit J/mol/K
F = 96.48534; % C/mol
Xi=Xo/exp(E*F/(R*T));