A numerical method based on the Jouyban-Acree model was presented for prediction of drug solubility in water-dioxane mixtures at various temperatures. The method requires drug solubility in mono-solvent systems, i.e. two data points for each temperature of interest. The mean percentage deviation (MPD) of predicted solubilities was calculated to show the accuracy of the predicted data and 27% was found as the average MPD for 36 data sets studied. The proposed numerical method reduced the number of required experimental data from five to two points and could also be extended to predict solubility at various temperatures.