A recently proposed mathematical model for turor angiogenesis consists of a coupled system of ordinary and partial differential equations, some of which are strongly convection dominated diffusion equations. A numerical method based on the use of characterics, which is mass conserving, can be used to effectively handle this feature. The model and the numerical method are presented.