function dydt = quasiperiodic_deriv ( t, y )

%*****************************************************************************80
%
%% quasiperiodic_deriv() returns the right hand side of quasiperiodic_ode().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    23 July 2020
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    real T, the current time.
%
%    real Y(4,1), the current state values.
%
%  Output:
%
%    real DYDT(4,1), the time derivatives of the current state values.
%
  dydt = zeros ( size ( y ) );

  dydt(1,:) = y(2,:);
  dydt(2,:) = y(3,:);
  dydt(3,:) = y(4,:);
  dydt(4,:) = - ( pi^2 + 1.0 ) * y(3,:) - pi^2 * y(1,:);

  return
end

