function rhs = Fode(t,y)

%
%  save file as Fode.m
%     - right-side forcing for ode45 solver
%

rhs = zeros(size(y));

rhs(1) =  y(2);
rhs(2) = -y(1);
