Tricritical points of mixtures are characterized by thermodynamic coordinates where three coexisting phases become identical. Here, we formulate and solve for the first time the tricritical-point calculation as a global optimization problem. Using finite-difference formulas and the ordinary criticality criteria, we develop novel numerical approximations for tricriticality conditions associated with the derivatives of the fourth and fifth orders with respect to the compositions, obtained from the Gibbs tangent plane criterion for phase stability. Such approximations simplified the numerical implementations of these two conditions, avoiding propagation errors related to floating-point arithmetic without requiring a code with excessive numerical precision, which usually leads to high computational costs. This paper deals only with ternary mixtures. Thus, the present formulation is described as a minimization problem with strict inequality constraints, having four variables (two independent mole fractions, temperature, and pressure). Using the Peng−Robinson equation of state together with the classical one-fluid van der Waals mixing rule, we observe that the novel procedure is able to determine the tricritical coordinates in a set of ternary mixtures formed by alkanes or alkanes and light gases, which were previously studied by other authors and generally have complex multiphase behavior. The proposed method is compared to previous ones, and the results presented here compare favorably with the experimental values, which are available in the literature. In addition, we present an unprecedented prediction using the Peng−Robinson equation. We also show an example which leads to a no prediction of the tricritical point of a polar mixture, whose prediction (with this equation of state) was considered possible in a previous work.