Se desarrolla un enfoque de cuasilinealización por pasos temporales para la forma mixta (o acoplada) de las ecuaciones de Richards de una y dos dimensiones. Para la solución numérica de la ecuación diferencial ordinaria (EDO) lineal para el caso 1D y elíptica para el caso 2D, obtenida después de esta semidiscretización, se utiliza un método de volumen finito para los problemas directos que surgen en cada nivel temporal. A continuación, proponemos una versión del método de descomposición para la solución numérica de la EDO inversa y problemas de contorno elípticos 2D. Se presentan resultados computacionales para algunos tipos de suelos y sus parámetros relacionados reportados en la literatura.