#include <iostream>
#include <cmath>
#include <vector>
#include <fstream>
#include <string>
#include <algorithm>

using namespace std;

int main(int argc, char* argv[]) {
    int ns = 300;
    int nh = 10;
    double AngularFrequency = 5.0;

    cout << "Samples: " << ns;
//    vector<double> samples;
//    const double pi = std::atan(1) * 4;
    const double pi = 3.141519;

    ofstream outFile;
    try {
        outFile.open("./triangle.txt");
    }
    catch(std::exception e) {
        cout << "Coult not open file" << endl;
    }

    for(int n = 0; n != ns; ++n) {
        double sum = 0.0;
        for(int k = 0; k != nh; ++k) {
            sum += pow(-1,k)*sin((2*k+1)*AngularFrequency*n)/sin(pow(2*k+1,2));
        }
//        cout << sum << endl;
//        samples.push_back(sum * 8 / (pi*pi));
        outFile << n << sum << endl;
    }

    cout << "done" << endl;
    // write samples to file


//    for_each(begin(samples), end(samples), [&](double d){
//        outFile << d << endl; });

    return 0;
}