This report describes a set of computer programs for NAL’s method of estimation of the vibration characteristics boundaries for flutter and divergence of an aeroelastic system subjected to stationary random noises. The programs listed here are those for the calculations of the autocoariance function of a time series and Jury’s stability parameters of the characteristic equations of fourth and eighth orders, and those for the estimation of the stability boundary, frequency and damping ratio of the system.