Computing Lyapunov constants for random recurrences with smooth coefficients

In recent years, there has been much interest in the growth and decay rates (Lyapunov constants) of solutions to random recurrences such as the random Fibonacci sequence x_n+1=±x_n±x_n−1. Many of these problems involve nonsmooth dynamics (nondifferentiable invariant measures), making computations hard. Here, however, we consider recurrences with smooth random coefficients and smooth invariant measures. By computing discretised invariant measures and applying Richardson extrapolation, we can compute Lyapunov constants to 10 digits of accuracy. In particular, solutions to the recurrence x_n+1=x_n+c_n+1x_n−1, where the {c_n} are independent standard normal variables, increase exponentially (almost surely) at the asymptotic rate (1.0574735537…)ⁿ. Solutions to the related recurrences x_n+1=c_n+1x_n+x_n−1 and x_n+1=c_n+1x_n+d_n+1x_n−1 (where the {d_n} are also independent standard normal variables) increase (decrease) at the rates (1.1149200917…)ⁿ and (0.9949018837…)ⁿ, respectively.