“…If non-linearities are taken into account, Itô calculus has been used to develop time-domain solution [8][9][10], while the concepts of Volterra and Wiener series expansion have been widely applied in the frequency domain [11]. Monte Carlo simulation procedures can be formulated in the frequency domain, applying the usual superposition harmonic technique [12] or its more sophisticated derivations [13,14], as well as in the time domain by building suitable recursive models such as AR, MA and ARMA [15,16].…”