function N=NGiros(f,z0,n)
%% 
% Programador: Fredy Vides
% Proyecto: ACRPkG/CICC-UNAH 2020
% Example:
% f=@(x)polyval([.1 .1 .1 .1 .1 .1 .5 1.8 -.7],x);
% N=NGiros(f)
%%
if nargin<2
n=14;
z0=0;
end
if n<=1, n=2;end
t=0:1/n:1;
z=exp(2*pi*i*t(1:n));
pf=polyfit(z,f(z),n);
dpf=polyder(pf);
dpf=@(t)(exp(2*pi*i*t).*polyval(dpf,exp(2*pi*i*t)))./(polyval(pf,exp(2*pi*i*t))-z0);
N=real(quadl(dpf,0,1,1e-14));
