Code for getting stable/equilibrium configuration of interacting Wigner Molecules starting from random initial configuration. This code will calculate the positional and orientational correlation between the interacting particles, along with Lindemann ratio and nearest neighbors to estimate the different topological defects like dislocation and disclination.