A semi-inverse variational method for generating the bound state energy eigenvalues in a quantum system: the Dirac Coulomb type-equation