% phase discriminator in Costas loop

alpha = .01;
beta = .0002;

b=[alpha+beta -alpha];
% b=[1];
a=[1 -1];

plotZTP(b,a,50,8000)