  We have written a {\it Mathematica} program that calculates the integrand
corresponding to any amplitude in the closed-time-path formulation of real time
statistical field theory. The program is designed so that it can be used by
someone with no previous experience with {\it Mathematica}. It performs the
contractions over the tensor indices that appear in real time statistical field
theory and gives the result in the 1-2, Keldysh or RA basis. We have used the
program to calculate the ward identity for the QED 3-point function, the QED
4-point function for two photons and two fermions, and the QED 5-point function
for three photons and two fermions. In real time statistical field theory,
there are seven 3-point functions, 15 4-point functions and 31 5-point
functions. We produce a table that gives the results for all of these
functions. In addition, we give a simple general expression for the KMS
conditions between $n$-point green functions and vertex functions, in both the
Keldysh and RA bases
