“…Now, we implement a new fractional pseudospectral method based on the FL and SL polynomials. In this method, the state and control variables of FHOC problem () are approximated as follows where and for j = 1, 2,…, N are unknown variables, and for k = 1, 2,…, N are arbitrary interpolating polynomials (they can be considered as classical Lagrange polynomials or other types of fractional Lagrange polynomials 12,13 ). By interpolation property, we get the following relations It is easy to show that …”