“…The application of fractional optimal control problems can be seen in engineering and physics and the aim of solving an optimal control problem is extremizing a cost function over an admissible set of control and state functions. Several numerical methods are applied to find an approximate solution to one-dimensional fractional optimal control problems, such as eigen functions method (Agrawal, 2008), rational approximation method (Tricaud and Chen, 2010), Legendre orthonormal basis method (Lotfi et al, 2013), Legendre operational technique (Bhrawy and Ezz-Eldien, 2016), Bernoulli polynomials method (Rabiei et al, 2018b), Hybrid of block-pulse functions and Bernoulli polynomials (Mashayekhi and Razzaghi, 2018), hybrid Chelyshkov functions (Mohammadi et al, 2018), fractional order Lagrange polynomials (Sabermahani et al, 2019), Adomian decomposition method (Alizadeh and Effati, 2018), Chebysheve collocation method (Rabiei and Parand, 2019), Grunwald-Letnikov, trapezoidal and Simpson fractional integral formulas (Salati et al, 2019), low dimensional approximations (Peng et al, 2019) and Spectral Galerkin approximation (Zhang and Zhou, 2019). But there are few researches devoted to two-dimensional problem especially in fractional area; for example, the authors in Nemati and Yousefi (2017) used the Ritz method to solve a class of these problems.…”