function dx = fingIPVeqns(t,x)

% This is the system of differential equations for the IPV model.  They are
% used by 'fingIPVint.m'

global p m a b c f g h  

dx = [  m*(1-p)-b*(x(3)+g*x(4)).*x(1)-m*x(1);  
        p*m*(1-a)-b*(x(3)+g*x(4))*f.*x(2)-m*x(2);
        b*(x(3)+g*x(4)).*x(1)-(c+m)*x(3);
        b*(x(3)+g*x(4))*f.*x(2)-(c/h+m)*x(4);
        c*(x(3)+x(4)/h)-m*x(5)+m*a*p ];