  We construct a program to calculate Feynman amplitudes at finite temperature
in the real time Keldysh formalism using the symbolic manipulation program {\it
Mathematica}. As an example, the usefulness of this program is demonstrated by
proving the finite temperature Ward identity for QED in a second order
effective theory.
