  Many stochastic continuous-state dynamical systems can be modeled as
probabilistic programs with nonlinear non-polynomial updates in non-nested
loops. We present two methods, one approximate and one exact, to automatically
compute, without sampling, moment-based invariants for such probabilistic
programs as closed-form solutions parameterized by the loop iteration. The
exact method applies to probabilistic programs with trigonometric and
exponential updates and is embedded in the Polar tool. The approximate method
for moment computation applies to any nonlinear random function as it exploits
the theory of polynomial chaos expansion to approximate non-polynomial updates
as the sum of orthogonal polynomials. This translates the dynamical system to a
non-nested loop with polynomial updates, and thus renders it conformable with
the Polar tool that computes the moments of any order of the state variables.
We evaluate our methods on an extensive number of examples ranging from
modeling monetary policy to several physical motion systems in uncertain
environments. The experimental results demonstrate the advantages of our
approach with respect to the current state-of-the-art.
