This report describes a new approach to find the global optimal solution for a system composing of M inputs and N outputs. The optimization process consists of defining an objective function in terms of the outputs, sending the value of the objective function to IMIGO, computing the next step change in the inputs with IMIGO, changing the inputs to the values computed by IMIGO, and repeating these steps until convergence is achieved as described here.