Seismic anisotropy of polycrystalline materials depends on the characteristics of microcracks as well as the crystallographic orientations of minerals. Here, we present a MATLAB-based software, AnisEulerSC (Anisotropy from Euler angles using Self-Consistent approximation), for modeling the anisotropic seismic properties of a polycrystalline aggregate with microcracks using the self-consistent approximation. In this program, several commands of MTEX (a MATLAB toolbox for texture analysis) are utilized to analyze constituent components of a rock sample based on experimental measurement data (e.g., electron backscatter diffraction) and to calculate and visualize seismic wave velocities and polarizations. We provide several examples of SC modeling to explore the effects of grain shapes on a polycrystalline aggregate and the effects of crack shapes and orientations on the anisotropic seismic properties of a polycrystalline aggregate.