In the current study, an effective numerical technique is presented for solving a class of three-dimensional fractional differential-integral equations via the three-dimensional block-pulse functions. The Caputo fractional derivative is applied to deal with the fractional derivative terms, and the differential operational matrices of 3D-BPFs are utilized to transform the original system. In the end of the paper, two numerical examples are presented to verify the proposed method.