Optimal Auxiliary Functions Method for Nonlinear Dynamical Systems