... An operational matrix to find numerical approximation of ψ-fractional differential equations (FDEs) is derived. This study extends the method to nonlinear FDEs by using quasi linearization technique to linearize the nonlinear problems. Findings The error analysis of the proposed method is discussed...