Numerical analysis of 2D MHD equilibrium with non-inductive plasma current in tokamaks