Abstract The dependence of eigenvalues of Dirac system with general boundary conditions is studied. It is shown that the eigenvalues of Dirac operators depend not only continuously but also smoothly on the coefficients, the boundary conditions, and the endpoints of the problem. Furthermore, the differential expressions of the eigenvalues as regards these parameters are given. The results obtained in this paper would provide theoretical support for the numerical calculations of eigenvalues of the corresponding problems.