function v = Dhandle(x)

global B0
global MVcount
global s

v = [x(1);B0*x(2:end)];
if s + 1 > 0
  v = -v; % The handle is -D in this case
end
MVcount = MVcount + 1;