#include <iostream>

using namespace std;

#include "util/util.h"

#define SAMPLE 50000
#define MU 2.5E-6
#define MU2 1E-8

int main(){
 
  double x_n(0), x_n1(0), x_n2(0), x_n3(0), x_n4(0);
  double a_n(rand_regularized(0, MU)), a_n1(rand_regularized(0, MU)), a_n2(rand_regularized(0, MU));
  double x(0);
  
  for(int i = 0; i < SAMPLE; i++){
    x_n = x_n1 * 3.462618155 + x_n2 * -4.410112265
        + x_n3 * 2.432313950 + x_n4 * -0.48481984
        + a_n - a_n1 * 1.83075548 + a_n2 * 0.83245797;
    x = x_n + rand_regularized(0, MU2);
    cout << x << endl;
    x_n4 = x_n3;
    x_n3 = x_n2;
    x_n2 = x_n1;
    x_n1 = x_n;
    a_n2 = a_n1;
    a_n1 = a_n;
    a_n = rand_regularized(0, MU);
  }
  
  return 0;
}