Abstract: In practical applications, equation of states (EOS) consists of several piecewise smooth surfaces, which leads to discontinuity at interface. As a traditional nonlinear iterative algorithm is applied to an energy equation with discontinuous EOS, it may lead to slow convergence and unphysical solutions. To overcome the difficulties, a nonlinear problem is designed, and a nonlinear iterative algorithm is proposed to solve the problem. The algorithm is fit for energy equations with discontinuous EOS of piecewise smooth functions. A parameter of energy change is defined in the algorithm so that it is unnecessary to know discontinuity position in advance. The algorithm calculates precisely net gain or leakage of energy, which can be used to assess influence of discontinuity in EOS. Typical numerical experiments verify that the algorithm converges stably, and gives physical solutions.