Quasilinearization method for finite systems of nonlinear RL fractional differential equations

In this paper the quasilinearization method is extended to finite systems of Riemann-Liouville fractional differential equations of order 0 < q < 1. Existence and comparison results of the linear Riemann-Liouville fractional differential systems are recalled and modified where necessary. Using upper and lower solutions, sequences are constructed that are monotonic such that the weighted sequences converge uniformly and quadratically to the unique solution of the system. A numerical example illustrating the main result is given.