function u = f_RK( t, x )

u(1) = -x(2);
u(2) = x(1);

end

Published with MATLAB® R2014a