In recent years, progress toward the classification of superintegrable systems with higher order integrals of motion has been made. In particular, a complete classification of all exotic potentials with a third or a fourth order integrals, and allowing separation of variables in Cartesian coordinates. All doubly exotic potentials with a fifth order integral have also been completely classified. It has been demonstrated how the Chazy class of third order differential equations plays an important role in solving determining equations. Moreover, taking advantage of various operator algebras defined as Abelian, Heisenberg, Conformal and Ladder case of operator algebras, we re-derived these models. These new techniques also provided further examples of superintegrable Hamiltonian with integrals of arbitrary order. It has been conjectured that all quantum superintegrable potentials that do not satisfy any linear equation satisfy nonlinear equations having the Painlevé property. In addition, it has been discovered that their integrals naturally generate finitely generated polynomial algebras and the representations can be exploited to calculate the energy spectrum. For certain very interesting cases associated with exceptional orthogonal polynomials, these algebraic structures do not allow to calculate the full spectrum and degeneracies. It has been demonstrated that alternative sets of integrals which can be build and used to provide a complete solution. This this allow to make another conjecture i.e. that higher order superintegrable systems can be solved algebraically, they require alternative set of integrals than the one provided by a direct approach.consisted in obtaining the determining equations (N=3) from a third order integral [48] and classifying all Hamiltonians allowing separation of variable in Cartesian coordinates with a third order integrals [49]. A connection with Painlevé transcendents was also established for the first time and demonstrated the significance of studying higher order superintegrability.Much recently, it has been recognized how the search can be narrowed to exotic potentials i.e. the potential do not satisfy any linear differential equation. All exotic potentials for integrals of fourth order (N=4) [67] and doubly exotic systems with integral of fifth order (N=5) [2] were obtained. These potentials can be rewritten in term of the first, second, third, fourth and fifth Painlevé transcendents. Superintegrable systems allowing separation of variables in parabolic coordinates for N=3 [88], as well as systems allowing separation of variables in polar coordinates for N=3 [94] and N=4 [29,30]. Models were obtained in terms of the sixth Painlevé transcendents. Some work have been done in regard of the case of Nth order integrals related to superintegrable systems with separation of variables in Cartesian coordinates [66] and [89]. In particular various alternative form for the integrals of motion have been obtained. Similar work for Nth order integrals have also been done in regard of polar coo...