This paper proposes a new recursive based polynomial approach for modeling a hysteresis of a piezostack actuator. The formulation is based on two main curves of hysteresis that are identified experimentally. In the proposed model, an adjusting function that relates the outer curve and its next minor is given. Its coefficients are obtained from the boundary conditions and experiment. From this adjusting function, a recursive-based polynomial formula for the cases of monotonic increasing and decreasing of input excitation are developed to predict the hysteresis. To evaluate the accuracy of the proposed new model, experiments with two different waveforms are carried out.