In this paper, we develop an iterative method of lines scheme for the numerical solution to time fractional Richards equation with implicit Neumann boundary conditions, which is effective tool describing a process flow through unsaturated media. A example provided show effectiveness presented different model parameters and inputs. The illustrated here can be applied other types various input fu...