  function yt=ode1p8(t,y)
%
% Global variables
  global a b;
%
% ODEs
  yt(1)=-a*y(1)+b*y(2);
  yt(2)= b*y(1)-a*y(2);
  yt=yt';
