The electron system of arbitrary degeneracy can be mapped onto a classical system where electrons of the same spin are assigned an additional interaction and the effect of degeneracy is taken into account through an imaginative temperature. We apply this method to electrons in quantum dots modeled as two-dimensional electron liquid confined in a finite domain by a harmonic potential. We analyze distribution functions by molecular dynamics instead of solving integral equations which are not so useful in the case without translational invariance.