The mathematical model is represented by a system of partial differential equations, based on the interactions among cancer cells, dendritic cells, CD4+ and CD8+ T cells, cytokines IL-12 and IL-2, GM-CSF, PD-1, PD-L1, PD-1-PD-L1 complex and anti-PD-1. The gene discussed is CSF2; the disease is cancer.