Tis research endeavors to introduce novel fractional pseudospectral methodologies tailored for addressing fractional optimal control problems encompassing inequality constraints and boundary conditions. Leveraging fractional Lagrange interpolation functions, we formulate diferential and integral pseudospectral matrices pivotal in discretizing fractional optimal control problems. Te stability of these matrices is ensured through the utilization of M¨untz–Legendre polynomials. Discretization of fractional optimal control problems entails employing integral and diferential fractional matrices, with collocation strategically positioned at both Gauss and fipped Radau-type points. Our results demonstrate that the proposed method is well suited for practical problems with larger domains. Furthermore, it proves efective for Bang-Bang type problems and ofers substantial benefts for problems with nonsmooth solutions. Comprehensive numerical evaluations on benchmark fractional optimal control problems substantiate the efectiveness of the devised pseudospectral methodologies, showcasing their commendable performance and potential for practical applications.