  In this short note, we study the gradient estimate of positive solutions to
Poisson equation and the non-homogeneous heat equation in a compact Riemannian
manifold (M^n,g). Our results extend the gradient estimate for positive
harmonic functions and positive solutions to heat equations.
