A variational pseudospectral method based on dual variational principles is proposed.
The proposed method is highly efficient due to the benefits the sparse and symmetric properties.
Solutions of high precisions can be obtained with few time nodes.
Extra costate estimations are not required.