“…Due to the nonlinearity, this differential equation cannot be solved in an elementary way, unless it is approximated that the column deflection is exceedingly small. Therefore, various methods have been proposed for solving it to realize high-accuracy analysis of columnar buckling behaviors under different mechanical conditions, including the geometry, loading conditions, and boundary conditions [2,. Most methods have been based on either the Bessel function, elliptic integral, or other nonelementary functions/integrals, ensuring high accuracy in the analysis of column buckling behavior even with large deflections; those based on generalized hypergeometric functions [29] and the singular perturbation method [31] are only a few to mention.…”