We study the energy density of shaped waves inside a quasi-1D disordered waveguide. We find that the spatial energy density of optimally shaped waves, when expanded in the complete set of eigenfunctions of the diffusion equation, is well described by considering only a few of the lowest eigenfunctions. Taking into account only the fundamental eigenfunction, the total internal energy inside the sample is underestimated by only 2%. The spatial distribution of the shaped energy density is very similar to the fundamental eigenfunction, up to a cosine distance of about 0.01. We obtain the energy density of transmission eigenchannels inside the sample by numerical simulation of the scattering matrix. Computing the transmission-averaged energy density over all transmission channels yields the ensemble averaged energy density of shaped waves. From the averaged energy density, we reconstruct its spatial distribution using the eigenfunctions of the diffusion equation. The results of our study have exciting applications in controlled biomedical imaging, efficient light harvesting in solar cells, enhanced energy conversion in solid-state lighting, and low threshold random lasers.