%first_rhs.m
function ydot = systems_rhs( t, y )   
ydot = 2*y-t;
return
