x is the simulated time series.
n is the number of the simulated points.
level is the noise standard deviation divided by the standard deviation of the
noise-free time series. We assume Gaussian noise with zero mean.
a, b, c, and s are the parameter
x0 is the initial values vector for x.