В работе рассматривается математическая модель иммунотерапии рака. Модель является модификации трехмерной модели иммунного ответа введением уравнения, описывающем динамику иммунотерапии. Итоговая модель представляет собой систему из четырех обыкновенных дифференциальных уравнений, описывающих динамику изменения раковых, здоровых, иммунных клеток и взаимодействие иммунных клеток с интерлейкином. В работе решается задача стабилизации количества раковых клеток с применением блочного подхода с использованием линеаризации обратной связью.