“…Based on the analysis, we developed an ordinary differential equations (ODEs) model to investigate the dynamic processes of the immune response during IL-2 treatment in our model ( Khailaie et al., 2020 ; Sontag, 2017 ). , , , , and denote the density of APC, TN, Treg d , Treg f , and NK cell, respectively, they all are in the unit mm −3 .…”