% §3.6 Summation Formulas: The Complex Case
% MATLAB Session p. 71

f=inline('z^(-4)','z');
Z=inline('1/2-i*t/2','t');
dZ = inline('-i/2','t');
F_Summation = inline('real(f(Z(t))*cot(pi*Z(t))*dZ(t)/2/i)','t','f','Z','dZ');
tol = 1e-16; level = 1; even = 'yes'; s = [];
for h=[0.3 0.15]
    s = [s; 
     TrapezoidalSum(F_Summation,h,tol,level,even,f,Z,dZ)];
end
s