The global population has encountered significant challenges throughout history due to infectious diseases. To comprehensively study these dynamics, a novel deterministic mathematical model, TCD I L 2 Z, is developed for the early detection and treatment of lung cancer. This model incorporates I L 2 cytokine and anti-PD-L1 inhibitors, enhancing the immune system’s anticancer response within five epidemiological compartments. The TCD I L 2 Z model is analyzed qualitatively and quantitatively, emphasizing local stability given the limited data-a critical component of epidemic modeling. The model is systematically validated by examining essential elements such as equilibrium points, the reproduction number ( R 0 ), stability, and sensitivity analysis. Next-generation techniques based on R 0 that track disease transmission rates across the sub-compartments are fed into the system. At the same time, sensitivity analysis helps model how a particular parameter affects the dynamics of the system. The stability on the global level of such therapy agents retrogrades individuals with immunosuppression or treated with I L 2 and anti-PD-L1 inhibitors admiring the Lyapunov functions’ applications. NSFD scheme based on the implicit method is used to find the exact value and is compared with Euler’s method and RK4, which guarantees accuracy. Thus, the simulations were conducted in the MATLAB environment. These simulations present the general symptomatic and asymptomatic consequences of lung cancer globally when detected in the middle and early stages, and measures of anticancer cells are implemented including boosting the immune system for low immune individuals. In addition, such a result provides knowledge about real-world control dynamics with I L 2 and anti-PD-L1 inhibitors. The studies will contribute to the understanding of disease spread patterns and will provide the basis for evidence-based intervention development that will be geared toward actual outcomes.