function xp=fitevictplasticitycomp(t,v); global m n c k G ed ug ue ae ag es s yp yh rs rd; %Define v(1) as p and v(2) as h, dominant fitness first and then %subordinate %Differentiate WD wrt p and WS wrt h xp=zeros(2,1); xp(1)= -(.5.*(1+yh)).*k.*(1-tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2).*(G.*(-1+n).*(-1+m)+(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+v(2).*rd.*yp+s.*yh.*yp+G.*v(2)-v(1).*rd+G-s)./(yh.*yp-1)+ed.*(-1+m))./((yh.*yp-1).*(-.5.*n-.5.*n.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))+1.))+(.5+.5.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(.5.*k.*(1-tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2).*n.*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+v(2).*rd.*yp+s.*yh.*yp+G.*v(2)-v(1).*rd+G-s)./(yh.*yp-1).^2-(.5-.5.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*(G.*yh+rd)./(yh.*yp-1))./(-.5.*n-.5.*n.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))+1.)-(.5.*(.5+.5.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)))).*(G.*(-1+n).*(-1+m)+(.5-.5.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+v(2).*rd.*yp+s.*yh.*yp+G.*v(2)-v(1).*rd+G-s)./(yh.*yp-1)+ed.*(-1+m)).*n.*k.*(1-tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2)./((-.5.*n-.5.*n.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))+1.).^2.*(yh.*yp-1))-(.5.*(1+yh)).*k.*(1-tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2).*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+v(2).*rd.*yp+s.*yh.*yp+G.*v(2)-v(1).*rd+G-s)./(yh.*yp-1).^2+(.5-.5.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-m).*(G.*yh+rd)./(yh.*yp-1); xp(2)= (.5.*(yp+1)).*k.*(1-tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2).*(1-ue).*(((.5+.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(G.*(-1+n).*(-1+m)+(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+v(2).*rd.*yp+s.*yh.*yp+G.*v(2)-v(1).*rd+G-s)./(yh.*yp-1)+ed.*(-1+m))./(-.5.*n-.5.*n.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))+1.)-(.5-.5*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+v(2).*rd.*yp+s.*yh.*yp+G.*v(2)-v(1).*rd+G-s)./(yh.*yp-1)).*ae-es)./(yh.*yp-1)+(.5+.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-ue).*((.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*yp.*rd./((yh.*yp-1).*(-.5.*n-.5.*n.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))+1.))-(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-m).*yp.*rd./(yh.*yp-1)).*ae+(.5.*(yp+1)).*k.*(1-tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2).*(1-ug).*(ag.*((.5+.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(G.*(-1+n).*(-1+m)+(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+v(2).*rd.*yp+s.*yh.*yp+G.*v(2)-v(1).*rd+G-s)./(yh.*yp-1)+ed.*(-1+m))./(-.5.*n-.5.*n.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))+1.)-(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k/(yh.*yp-1))).*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+v(2).*rd.*yp+s.*yh.*yp+G.*v(2)-v(1).*rd+G-s)./(yh.*yp-1))+(-v(1).*yh+v(2)).*c./(yh.*yp-1)+(v(2).*yp-v(1)).*rs./(-v(1).*yh+v(2)))./((yh.*yp-1).*(ag.*ug-ag-ug))-(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-ug).*(ag.*((.5+.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*yp.*rd./((yh.*yp-1).*(-.5.*n-.5.*n.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))+1.))-(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-m).*yp.*rd./(yh.*yp-1))+c./(yh.*yp-1)+yp.*rs./(-v(1).*yh+v(2))-(v(2).*yp-v(1)).*rs./(-v(1).*yh+v(2)).^2)./(ag.*ug-ag-ug); %Without rd added to WDE equation %xp(1)= -(.5.*(1+yh)).*k.*(1-tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2).*(G.*(-1+n).*(-1+m)+(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+s.*yh.*yp+G.*v(2)+G-s)./(yh.*yp-1)+ed.*(-1+m))./((yh.*yp-1).*(-.5.*n-.5.*n.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))+1.))+(.5+.5.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(.5.*k.*(1-tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2).*n.*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+s.*yh.*yp+G.*v(2)+G-s)./(yh.*yp-1).^2-(.5-.5.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*G.*yh./(yh.*yp-1))./(-.5.*n-.5.*n.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))+1.)-(.5.*(.5+.5.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)))).*(G.*(-1+n).*(-1+m)+(.5-.5.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+s.*yh.*yp+G.*v(2)+G-s)./(yh.*yp-1)+ed.*(-1+m)).*n.*k.*(1-tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2)./((-.5.*n-.5.*n.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))+1.).^2.*(yh.*yp-1))-(.5.*(1+yh)).*k.*(1-tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2).*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+v(2).*rd.*yp+s.*yh.*yp+G.*v(2)-v(1).*rd+G-s)./(yh.*yp-1).^2+(.5-.5.*tanh((v(2)*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-m).*(G.*yh+rd)./(yh.*yp-1); %xp(2)= (.5.*(yp+1)).*k.*(1-tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2).*(1-ue).*(((.5+.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(G.*(-1+n).*(-1+m)+(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+s.*yh.*yp+G.*v(2)+G-s)./(yh.*yp-1)+ed.*(-1+m))./(-.5.*n-.5.*n.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))+1.)-(.5-.5*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-m).*(-G.*v(1).*yh-G.*yh.*yp +v(2).*rd.*yp+s.*yh.*yp+G.*v(2)-v(1).*rd+G-s)./(yh.*yp-1)).*ae-es)./(yh.*yp-1)-(.5+.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-ue).*(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-m).*yp.*rd.*ae./(yh.*yp-1)+(.5.*(yp+1))*k*(1-tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1)).^2).*(1-ug).*(ag.*((.5+.5*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(G.*(-1+n).*(-1+m)+(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*n.*(1-m).*(-G*v(1).*yh-G.*yh.*yp+s.*yh.*yp+G.*v(2)+G-s)./(yh.*yp-1)+ed.*(-1+m))./(-.5.*n-.5.*n*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k/(yh.*yp-1))+1.)-(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-m).*(-G.*v(1).*yh-G.*yh.*yp+v(2).*rd.*yp+s.*yh.*yp+G.*v(2)-v(1).*rd+G-s)./(yh.*yp-1))+(-v(1).*yh+v(2)).*c./(yh.*yp-1)+(v(2).*yp-v(1)).*rs./(-v(1).*yh+v(2)))./((yh.*yp-1).*(ag.*ug-ag-ug))-(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-ug).*(-ag.*(.5-.5.*tanh((v(2).*yp-v(1).*yh+v(2)-v(1)).*k./(yh.*yp-1))).*(1-m).*yp.*rd./(yh.*yp-1)+c./(yh.*yp-1)+yp.*rs/(-v(1).*yh+v(2))-(v(2).*yp-v(1)).*rs./(-v(1).*yh+v(2)).^2)./(ag.*ug-ag-ug); return;