function dx = fingeqns(t,x)

% This is the system of differential equations for the OPV model in Eichner
% and Hadeler (see paper for reference).  This function is used by
% 'fingint.m'.
global p m a b c d  

dx = [(1-p)*m-a*x(1).*x(3)-b*x(1).*x(2)-m*x(1);
    p*m+b*x(1).*x(2)-c*x(2)-m*x(2);
    a*x(1).*x(3)-d*x(3)-m*x(3);
    d*x(3)+c*x(2)-m*x(4)];