#include <iostream>

using namespace std;

#include "util/util.h"

#define SAMPLE 50000
#define MU 2.9E-6
#define MU2 1E-3

int main(){
 
  double x_n(0), x_n1(0), x_n2(0);
  double a_n(rand_regularized(0, MU)), a_n1(rand_regularized(0, MU));
  double x(0);
  
  for(int i = 0; i < SAMPLE; i++){
    x_n = (x_n1 * 2 - x_n2) * (1.0 - 8E-6) + a_n + 0.286 * a_n1;
    x = x_n + rand_regularized(0, MU2);
    cout << x << endl;
    x_n2 = x_n1;
    x_n1 = x_n;
    a_n1 = a_n;
    a_n = rand_regularized(0, MU);
  }
  
  return 0;
}