Potential-density basis sets in axisymmetric coordinates