function ydot = sect_5_3_ode_fn_part_c(t,y)
% SECT_5_3_ODE_FN_PART_C - solve the ODE's from a Chandrasekhar-type system
%
% Written by:
% -- 
% John L. Weatherwax                2005-08-14
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

% some parameters:
omega_n = 10;
zeta = 0.3;
q = 10;

% extract the components:
y  = y(:); 
y1 = y(1);
y2 = y(2);
k11 = y(3);
k12 = y(4);
k22 = y(5);

ydot = zeros(5,1);

ydot(1) = -k11*y1 + (1. - k12)*y2;
ydot(2) = ( -omega_n^2 - k12 )*y1 + (- 2*zeta*omega_n - k22) * y2;
ydot(3) = y1^2;
ydot(4) = y1*y2;
ydot(5) = y2^2;