We study a model of interacting vortices in a type II superconductor. In the weak coupling limit, we constructed a mean-field theory which allows us to accurately calculate the vortex density distribution inside a confining ...