Zeros and Their Distribution for Leonardo Polynomials

Abstract

We investigate the zeros of Leonardo polynomials and their distribution in the complex plane. We show that the structure of the zeros differs between even and odd degrees. For even degrees, all zeros are simple and purely imaginary, with imaginary parts lying in \((-2,2)\). For odd degrees, there is exactly one positive real zero, whereas no negative real or purely imaginary zeros occur. We obtain two-sided bounds for the unique positive real zero, prove that it decreases strictly to zero as the degree increases through odd values, and determine its asymptotic behavior. Using the relation between Leonardo and Fibonacci polynomials together with a complex trigonometric parametrization, we transform the zero equation into a family of fixed-point problems. Suitable invariant disks are constructed, and the corresponding mappings are shown to be contractions. This yields a constructive description of the complete zero set, proves the simplicity of all zeros, and provides geometrically convergent iterations for their numerical approximation. The method also provides explicit control over the approximation error. Finally, we establish a degree-dependent vertical strip containing all zeros and show that its width tends to zero as the degree increases. Consequently, the zeros become uniformly concentrated near the imaginary axis. Numerical examples illustrate the convergence of the proposed method, the associated error estimates, and the resulting zero distributions.