%%  ODE for closed-cell gonadotroph model, Figures 5.12 and 5.13   
%%  Notation implies Li and Rinzel model
%%  Carol Lucas

function xdot=lrclosefun(t,x,flag,I,L,Ve)


fi=0.01;Vi=4;P=26640;
Ki=1.0;Ka=0.4;
Ke=0.2;A=0.5;Kd=0.4;
sigma=0.185;Ct=2;

C=x(1,1);
h=x(2,1);
Ce = (Ct - C)/sigma;

xdot(1,1) = fi/Vi*( (L + P*( (I*C*h)/( (I+Ki)*(C+Ka) ) )^3 )*(Ce - C) - Ve*C*C/(Ke*Ke+C*C) );
xdot(2,1) = A*( Kd - (C + Kd)*h );