localized-state ensemble (LSE) model

how can i solve the equation in the localized-state ensemble (LSE) model by MATLAB software ?

how can i do the resolution of parameter x?

How can i fetched parameters from formula by Matlab ?

E(T) and T are known .

E1= E0 – (alfa* T^2)/(theta+T)

E2= E0 – (alfa* T^2) / (theta+T) – x * KB * T

(x * e^x)= [(sigma /(KB * T))^2 – X ] * (τaur/τautr ) * e^(((E0 – Ea)) ⁄ KB * T)

where E1 and E2 are the energy , and they are related to each other the diffrent is in E2 we have the part of x incloud ,, and x change each time when T change , it is a small change but it apear in the figer

T is the temperature and it is change from 0 to 300 ,, x is the temperature-dependent coefficient and i have a condition on x which is : 0<X < (sigma/(KB.*T(i)))^2

that x is larger than zereo and less than (sigma/(KB*T(i)))^2