A model of the cardiovascular system coupling cell, hemodynamics, and autonomic nerve control function is proposed for analyzing heart mechanics. We developed a comprehensive cardiovascular model with multi-physics and multi-scale characteristics that simulates the physiological events from membrane excitation of a cardiac cell to contraction of the human heart and systemic blood circulation and ultimately to autonomic nerve control. A lumped parameter model is used to compute the systemic and pulmonary circulations interacting with the cardiac cell mechanism. For autonomic control of the cardiovascular system, we used the approach suggested by Heldt et al. [2002. Computational modeling of cardiovascular response to orthostatic stress. J. Appl. Physiol. 92, 1239-1254] (Heldt model), including baroreflex and cardiopulmonary reflexes. We assumed sympathetic and parasympathetic pathways for the nerve control system. The cardiac muscle response to these reflex control systems was implemented using the activation-level changes in the L-type calcium channel and sarcoplasmic/endoplasmic reticulum calcium ATPase function based on experimental observations. Using this model, we delineated the cellular mechanism of heart contractility mediated by nerve control function. To verify the integrated method, we simulated a 10% hemorrhage, which involves cardiac cell mechanics, circulatory hemodynamics, and nerve control function. The computed and experimental results were compared. Using this methodology, the state of cardiac contractility, influenced by diverse properties such as the afterload and nerve control systems, is easily assessed in an integrated manner.