/******************************************************************************
 
                              Online C++ Compiler.
               Code, Compile, Run and Debug C++ program online.
Write your code in this editor and press "Run" button to compile and execute it.
 
*******************************************************************************/
 
#include <stdint.h>
#include <stdbool.h>
#include <stdio.h>
#include <math.h>
 
using namespace std;
 
#define PI				3.14159265358979323846
 
static double Sinc(double x)
{
	if (x == 0.0) return 1.0;
	return sin(PI * x) / (PI * x);
}
 
static double SymmetricBlackmanHarris(double i, int n)
{
	i += (double)n / 2.0;
	double factor = i * 2.0 * PI / (double)n;
	const double a0 = 0.35875;
	const double a1 = 0.48829;
	const double a2 = 0.14128;
	const double a3 = 0.01168;
	return a0 - a1 * cos(1.0 * factor) + a2 * cos(2.0 * factor) - a3 * cos(3.0 * factor);
}
 
static double SymmetricBlackman(double i, int n)
{
	i += (double)n / 2.0;
	double factor = i * 2.0 * PI / (double)n;
	// Constants from https://e...content-available-to-author-only...a.org/wiki/Window_function#Blackman_window
	const double a0 = 0.42;
	const double a1 = 0.5;
	const double a2 = 0.08;
	return a0 - a1 * cos(1.0 * factor) + a2 * cos(2.0 * factor);
}
 
static double SymmetricHamming(double i, int n)
{
	i += (double)n / 2.0;
	double factor = i * 2.0 * PI / (double)n;
	const double a0 = 25.0 / 46.0;
	const double a1 = 1.0 - a0;
	return a0 - a1 * cos(1.0 * factor);
}
 
static double SymmetricHann(double i, int n)
{
	i += (double)n / 2.0;
	double factor = i * 2.0 * PI / (double)n;
	const double a0 = 0.5;
	const double a1 = 1.0 - a0;
	return a0 - a1 * cos(1.0 * factor);
}
 
int main()
{
	static constexpr int FIRN = 52;
	float tblF[FIRN];
    const double cutoff = 1.0;
 
    printf("Hello World\n");
 
	double sum = 0;
	for (size_t i = 0; i < FIRN; i++)
	{
		double t = double(i) - double(FIRN / 2.0) + 1.0 - 0.5;
		sum += tblF[i] = SymmetricBlackmanHarris(t, FIRN) * cutoff * Sinc(cutoff * t);
	}
	// Normalize
	for (size_t i = 0; i < FIRN; i++)
		tblF[i] *= 1.0 / sum;
 
	for (size_t i = 0; i < FIRN; i++)
	    printf("%.20g\n", tblF[i]);
 
    return 0;
}
 

LyoqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKgoKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgT25saW5lIEMrKyBDb21waWxlci4KICAgICAgICAgICAgICAgQ29kZSwgQ29tcGlsZSwgUnVuIGFuZCBEZWJ1ZyBDKysgcHJvZ3JhbSBvbmxpbmUuCldyaXRlIHlvdXIgY29kZSBpbiB0aGlzIGVkaXRvciBhbmQgcHJlc3MgIlJ1biIgYnV0dG9uIHRvIGNvbXBpbGUgYW5kIGV4ZWN1dGUgaXQuCgoqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqLwoKI2luY2x1ZGUgPHN0ZGludC5oPgojaW5jbHVkZSA8c3RkYm9vbC5oPgojaW5jbHVkZSA8c3RkaW8uaD4KI2luY2x1ZGUgPG1hdGguaD4KCnVzaW5nIG5hbWVzcGFjZSBzdGQ7CgojZGVmaW5lIFBJCQkJCTMuMTQxNTkyNjUzNTg5NzkzMjM4NDYKCnN0YXRpYyBkb3VibGUgU2luYyhkb3VibGUgeCkKewoJaWYgKHggPT0gMC4wKSByZXR1cm4gMS4wOwoJcmV0dXJuIHNpbihQSSAqIHgpIC8gKFBJICogeCk7Cn0KCnN0YXRpYyBkb3VibGUgU3ltbWV0cmljQmxhY2ttYW5IYXJyaXMoZG91YmxlIGksIGludCBuKQp7CglpICs9IChkb3VibGUpbiAvIDIuMDsKCWRvdWJsZSBmYWN0b3IgPSBpICogMi4wICogUEkgLyAoZG91YmxlKW47Cgljb25zdCBkb3VibGUgYTAgPSAwLjM1ODc1OwoJY29uc3QgZG91YmxlIGExID0gMC40ODgyOTsKCWNvbnN0IGRvdWJsZSBhMiA9IDAuMTQxMjg7Cgljb25zdCBkb3VibGUgYTMgPSAwLjAxMTY4OwoJcmV0dXJuIGEwIC0gYTEgKiBjb3MoMS4wICogZmFjdG9yKSArIGEyICogY29zKDIuMCAqIGZhY3RvcikgLSBhMyAqIGNvcygzLjAgKiBmYWN0b3IpOwp9CgpzdGF0aWMgZG91YmxlIFN5bW1ldHJpY0JsYWNrbWFuKGRvdWJsZSBpLCBpbnQgbikKewoJaSArPSAoZG91YmxlKW4gLyAyLjA7Cglkb3VibGUgZmFjdG9yID0gaSAqIDIuMCAqIFBJIC8gKGRvdWJsZSluOwoJLy8gQ29uc3RhbnRzIGZyb20gaHR0cHM6Ly9lLi4uY29udGVudC1hdmFpbGFibGUtdG8tYXV0aG9yLW9ubHkuLi5hLm9yZy93aWtpL1dpbmRvd19mdW5jdGlvbiNCbGFja21hbl93aW5kb3cKCWNvbnN0IGRvdWJsZSBhMCA9IDAuNDI7Cgljb25zdCBkb3VibGUgYTEgPSAwLjU7Cgljb25zdCBkb3VibGUgYTIgPSAwLjA4OwoJcmV0dXJuIGEwIC0gYTEgKiBjb3MoMS4wICogZmFjdG9yKSArIGEyICogY29zKDIuMCAqIGZhY3Rvcik7Cn0KCnN0YXRpYyBkb3VibGUgU3ltbWV0cmljSGFtbWluZyhkb3VibGUgaSwgaW50IG4pCnsKCWkgKz0gKGRvdWJsZSluIC8gMi4wOwoJZG91YmxlIGZhY3RvciA9IGkgKiAyLjAgKiBQSSAvIChkb3VibGUpbjsKCWNvbnN0IGRvdWJsZSBhMCA9IDI1LjAgLyA0Ni4wOwoJY29uc3QgZG91YmxlIGExID0gMS4wIC0gYTA7CglyZXR1cm4gYTAgLSBhMSAqIGNvcygxLjAgKiBmYWN0b3IpOwp9CgpzdGF0aWMgZG91YmxlIFN5bW1ldHJpY0hhbm4oZG91YmxlIGksIGludCBuKQp7CglpICs9IChkb3VibGUpbiAvIDIuMDsKCWRvdWJsZSBmYWN0b3IgPSBpICogMi4wICogUEkgLyAoZG91YmxlKW47Cgljb25zdCBkb3VibGUgYTAgPSAwLjU7Cgljb25zdCBkb3VibGUgYTEgPSAxLjAgLSBhMDsKCXJldHVybiBhMCAtIGExICogY29zKDEuMCAqIGZhY3Rvcik7Cn0KCmludCBtYWluKCkKewoJc3RhdGljIGNvbnN0ZXhwciBpbnQgRklSTiA9IDUyOwoJZmxvYXQgdGJsRltGSVJOXTsKICAgIGNvbnN0IGRvdWJsZSBjdXRvZmYgPSAxLjA7CgogICAgcHJpbnRmKCJIZWxsbyBXb3JsZFxuIik7CgoJZG91YmxlIHN1bSA9IDA7Cglmb3IgKHNpemVfdCBpID0gMDsgaSA8IEZJUk47IGkrKykKCXsKCQlkb3VibGUgdCA9IGRvdWJsZShpKSAtIGRvdWJsZShGSVJOIC8gMi4wKSArIDEuMCAtIDAuNTsKCQlzdW0gKz0gdGJsRltpXSA9IFN5bW1ldHJpY0JsYWNrbWFuSGFycmlzKHQsIEZJUk4pICogY3V0b2ZmICogU2luYyhjdXRvZmYgKiB0KTsKCX0KCS8vIE5vcm1hbGl6ZQoJZm9yIChzaXplX3QgaSA9IDA7IGkgPCBGSVJOOyBpKyspCgkJdGJsRltpXSAqPSAxLjAgLyBzdW07CgoJZm9yIChzaXplX3QgaSA9IDA7IGkgPCBGSVJOOyBpKyspCgkgICAgcHJpbnRmKCIlLjIwZ1xuIiwgdGJsRltpXSk7CgogICAgcmV0dXJuIDA7Cn0K