Abstract: We consider two-dimensional "artificial atoms" confined by an axially symmetric potential $V(\rho)$. Such configurations arise for systems effectively restricted to a 2D layer, such as electrons in quantum dots, in islands on liquid helium and graphene, and in planar molecules, as well as for cold ions in "pancake" traps. Using the Thomas-Fermi method, we present the first full analytic solution describing the density distribution, energy, and other parameters for any form of $V(\rho)$. An essential and nontrivial aspect of the problem is that the 2D density of states must be properly combined with 3D electrostatics. The solution turns out to have a universal form, with scaling parameters $\rho^2/R^2$ and $R/a_B^*$ ($R$ is the dot radius and $a_B^*$ is the effective Bohr radius).