The study aims to develop a reliable model using gene-expression programming (GEP) technique for estimating the unconfined compressive strength (UCS) of soil stabilization by cement and fly ash. The model considered the effects of several parameters, including the fly ash characteristics such as calcium oxide (CaO) content, CaO/SiO2 ratio, and loss of ignition. The research results show that the proposed model demonstrates superior performance with a high correlation coefficient (R > 0.955) and low errors. Therefore, the model could be confidently applied in practice for a variety of fly ash qualities. Besides, the parametric study was conducted to examine the effect of fly ash characteristics on the strength of soil stabilization. The study indicates that if the fly ash contains a high amount of calcium oxide, the strength of fly ash stabilized soil is significant. In addition, fly ash could be used in combination with cement to increase the strength of the mixture. A fly ash replacement ratio is suggested from 0.19 to 0.35, corresponding to the total binder used from 10% to 30%. The research findings could help engineers in optimizing the fly ash proportion and estimating the UCS of soil stabilization by cement and fly ash.