A novel study for solving systems of nonlinear fractional integral equations

In this study, we explore the solution of a nonlinear system of fractional integro-differential equations based on the operational matrix method. We have modified the operational matrix method to accommodate such systems and have streamlined the resulting algebraic system in a practical manner. Instead of dealing with a system of Formula: see text equations, we have simplified the problem to finding the solution for a set of Formula: see text nonlinear equations. Here, m represents the number of block pulse functions, which is typically a large positive integer. This approach significantly reduces computational costs, especially for large values of m. We provide proofs for both the existence and uniqueness of the solution for this system. Additionally, we investigate two types of stability: Ulam–Hyers stability and generalized Ulam–Hyers stability. We have evaluated the effectiveness of the proposed method by applying it to three different problems and comparing our results with existing literature. The outcomes of these tests highlight the efficiency and efficacy of our proposed method.