EN
The bound state solution of Coulomb Potential in the Dirac equation is calculated for a position dependent mass function M(r) within the framework of the asymptotic iteration method (AIM). The eigenfunctions are derived in terms of hypergeometric function and function generator equations of AIM.