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