%
%      ODE_rhs
%

  function dy = ODE_rhs(t,y,A_mat,G_vec);

  dy = A_mat*y + G_vec*sin(3*pi*t);

%