The solubility of ivermectin in high pressure carbon dioxide (CO2) at 338, 328, 318, and 308 K temperatures was investigated with and without an ethanol. Several methods were employed to model the data obtained from the experiments. ivermectin exhibited solubility values ranging from 0.130 × 10−5 to 1.760 × 10−4 in the binary system and 1.880 × 10−5 to 1.237 × 10−4 in the ternary system. The outcome indicated that the addition of ethanol significantly increased the mole fraction of ivermectin in CO2. The highest impact of ivermectin solubility was found in the ivermectin -Ethanol- CO2 system at 338 K and 12 MPa, which was about 14.35 times greater than binary system under the same conditions. The models proposed by Jouyban et al. and Sodeifian-Sajadian models showed the best correlation with average absolute relative deviation (AARD%) and Akaike information criterion (AICc) for binary and ternary approaches, respectively. The results confirmed that rapid expansion supercritical solution (RESS) and gas antisolvent (GAS) method can be apply for producing nanoparticles at the suitable ranges of the drug solubility.