#include using namespace std; #include "util/util.h" #define SAMPLE 50000 #define MU 6E-1 #define MU2 5E-1 int main(){ double x_n(0), x_n1(0); double a_n(rand_regularized(0, MU)); double x(0); for(int i = 0; i < SAMPLE; i++){ x_n = x_n1 * (1. - 3E-3) + a_n; x = x_n + rand_regularized(0, MU2); cout << x << endl; x_n1 = x_n; a_n = rand_regularized(0, MU); } return 0; }