It is proposed an engineer oriented method for the optimization of subsonic aerodynamic shape by using the metamodeling technique. The method is based on using the commercial CAD and CFD software as well as the original code EDAOpt developed by authors for design of experiments, approximation and optimization. The method is demonstrated by the solution of a test problem for an airfoil subsonic shape optimization. The developed method essentially decreases necessary computational resources.