The objective of this work is the study of numerical approximation for solving the problem of single phase Darcy flow. The mathematical model is given in three space dimension by one scalar partial differential equation of elliptic nature. Mixed and hybrid finite element method is used for discretization of the model, resulting in an algebraic linear system. We present good approximations of the three-dimensional problem at hand.