Many physical problems need a multidimensional description and involve high dimensional spaces. Standard discretization techniques often lead to an excessive computation time. To solve this problem, we develop in this paper an empirical interpolation decomposition (EID) for multivariate functions. This method provides an approximate representation of a given function in separate form. Error estimates of the developed EID are derived and some properties are given.