function cdot = ex76(z, c)
% 
% this file defines Eq.(7.6.10), for both the particular solution and
% the homogeneous solutions;
% c(1) = Cp, c(2) = DCp;
% c(3) = Ch1, c(4) = DCh1;
% c(5) = Ch2, c(6) = DCh2.
%

global Cn I;

% parameters
%
Pe = 10.0;
R = 0.5;

% for particular solution
%
cdot(1) = c(2);
cdot(2) = 2*Pe*R*Cn(I)*c(1) - Pe*c(2) - Pe*R*Cn(I)^2;

% for homogeneous solution 1
%
cdot(3) = c(4);
cdot(4) = 2*Pe*R*Cn(I)*c(3) - Pe*c(4);

% for homogeneous solution 2
%
cdot(5) = c(6);
cdot(6) = 2*Pe*R*Cn(I)*c(5) - Pe*c(6);
cdot = cdot(:);
