Calculations of multiquark functions in effective models of strong interaction

In this paper we present our results of the investigation of multiquark equations in the Nambu-Jona-Lasinio model with chiral symmetry of SU(2) group in the mean-field expansion. To formulate the mean-field expansion we have used an iteration scheme of solution of the Schwinger-Dyson equations with the fermion bilocal source. We have considered the equations for Green functions of the Nambu-Jona-Lasinio model up to third step for this iteration scheme. To calculate the high-order corrections to the mean-field approximation, we propose the method of the Legendre transformation with respect to the bilocal source, which allows effectively to take into account the symmetry constraints related with the chiral Ward identity. We discuss also the problem of calculating the multiquark functions in the mean-field expansion for Nambu-Jona-Lasinio-type models with other types of the multifermion sources.