y=myfmdemod(x)

function y=myfmdemod(x)
% y = myfmdemod(x)
% fm demodulate the signal

% Author: Eric Spotted Elk, 1 May 2005

p=unwrap(angle(x))/(2*pi);
g=[-0.0012787971436, 0.005141207704, -0.012946369848,...
   0.02740789942, -0.05181317558, 0.09089547374, ...
   -0.15254646724, 0.2535622556, -0.443341752, ...
   0.97056403, 0, -0.97056403, ...
   0.443341752, -0.2535622556, 0.15254646724, ...
   -0.09089547374, 0.05181317558, -0.02740789942, ...
   0.012946369848, -0.005141207704, 0.0012787971436];
y=filter(g,1,p);

%recently added 4/30/2004
%trim off filter settling
y=y(length(g):end);