#!/usr/bin/env python
from __future__ import print_function
import os, sys, math

path = os.path.dirname(__file__)
path = os.path.dirname(path)
sys.path.insert(0, path)
if len(sys.argv) > 1:
    sys.stdout = open(sys.argv[1], "w")

import liba  # type: ignore

Ts = 1.0
f_0 = 1.0
f_e = 1000.0
c = (f_e - f_0) / Ts
T = 1.0 / 10000
t = 0
x = []
while t <= Ts:
    x.append(t)
    t += T
hpf = liba.hpf(10, 0.01)
for i in range(len(x)):
    x1 = math.cos(2 * math.pi * (f_0 * x[i] + 0.5 * c * x[i] * x[i]))
    print("%g,%g,%g" % (x[i], x1, hpf(x1)))

hpf.zero()
hpf.alpha
hpf.output
hpf.input