The Mathematica notebooks calculate mean values of the energy density of the quantum fields under curved background. The background of the uniform isotropic expanding universe is implied. The notebooks densitySec2.nb and densitySec3.nb are related to the methods described in the second and third sections of the paper respectively. These notebooks give the same results.

Operator of the energy density acts to the state of a system taken in the form of the wave packet. Then mean value is culculated by the integration. Only one k-mode of the quantum field is considered. 

These codes work on Mathematica 9,10,11,12 versions.