function f=AG(x)
%algunas opciones: 
%op = gaoptimset('PopulationSize',100,'InitialPopulation', [-10:1:0]','PlotFcns',{@grafico},'Generations',40);
%llamada al algoritmo:
%[xo,fmin,flag,opciones]=ga(@AG,1,[],[],[],[],-20,20,[],op)

f=(x+10).*(x+6).*(x+5).*(x+1).*(x-7).*(x-10)/1000;