function zdot = phase(s,z)

global n phi

eiphi = exp(i*phi); 
t = 1; f = t;
for k=1:n-1, t = t*z/k; f = f+t; end;
zdot = eiphi*abs(f)/f;

