In this paper, the operational matrix of Euler functions for fractional derivative of order β in the Caputo sense is derived. Via this matrix, we develop an efficient collocation method for solving nonlinear fractional Volterra integro-differential equations. Illustrative examples are given to demonstrate the validity and applicability of the proposed method, and the comparisons are made with the existing results.