function udot=ex1_ode(t,u)
% udot=ex1_ode(t,u)
% computes the right side of the ODE du/dt=sin(t)-u
% t,u are scalars
% udot is value of du/dt

% mms  $Id: $

udot=sin(t)-u;