In this paper the progression of tumor angiogenic factor in a tissue is obtained by the numerical solution of a non-linear three dimensional mathematical model using Method of Lines . This method converts a partial differential equation into a system of ordinary differential equations together with the initial and boundary conditions. As the number of lines is increased, the accuracy of the method increases. The method provides very accurate numerical solutions for linear and non-linear problems in comparison with other existing methods. We provide Matlab-generated figures which present the progression of tumor angiogenic factor in porous media and explain the biological importance of this progression. The computer codes are also given.