As mentioned earlier, VEGA [5] is the first GA used to approximate the Pareto-optimal set by a set of nondominated solutions. In VEGA, population _ is randomly divided into K equal sized sub-populations; _1 , _2,…, _K. Then, each solution in subpopulation _ is assigned a fitness value based on objective function_. Solutions are selected from these subpopulations using proportional selection for crossover and mutation. Crossover and mutation are performed on the new population in the same way as for a single objective GA.