Abstract Treating specific tissues without affecting other regions is a difficult task. It desirable to target the particular tissue where chemical has its biological effect. To study this phenomenon computationally, in work we numerically mathematical model which written as nonlinear system composed by three parabolic partial differential equations. The variables involved are concentration of ...