function dydt = rhs1(t,y)
% rhs1   Evaluate right hand side of dy/dt = t - 2*y 
dydt = t - 2*y;
