The purpose of the paper is to develop a model of iron wolf howling. The recorded live sounds of wolf howling are used in order to obtain parameters of the model. Spectral density of the recorded sounds is investigated by means of MUSIC method. A sum of quasipolynomials of the second order is applied as a model of wolf howling. The parameters of the model are estimated using short intervals of records (pitches). The parameters of the model are estimated by the modified Levenberg method. The final model is presented in a state space form.