Abstract
We apply two mathematical techniques, specifically, the unified solver approach and the exp(-phi(xi))-expansion method, for constructing many new solitary waves, such as bright, dark, and singular soliton solutions via the fractional Biswas-Milovic (FBM) model in the sense of conformable fractional derivative. These solutions are so important for the explanation of some practical physical problems. Additionally, we study the stochastic modeling for the fractional Biswas-Milovic, where the parameter and the fraction parameters are random variables. We consider these parameters via beta distribution, so the mathematical methods that were used in this paper may be called random methods, and the exact solutions derived using these methods may be called stochastic process solutions. We also determined some statistical properties of the stochastic solutions such as the first and second moments. The proposed techniques are robust and sturdy for solving wide classes of nonlinear fractional order equations. Finally, some selected solutions are illustrated for some special values of parameters.