function dpz = dpdz (t)
 T = 10.0;
 w = 2*pi/T;
 dpz = cos(w*t);
 #dpz = -1.;
endfunction