We construct a program to calculate Feynman amplitudes at finite temperature in the real time Keldysh formalism using the symbolic manipulation program 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.