An important approach towards understanding the cancer dynamics is the modeling of angiogenesis process. There have been several attempts to model this process. Among them angiogenesis models with time delays, caused by the physical distance between the tumor and the vessel, are the most realistic ones. Recent studies have suggested that those delays can cause oscillatory behavior in the angiogenesis process. In this work we employed piecewise linear hybrid systems with delay on the piecewise constant part. Our approach is based on piecewise linearization of the system behavior where the delays occur at threshold crossings and state transitions. Piecewise linear systems with a single threshold for each variable are useful in approximating and modeling the dynamical systems especially when the model might need to be calibrated by the observations. Therefore, we used piecewise linear systems where the delays are introduced in piecewise constant part of the equations. Our approach allows tractable approximation of the angiogenesis process with possible advances of incorporating more variables, involving the effect of some possible external inputs, and possible adjustment or correction of parameters by observations. (C) 2009 Elsevier Ltd. All rights reserved.