Experimental solubilities of diazepam in binary and ternary solvents of polyethylene glycols 200 and 400 with N-methyl pyrrolidone and water at T = 298.2 K are reported. The Jouyban–Acree model was used to fit solubility data of diazepam in the binary and ternary solvent mixtures (106 data points) in which the overall mean relative deviations (OMRD %) is 13.1 % and the prediction OMRD % is 31.7 %. The combined version of the Jouyban–Acree model with Hansen solubility parameters was used for fitting and predicting the solubility data and the OMRDs % are 10.0 and 20.8 %, respectively. Also, the previously proposed trained versions of the Jouyban–Acree model were used for predicting the reported data in this work and all results are listed in the tables. The density of the solute-free solvent mixtures were measured and employed to calculate the constants of the Jouyban–Acree model and then the densities of the saturated solutions were predicted.