A precise numerical model based on Crank-Nicolson difference scheme was proposed for calculating water distribution of unsaturated soil under rainfall infiltration.This model has two-order accuracy both in time and space,and by using iterative method in every time step,this model also has higher accuracy.An example of this model was given in this paper to simulate the distribution of water in unsaturated soil under a whole rainfall infiltration course.