The slow and fast system approximations are identified by a least squares method. The identifiers are coordinated by a Uzawa algorithm to recover the parameters of the global system. This is the implementation of the method described in this paper. Note: the Matlab Control Toolbox is required to switch from classical representation to zero-pole representation.