“…For recent advances in the pseudospectral methods, see, for example, Gong, Ross, and Fahroo (2016) ;Tang, Liu, and Hu (2016). Pseudospectral methods approximate the state and control variables using interpolating polynomials with specific collocation points such as Legendre-Gauss-Lobatto(LGL), Legendre-Gauss(LG) (Mehrpouya, Shamsi, & Azhmyakov, 2014) and Legendre-Radau(LR) points (Garg et al, 2010). Then, by collocating the state equations and path constraints and using the differentiation matrix, the problem is transcribed to a non-linear programming problem(NLP), which can be solved by the well-developed parameter optimization algorithms.…”