function dy = dfun(x,y)
%  dy = dfun(x,y)
%  Derivative function for DE: y' = -y
%  
dy = -y;