%
% function for fig14-1
%
function dcdt = fun14_1( t, c )
%
global rate
%
dcdt = - rate * c / ( c + 5 );
%
