We study a simple first-order nonnegative bilinear time-series model and give conditions under which the model is stationary. The probability density function of the stationary distribution (when it exists) is found. We also discuss the tail behaviour of the stationary distribution and calculate the probability density function by a numerical method. Simulation is used to check the calculation.