%Kernal 1
Fs = 33600
n = 1:1:33600;
a = floor(n-(Fs/2));
c = (1200*pi)/Fs; % c is greek letter chi
g = pi/Fs; % g is greek symbol gamma
% To evaluate the function, you need to consider the case where denominitor
% is zero
% Alternatively, you can use sinc function. doc sinc for details
k1 = zeros(size(a));
idx = ~(a==0);
k1(~idx) = c./g;
k1(idx) = (sin(a(idx)*c))./(a(idx)*g); % ./
plot(n,k1)
