Two algorithms of computing stress increment by using the elasto-plasticity constitutive model are firstly formulated, which are the Euler integration method and the radial ret method.Hill’48 anisotropic yield criterion is used. The Euler integration method can not obtain more accurate computation of the stress increment as the radial ret method unless enough subintervals are taken,by which the Euler integration method will take excessive computing time. Without decreasing any accuracy, the radial ret method can save much time. Finally, a square cup deep drawing from NUMISHEET’93 benchmarks is simulated with a self-developed code SheetForm in order to investigate the accuracy and efficiency of the radial ret method.