Self-consistent hybrid neoclassical-magnetohydrodynamic simulations of axisymmetric plasmas