function [b,a] = calcCoeffs(Fs, F0)
    alpha = 2*pi*F0/Fs;
    a = [1, -2*cos(alpha), 1];
    b = [1,-cos(alpha), 0];
end