  The two-dimensional Dirac equation has been widely used in graphene physics,
the surface of topological insulators, and especially quantum scarring.
Although a numerical approach to tackling an arbitrary confining problem was
proposed several years ago, several fundamental issues must be thoroughly
understood and solved. In this work, we conceal and address these challenges
and finally develop a complete method, validated by comparison with analytical
results.
