In this paper, we discuss a Mathematica program that we have written that calculates the integrand for amplitudes in the closed-time-path formulation of real-time finite-temperature field theory. The program is designed to be used by someone with no previous experience with Mathematica. It performs contractions over tensor indices that appear in real-time finite-temperature field theory and gives the result in the 1-2, Keldysh or R/A basis. As an illustration of this program, we discuss the calculation of all 3-point ward identities in finite-temperature quantum electrodynamics with full vertices. PACS Nos.: 11.10.Wx,11.15.-q