SUMMARYAn approximation based on multiple function and gradient information is developed using Hermite interpolation concepts. The goal is to build a high-quality approximation for complex and multidisciplinary design optimization problems employing analysis such as aeroservoelasticity, structural control, probability, etc. The proposed multidimensional approximation utilizes exact analyses data generated during the course of iterative optimization. The approximation possesses the property of reproducing the function and gradient information of known data points. The accuracy of the new approach is compared with linear, reciprocal and other standard approximations. Because the proposed algorithm uses more data points, its efficiency has to be compared in the context of iterative optimization.