vector<double> a;
    a[0] = 1;
    a[1] = 3;
    a[2] = 5;
    a[3] = 10;
    cout << moment(1,a, false)  << endl;