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

% mms  $Id: $

udot=1000*(sin(t)-u);