Se desarrolló un modelo simple de simulación para predecir la temperatura del suelo bajo cubiertas plásticas sobre la base de la solución numérica de la ecuación de flujo de calor en el suelo. Ella requiere definir las condiciones de borde, que en la superficie son resueltas mediante el balance de radiación y energía. También considera un coeficiente de reparto de calor y el albedo de la superficie del suelo cubierto con plástico. Este parámetro es variable de acuerdo a la hora del día. Así se determinó la temperatura superficial del suelo, suponiendo que era la misma del plástico. Para iniciar el proceso de simulación y determinar la temperatura en el perfil del suelo se recurrió a una función periódica que combina la solución analítica y el análisis de Fourier con un número de elementos k =4. La temperatura a 30 cm de profundidad se consideró invariable para la escala del modelo y estuvo corroborada experimentalmente. Se elaboró un programa computacional que pudiera operar a diferentes lapsos de tiempo y espaciamientos de profundidad, de acuerdo a la difusividad térmica del suelo. La hipótesis de que el plástico está en contacto directo con el suelo permitió que el modelo produjera resultados aceptables al contrastarlo con mediciones reales. Sin embargo, fue forzoso modificar el valor de la emisividad, por lo que se requiere mayor investigación para aclarar este punto.