Eigenvalue problems with eigenparameter appearing in the boundary conditions usually have complicated characteristic determinant where zeros cannot be explicitly computed. In this paper we use the derivative sampling theorem ‘Hermite interpolations’ to compute approximate values of the eigenvalues of Dirac systems with eigenvalue parameter in one or two boundary conditions. We use recently deri...