  function yt=ode1p7(t,y)
%
% Global variables
  global lambda alpha ncall;
%
% ODE
  yt(1)=lambda*exp(-alpha*t)*y(1);
%
% Increment counter for derivative evaluations
  ncall=ncall+1;
