A computational model is presented to calculate the ground state energy of neutral and charged excitons confined in semiconductor quantum dots. The model is based on the variational Quantum Monte Carlo method and effective mass Hamiltonians. Through an iterative Newton–Raphson process, minimizing the local energy, and (optional) parallelization of random walkers, fast and accurate estimates of both confinement and Coulomb binding energies can be obtained in standard desktop computers. To illustrate the reach of the model, we provide Fortran programs and illustrative calculations for colloidal CdSe nanoplatelets with large lateral dimensions and dielectric confinement, where electronic correlations are strong. The results compare well with exact variational calculations and largely outperform configuration interaction calculations in computational efficiency. Program summaryProgram title: vqmc-emaCPC Library link to program files:https://doi.org/10.17632/dtpyc2pffg.1Licensing provisions: GPLv3Programming language: Fortran 90Nature of problem: Calculation of either exciton or trion ground state energy and wave function in a cuboidal semiconductor nanoplatelet. Hard wall quantum confinement coexists with dielectric confinement (self-energy and Coulomb polarization terms).Solution method: Variational quantum Monte Carlo with effective mass Hamiltonians, integrated into a Newton–Raphson solver. OMP parallelization library can be (optionally) linked.