In this paper, the results of the investigation of multi-quark equations in the Nambu–Jona-Lasinio (NJL) model in the mean-field expansion are presented. The multi-quark functions have been considered up to the third order of expansion. One of the purposes of our computations is the study of corrections of higher orders to parameters of the model. The important problem of the application of the NJL model is regularization. We compare the NJL model with 4-dimensional cutoff regularization and the dimensionally analytical regularization. We also discuss so-called “predictive regularization” in the NJL model, and a modification of this regularization, which is free of the Landau pole, is proposed. To calculate the high-order corrections, we use the Legendre transform method in the framework of bilocal-source formalism, which allows one effectively to take into consideration the symmetry constraints. A generalization of the mean-field expansion for other types of multi-quark sources is also discussed.