BackgroundCancer, a complex and deadly health concern today, is characterized by forming potentially malignant tumors or cancer cells. The dynamic interaction between these cells and their environment is crucial to the disease. Mathematical models can enhance our understanding of these interactions, helping us predict disease progression and treatment strategies.MethodsIn this study, we develop a fractional tumor-immune interaction model specifically for lung cancer (FTIIM-LC). We present some definitions and significant results related to the Caputo operator. We employ the generalized Laguerre polynomials (GLPs) method to find the optimal solution for the FTIIM-LC model. We then conduct a numerical simulation and compare the results of our method with other techniques and real-world data.ResultsWe propose a FTIIM-LC model in this paper. The approximate solution for the proposed model is derived using a series of expansions in a new set of polynomials, the GLPs. To streamline the process, we integrate Lagrange multipliers, GLPs, and operational matrices of fractional and ordinary derivatives. We conduct a numerical simulation to study the effects of varying fractional orders and achieve the expected theoretical results.ConclusionThe findings of this study demonstrate that the optimization methods used can effectively predict and analyze complex phenomena. This innovative approach can also be applied to other nonlinear differential equations, such as the fractional Klein–Gordon equation, fractional diffusion-wave equation, breast cancer model, and fractional optimal control problems.