A finite element program for elastic-plastic stress analysis based on a new approximation of the constitutive equations is presented. A brief introduction to the theory of this new procedure and a detailed explanation of the subroutines that are essential for this method are also described.