In this paper, we give sufficient conditions on a measurable function \(p:(0,\infty )^{n}\rightarrow [1,\infty )\) in order that harmonic analysis operators (maximal operators, Riesz transforms, Littlewood–Paley functions, and multipliers) associated with \(\alpha \) -Laguerre polynomial expansions are bounded on the variable Lebesgue space \(L^{p(\cdot )} ((0,\infty )^{n}, \mu _\alpha )\) , where \(\mathrm{{d}}\mu _\alpha (x)=2^{n}\prod _{j=1}^{n} \frac{x_{j}^{2\alpha _{j}+1} e^{-x_{j}^{2}}}{\Gamma (\alpha _{j}+1)} \mathrm{{d}}x\) , being \(\alpha =(\alpha _{1}, \dots , \alpha _{n})\in [0,\infty )^{n}\) and \(x=(x_1,\dots ,x_n)\in (0,\infty )^{n}\) .