The proposed paper is devoted to presenting and analyzing a kinetic model describing the development of autoimmune disorders. The proposed model is a nonlinear system of differential equations that considers the biological activity of the interacting populations. The main characteristics of autoimmune diseases are taken into account. Preliminaries to the research area are provided. The modeling problem is discretized and solved approximately. The numerical results illustrate typical outcomes of autoimmune diseases.
