Stall induced oscillations are likely instabilities in modern wind engineering systems. For their robust design, it is crucial to take parametric and load uncertainties into consideration. Presently, study of a stall flutter system under the influence of random gust is undertaken using the probability density evolution method, and qualitative changes in the response output are presented. This probability based approach when compared with a spectral tool called polynomial chaos expansion, shows better accuracy for the same level of discretisation. © 2015 Elsevier Ltd.