Coronary flow is different from the flow in other parts of the arterial system because it is influenced by the contraction and relaxation of the heart. To model coronary flow realistically, the compressive force of the heart acting on the coronary vessels needs to be included. In this study, we developed a method that predicts coronary flow and pressure of three-dimensional epicardial coronary arteries by considering models of the heart and arterial system and the interactions between the two models. For each coronary outlet, a lumped parameter coronary vascular bed model was assigned to represent the impedance of the downstream coronary vascular networks absent in the computational domain. The intramyocardial pressure was represented with either the left or right ventricular pressure depending on the location of the coronary arteries. The left and right ventricular pressure were solved from the lumped parameter heart models coupled to a closed loop system comprising a three-dimensional model of the aorta, three-element Windkessel models of the rest of the systemic circulation and the pulmonary circulation, and lumped parameter models for the left and right sides of the heart. The computed coronary flow and pressure and the aortic flow and pressure waveforms were realistic as compared to literature data.