We present a Legendre pseudospectral method for directly estimating the costates of the Bolza problem encountered in optimal control theory. The method is based on calculating the state and control variables at the Legendre‐Gauss‐Lobatto (LGL) points. An Nth degree Lagrange polynomial approximation of these variables allows a conversion of the optimal control problem into a standard nonlinear programming (NLP) problem with the state and control values at the LGL points as optimization parameters. By applying theKarush ‐Kuhn‐Tucker (KKT) theorem to the NLP problem, we show that the KKT multipliers satisfy a discrete analog of the costate dynamics including the transversality conditions. Indeed, we prove that the costates at the LGL points are equal to the KKT multipliers divided by the LGL weights. Hence, the direct solution by this method also automatically yields the costates by way of the Lagrange multipliers that can be extracted from an NLP solver. One important advantage of this technique is that it allows a very simple way to check the optimality of the direct solution. Numerical examples are included to demonstrate the method.