function [ dh ] = DIFF_H(t,x,y,u)
%
% Written by:
% -- 
% John L. Weatherwax                2006-07-06
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

dh = 1.d0;