I show the first notions of how to go about calculating thermal expectation values in quantum Monte Carlo. I work out an example for the particle number operator and show expressions for the one-body density matrix. I assume that the field configurations are given. More on how to obtain those samples next time!