The problem of finding analytical solutions for time-dependent or time-independent heat equations, especially for solid-state laser media, has required a lot of work in the field of thermal effects. However, to calculate the temperature distributions analytically, researchers often have to make some approximations or employ complex methods. In this work, we present full analytical solutions for anisotropic time-dependent heat equations in the Cartesian coordinates with various source terms corresponding to various pumping schemes. Moreover, the most general boundary condition of Robin (or impedance boundary condition), corresponding to the convection cooling mechanism, was applied. This general condition can be easily switched to constant temperature and thermal insulation as special cases. To this end, we first proposed a general approach to solving time-dependent heat equations with an arbitrary heat source. We then applied our approach to explore the temperature distribution for three cases: steady-state pumping or long transient, single-shot pumping or short transient, and repetitively pulsed pumping. Our results show the possibility of an easier and more accurate approach to analytical calculations of the thermal dispersion, thermal stresses (strains), thermal bending, thermal phase shift, and other thermal effects.