The electron density of monoclinic paracetamol was derived from high-resolution X-ray diffraction at 100 K. The Hansen-Coppens multipole model was used to refine the experimental electron density. The topologies of the electron density and the electrostatic potential were carefully analyzed. Numerical and analytical procedures were used to derive the charges integrated over the atomic basins. The highest charge magnitude (-1.2 e) was found for the N atom of the paracetamol molecule, which is in agreement with the observed nucleophilic attack occurring in the biological media. The electric field generated by the paracetamol molecule was used to calculate the atomic charges using the divergence theorem. This was simultaneously applied to estimate the total electrostatic force exerted on each atom of the molecule by using the Maxwell stress tensor. The interaction electrostatic energy of dimers of paracetamol in the crystal lattice was also estimated.