In this letter, complex operational matrices of Euler functions and their interesting properties are obtained to provide a novel method for solving linear complex differential equations under mixed initial conditions. Convergence conditions of this method are studied in depth and numerical experiments show the efficiency of this method. In addition, reasonable numerical results are obtained by selecting small number of basis functions.