The performance modeling study of Delay-Tolerant Network routing, in general, assumes the nodes to be homogeneous (in terms of features such as the coverage range) and uncompromised (in terms of forwarding messages). However, in realistic settings this may not be the case. The routing performance modeling of such realistic scenarios that involve multifariously-featured and egotistic nodes would be interesting and insightful. To this end, in this paper, we analytically model the routing behavior of such nodes using Ordinary Differential Equations for two different routing protocols namely, Epidemic Routing and Two-Hop Routing. Furthermore, we also study the degradation in the routing performance caused by an increase in the fraction of selfish nodes present in the heterogeneous node population. The proposed analytical model is validated via extensive simulations. © 2013 Springer Science+Business Media New York.