function dx = zdv(t,x)
% function dx = zdv(t,x)

dx = zeros(size(x));