An improved method for transfer function synthesis from measured frequency response data is presented. The improvement is due to a new scheme of matrix adaptation that eliminates measurement dependent inaccuracies step by step. Arbitrary coefficients can be set to any preassigned values, so that account can be taken of previous knowledge of the investigated system. An extension to integrating systems is developed. The efficiency of the proposed method is demonstrated for an example taken from the literature.