a b s t r a c tMost real-life problems are computationally burdensome and time consuming, so we have to rely on approximate functions to represent the original functions with a specific level of accuracy. In this study, we propose a new method for approximating multivariate functions based on the diffeomorphic modulation under observable response preserving homotopy (D-MORPH) algorithm. D-MORPH is a regression technique that was originally developed for solving differential equation. Two distinct approaches using the proposed method are described, where one operates over the whole domain and the other operates by sub-dividing the domain into a finite number of sub-domains. This technique is based on a cost/objective function, which depends on a weight matrix. We also introduce a new weight matrix, which dramatically reduces the prediction error and improves the prediction accuracy. The potential of the proposed approach for approximating multivariate functions is illustrated using eight mathematical functions and two practical problems. Furthermore, a comparative assessment with other methods is provided to demonstrate the elegance of the proposed method.