A computer model and numerical method for calculating left epicardial coronary blood flow has been developed. This model employs a finite-branching geometry of the coronary vasculature and the one-dimensional, unsteady equations for flow with friction. The epicardial coronary geometry includes the left main and its bifurcation, the left anterior descending and left circumflex coronary arteries, and a selected number of small branches. Each of the latter terminate in an impedance, whose resistive component is related to intramyocardial compression through a linear dependence on left ventricular pressure. The elastic properties of the epicardial arteries are taken to be non-linear and are prescribed by specifying the local small-disturbance wave speed. The model allows for the incorporation of multiple stenoses as well as aorto-coronary bypasses. Calculations using this model predict pressure and flow waveform development and allow for the systematic investigation of the dependence of coronary flow on various parameters, e.g., peripheral resistance, wall properties, and branching pattern, as well as the presence of stenoses and bypass grafts. Reasonable comparison between calculations and earlier experiments in horses has been obtained.