* 장소 : 25동 114호 In classical finite volume methods for the Darcy flow, numerical flux on a face is expressed by numerical pressure values around the face whereas numerical solutions with mixed methods do not have the property. This local dependency of numerical flux on pressure values reflects numerical constitutive law which can be preferred from physical viewpoint. A mixed multipoint flux finite element method (MFMFE) is a variational crime of mixed method which has the local dependency but only the low order method was known. In this work, we propose MFMFE methods for all higher orders in rectangular and hexahedral meshes. A key of the development is construction of a new family of mixed finite elements fulfilling several desired properties.