Numerical resolution of a reinforced random walk model arising in haptotaxis

In this paper we study the numerical resolution of a reinforced random walk model arising in haptotaxis and the stabilization of solutions. The model consists of a system of two differential equations, one parabolic equation with a second order non-linear term (haptotaxis term) coupled to an ODE in a bounded two dimensional domain. We assume radial symmetry of the solutions. The scheme of resolution is based on the application of the characteristics method together with a finite element one. We present some numerical simulations which illustrate some features of the numerical stabilization of solutions.