In this work, the solubility of bezafibrate in 16 pure solvents including 9 alcohols (methanol, ethanol, n-propanol, iso-propanol, n-butanol, iso-butanol, sec-butanol, 1-pentanol, and iso-pentanol), 4 alkyl acetate (ethyl acetate, n-propyl acetate, n-butyl acetate, and methyl propionate), 2 ketones (methyl ethyl ketone and cyclohexanone) and one nitrile (acetonitrile) and in binary mixed solvents (methanol + acetonitrile) at temperature of 283.15-323.15 K was determined. The solubility of bezafibrate in pure solvents was simulated using modified Apelblat, Van't Hoff, Wilson, and NRTL model. That in binary mixed solvents was simulated using the modified Apelblat and Apelblat-Jouyban-Acree model. The thermodynamic mixing properties (Delta(mix)G, Delta H-mix, Delta S-mix, gamma(infinity)(1), and H-1(E),infinity) of the bezafibrate-pure solvent systems were evaluated based on the solubility data and Wilson model.