The 1.5 D code program of numerical simulation of nonlinear nonstationary processes in the klystron-type devices based on the nonstationary L.A. Vainshtein’s theory of cavity excitation and the «particle-in-cell» method for modeling of the electron beam dynamics is developed. The results of numerical simulation of the basic oscillation modes of the fourcavity klystron oscillator with the external delayed feedback are presented.