In this work, the vapor-liquid-liquid equilibrium (VLLE) of the water/n-pentane, water/n-hexane, water/n-octane, and water/n-decane binary systems is calculated by computer simulation using the NVT-Gibbs ensemble (in the version of three simulation boxes) combined with the configurational bias Monte Carlo method. The combination of both methods, the molecular potential models used, and the simulation details allowed us to calculate the triphasic equilibrium properties of the systems studied: the densities of the three phases in equilibrium, their compositions, and potential energies. In previous works, these simulations were not carried out at a temperature range nor water/n-alkanes systems simulated in this work, probably because they are highly nonideal systems; so, to the best of our knowledge, this is the first time that this phenomenon is studied in detail. The results from VLLE simulations of the water/n-pentane system for temperatures from 343.2 to 435 K, the water/n-hexane system for temperatures from 373.11 to 473.15 K, the water/n-octane system for temperatures from 310.9 to 500 K, and for the water/n-decane system for temperatures from 374.15 to 525 K are reported here. The temperature range was selected in concordance with the experimental data available for an adequate study of the VLLE simulation results. The subcritical densities (vapor and liquid rich in n-alkane phases) at various temperatures fit well with the scaling law and the law of rectilinear diameters, allowing the estimation of upper critical end point temperature and density of the VLLE. The simulation results show a good prediction with experimental data reports in the literature.