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