function xdot=klfunPo(t,x,flag,cinit,cfinal,tjump)
%%  Reduced Model
%%  Called by F5_6A and F5_6A_Full_Reduced
%%  cfinal is the height of the step jump
%%  tjump is the time of the jump
%%  Carol Lucas


kap=1500; kam=38.4; kbp=1500; kbm=385.9; kcp=1.75; kcm=0.1;
Ka=(kam/kap)^(1/4);Kb=(kbm/kbp)^(1/3);Kc=kcm/kcp;

c = cinit + (cfinal - cinit)*(t>=tjump);

winf=(1+(Ka/c)^4+(c/Kb)^3)/(1+(1/Kc)+(Ka/c)^4+(c/Kb)^3);
tau=winf/kcm;

Po=x(1,1);

w=(1+(Ka/c)^4+(c/Kb)^3)*Po/(1+(c/Kb)^3);
dwdt=(winf-w)/tau;

xdot(1,1)=dwdt*(1+(c/Kb)^3)/(1+(Ka/c)^4+(c/Kb)^3);