An incompressible variational ideal ballooning mode equation is discretized with the COOL finite element discretization scheme using basis functions composed of variable order Legendre polynomials. This reduces the second order ordinary differential equation to a special block pentadiagonal matrix equation that is solved using an inverse vector iteration method. A benchmark test of BECOOL (Ballooning Eigensolver using COOL finite elements) with second order Legendre polynomials recovers precisely the eigenvalues computed by the VVBAL shooting code. Timing runs reveal the need to determine an optimal lower order case. Eigenvalue convergence runs show that cubic Legendre polynomials construct the optimal ballooning mode equation for intensive computations.