Calculated alpha values for the flare frequency distribution within the energy range of 10³³ to 10³⁷ erg using the MCMC method. The green error bars represent the statistical uncertainty, derived as the square root of the total number of flares within each energy range. The dashed vertical red line represents that the fitting energy starts from 10³³ ergs. The grey lines represent predictive samples from posterior distributions, while the red line represents the mean trend determined by calculating the average fitted slope and its uncertainty. For the entire sample of flaring stars, the alpha value is found to be approximately 1.22.